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In  this  report  a  single  longitudinal  tracker  is  developed  for  the 
FPCC  aircraft  for  three  different  flight  conditions.  The  method  used  is 
the  singular  perturbation  method  applied  to  fast-sampling, 
output-feedback  digital  control.  Each  flight  condition  has  three  command 
modes:  positive  pitch  pointing,  vertical  translation  and  straight  climb. 
A  sensitivity  study  is  performed  to  validate  the  design  and  illustrate 
design  parameter  influences  on  system  response. 

A  computer-aided-design  program,  MULTI,  is  developed  to  assist  in 
the  iterative  design  process.  The  program  is  fully  interactive, 
user-oriented,  and  provides  error  protection.  The  program  allows 
complete  design  and  simulation  of  three  types  of  control  law  designs: 
known-regular  plants,  known-irregular  plants,  and  unknown  plants.  The 
report  contains  a  brief  but  complete  summary  of  each  of  these  control  law 
design  methods. 

A  user's  manual  and  a  programmer's  manual  are  provided  for 
further  development  of  the  program. 
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Digital  flight  control  systems  (DFCS)  for  aircraft  can  provide 
more  flexibility  and  ability  to  perform  widely  varying  taaka  than  analog 
flight  control  systems  (AFCS).  As  an  example,  the  space  shuttle 
currently  uses  a  digital  fly-by-vire  (FBW)  system  to  control  its  motion 
over  a  large  range  of  speeds  and  altitudes. 

Flight  propulsion  control  coupling  (FFCC)  is  one  of  a  seiies  of 
test  studies  conducted  by  the  Flight  Dynamics  Lab  to  investigate  the 
possibilities  of  integrating  one  or  more  subsystems  on  an  aircraft  to 
provide  greater  performance.  These  subsystems  include  guidance  and 
navigation,  propulsion,  flight  control,  fire  control  and  others.  By 
integrating  these  subsystems  together,  pilot  workload  is  reduced,  mission 
performance  is  optimized  and  survivability  is  enhanced.  The  FPCC  "paper" 
aircraft  used  for  this  thesis  is  a  control-configured  vehicle  (CCV)  and 
can  provide  direct  lift  and  sideforce  as  well  as  the  normal  modes  of 
control.  This  is  accomplished  by  the  use  of  jet  flaps  and  horizontal  and 
vertical  canards.  The  aircraft  is  also  statically  unstable  to  provide 
greater  pitch  response. 

The  final  key  to  the  solution  of  combining  a  DFCS  with  a  CCV 
aircraft  is  a  method  to  integrate  the  controls  to  achieve  the  performance 
desired.  The  method  studied  and  presented  in  this  thesis  is  based  on  the 
latest  work  by  Professor  Brian  Porter  of  the  University  of  Salford, 
England.  He  is  currently  under  contract  to  the  Flight  Dynamics  Lab  and 
his  latest  work  concerns  the  application  of  singular  perturbation  methods 


Co  multivariable  digital  controls  (Refs.  1,  2  and  3).  By  using  output 
feedback,  fast  sampling,  and  a  proportional-plus-integral  controller 
network,  fast  tracking  coupled  with  good  disturbance  rejection  is 
achieved.  Futhermore,  the  controller  is  robust  enough  to  handle  plant, 
input,  and  output  parameter  changes  (such  as  would  occur  at  various 
flight  conditions)  and  still  return  good  performance. 

STATEMEHT  OF  PROBLEM 

The  thesis  is  broken  into  two  separate  and  distinct  tasks: 

(1)  Develop  an  interactive,  user-oriented  computer  program  based 
on  the  algorithms  furnished  by  Professor  Porter  (Ref.  4).  It  includes 
the  three  design  methods  outlined  in  Chapter  II  and  contains  a  flexible 
plotting  package  to  provide  a  quick  graphical  interpretation  of  the 
results.  Finally,  the  program  is  flexible  enough  to  allow  a  meaningful 
iterative  design  process  and  the  ability  to  test  a  single  control  law 
design  over  a  wide  range  of  flight  conditions. 

(2)  Using  the  computer  program,  named  MULTI,  design  a  digital 
flight  control  system  for  the  longitudinal  axis  to  track  three  different 
command  vectors:  positive  pitch  pointing,  vertical  translation,  and 
straight  climb.  The  design  is  based  upon  the  transonic  flight  condition, 
Mach  0.9,  h  *  30,000  ft.  A  robustness  check  is  made  by  applying  the 
newly  developed  control  law  to  two  other  flight  conditions,  subsonic, 
Mach  0.6,  0  ft.,  and  supersonic,  Mach  2.3,  40,000  ft. 

APPROACH 

A  top  down  approach  is  employed  and  modularity  is  stressed  in  the 
development  of  MULTI.  The  list  of  program  requirements  and  programming 
techniques  are  given  in  Chapter  III.  The  core  of  algorithms  supplied  by 
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Professor  Porter  are  adapted  for  use  on  the  CDC  6600  computer  at  AFIT  and 
provide  an  iterative  design  tool. 

The  equations  of  motion  for  the  FPCC  aircraft  are  derived  by  the 
use  of  the  FPCCSIM  program  provided  to  the  Flight  Dynamics  Lab  by 
Honeywell  (Refs.  5  and  6).  Three  flight  conditions  are  examined, 
subsonic,  transonic  and  supersonic,  and  a  reduced  4-state,  3  input,  3 
output  model  is  derived  for  each. 

The  program  MULTI  is  used  to  develop  a  longitudinal  tracker 
control  law  for  the  transonic  flight  condition.  The  control  law  is  then 
applied  to  the  subsonic  and  supersonic  flight  conditions  to  test  the 
robustness  of  the  design.  Finally,  a  sensitivity  study  is  performed  by 
changing  the  individual  design  parameters,  both  to  validate  the  original 
design  and  to  provide  guidance  to  future  designers  as  to  the  influence  of 
the  individual  design  parameters. 

REVIEW  OF  LITERATURE 

The  purpose  of  this  thesis  was  to  obtain  a  working  knowledge  of 
Porter's  latest  work.  The  review  of  literature  consisted  of  the  series 
of  papers  (previously  cited)  dealing  with  singular  perturbation  methods 
in  multi-variable  output  feedback  control.  Other  literature  included 
topics  such  aa  transmission  zeros  (Ref.  7)  and  the  course  material 
presented  in  EE  7.08  -  Linear  Multivariable  Feedback  Systems. 


CHAPTER  II. 


DESIGH  METHODOLOGY  USING  SINGULAR  PETURBATION  METBODS 
IHTRODU  CTI ON 

The  theory  presented  here  is  an  overview  of  Dr.  Brian  Porter's 
papers  (Refs.  1,  2,  3)  presented  during  a  series  of  lectures  at  AF1T  in 
May,  1981. 

Dr.  Porter's  latest  work  concentrates  on  the  application  of 


singular-perturbation  methods  for  multi-variable  control.  Output 
feedback  is  u6ed  to  regulate  the  output  and  track  the  desired  input 
command  vector.  The  state  equations  are: 


x  ■=  Ax  ♦  Bu 

(1) 

y  «  Cx 

(2) 

where : 

x  ■  state  vector 

u  *  input  vector 

y  ■  output  vector 

The 

structure  of  the  controller  is  a 

proportional-plus-integral 

controller . 

Figure  1  shows  it6  structure  in 

the  overall  system  block 

diagram. 
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The  controller  may  either  be  analog  or  digital.  This  thesis 
concentrates  on  the  digital  case.  Samplers  are  assumed  to  be  placed  in 
the  feedback  loop  and  at  the  command  input  vector,  x.  A  sample-and-hold 
device  is  assumed  to  be  placed  after  the  gain  parameter  to  maintain 
piece-wise  constant  inputs  to  the  plant.  The  design  method  of  this 
thesis  does  not  require  a  z-domain  analysis. 

The  controller  matrices,  Kq  and  Kj  are  computed  by  using  the 
first  Markov  parameter  of  the  continuous  plant  as  their  basis.  The 
concept  of  the  first  Markov  parameter  is  important  and  i6  developed  in 
each  of  the  following  sections.  It  is  unique  for  each  of  the  3  types  of 
plants  studied  :  Known  Plants  (regular  or  irregular)  and  Unknown  Plants. 


PIOWH  PLAMTS 


The  category  of  plants  consists  of  those  that  have  a  known 
detailed  state  space  description.  Equations  (1)  and  (2)  are  constrained 
to  the  following  format: 


Xj(t) 

A11 

A12 

XctT 

x2(t) 

_A21 

A22. 

X2Ct ) 

u(t) 


y  <t)  «  icj,  c2i 


Xj(t) 

x,(t) 
.  ^  — 


(3) 

(4) 


where 

n  “  number  of  states 

ft  ”  number  of  inputs  ■  number  of  outputs 
Xj  “  state  vector  partition;  length  (n-I) 
x2  “  atate  vector  partition;  length  (jf) 
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Ajj  ■  A  matrix  partition;  sixe  (n-|)  X  (n-J|) 

Ajj  *  A  matrix  partition,  sixe  (n-|)  X  | 

Ajj  “  A  matrix  partition,  aixe  }  X  (n-J) 

Ajj  “  A  matrix  partition,  size  J>  X  JP 

Bj  “  non-zero  B  matrix  partition;  full  rank;  size  J  X  $ 

Cj  “  C  matrix  partition;  size  I  X  (n-|) 

Cj  “  C  matrix  partition;  size  }  X  ^ 

If  the  6tate  equation  is  not  in  the  form  of  Eq.  (3),  this  form 
can  be  obtained  by  a  transformation.  However,  it  is  not  necessary  to 
perform  this  transformation  to  synthesize  the  control  law.  The  control 
law  is: 

u(kT)  -  f£xoe(kT)  +  XjZtkT)^  (5) 

where 

f  ”  sampling  frequency 
e  “  error  vector  ■>  v  -  y 
z  “  integral  error  vector 
XQ  *  proportional  error  controller  matrix 
Xj  “  integral  error  controller  matrix 
kT  -  sampling  time  interval 
v  “  constant  cownd  vector 

The  closed-loop  system  augmented  with  integrators,  in  order  to 
obtain  the  control  law  of  Eq.  (5),  tracks  a  constant  command  input,  that 
is, 

lim  e(t)  “  0 
t  — ►  «o 

However,  this  is  possible  only  if  the  matrix  pair  (A,B)  is 
controllable  and  the  system  augmented  with  the  integrators  is 
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controllable.  The  latter  condition  is  satisfied  if 


Note,  to  satisfy  the  controllability  requirement,  the  Bj  matrix 
must  have  full  rank.  This  requirement  can  be  satisfied  two  ways.  One 
way  is  to  combine  the  servo  dynamics  into  the  equations.  This  increases 
the  number  of  states  and  complicates  the  design  procedure.  The  second 
way  is  to  combine  two  inputs  that  are  functionally  similar.  This  was 
done  in  this  thesis  by  combining  the  jet  f laps  and  maneuver  flaps,  both 
of  which  provided  large  lift  with  minimal  pitching  moment.  This  approach 
was  taken  since  the  aircraft  studied  has  a  rank-deficient  matrix. 
This  is  discussed  in  Chapter  IV. 

For  most  aircraft,  the  first  partitioned  set  of  states,  x^ , 
usually  consists  of  kinematic  equations  that  have  no  control  inputs,  such 
as  6  -  q,  -  p  and  f* *  r .  The  second  partition,  x2>  consists  of  states 
which  have  control  inputs  through  the  B^  matrix.  Note  that  no  provision 
for  a  feed  forward  (D)  matrix  i6  included.  The  presence  of  a  D  matrix 
(caused  by  choosing  accelerations  a6  outputs)  can  be  circumvented  by 
including  the  servo  outputs  (surface  deflections)  a6  states. 

Once  the  plant  is  in  the  required  format,  its  first  Markov 
parameter  is  determined.  If  C2  is  full  rank,  the  plant  is  described  as 
"Regular"  and  the  matrix  C2  ■  B2  has  full  rank  and  is  the  first  Markov 
parameter.  Note  that  CB  *  C^Bj. 

If  Cj  is  rank  deficient,  the  plant  is  described  as  "Irregular" 
and  C2  must  be  augmented  to  make  it  full  rank.  Design  methodologies  for 
"Regular"  and  "Irregular"  plants  are  described  below. 
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(A)  REGULAR  PLANTS 


If  the  plant  is  regular,  then  the  first  Markov  parameter,  C^B^, 
is  invertible.  Thi6  is  the  basis  for  computing  the  controller  matrices: 


«^(C2B2)"1  £ 

(6) 

(sV1^ 

(7) 

where 


«  proportional  error  control  matrix 
«  integral  error  controller  matrix 
-  proportion  of  direct  and  integral  error  control 
*  first  Markov  parameter 

«  output  weighting  matrix  (selected  as  a  diagonal 
matrix  to  reduce  intercoupling) 


For  regular  plants  the  error  is  the  difference  between  the 
command  input  vector  v(t)  and  the  output  vector  y(t),  that  is, 
e(t)  -  v(t)  -  y(t)  (8) 

Since  the  complete  system  is  augmented  with  a  vector  integrator, 
with  a  constant  command  vector  the  steady-state  value  of  the  error  vector 
is  zero. 

eH(l)  -  0  (9) 

Porter  demonstrates  that  the  closed  loop  transfer-function 
matrix,  GQO»  assumes  the  asymptotic  form, T  (A)>  as  the  perturbation 
parameter,  ^8amp*  approaches  zero  (or  f  )  .  r  (A)  may  be 

partitioned  into  two  components  given  by: 

r  (A) -fo)  ♦Poo  (10) 

where 
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0 


r<A)  -  U~\,0)  lAi*  -  I*  -  tao) 

ra>-c2i>iL  -  iL  ♦  b2koc2]_1b2ko  (i2) 

and 

(13) 

Since  A  is  block  structured,  the  roots  of  Eq .  (11)  can  be 
o 

determined.  They  are  finite  and  are  identified  with  the  "slow  modes". 
They  are  divided  into  two  subsets 


\  0 


lUC2i\\  A11  -  A12C2'1C1 


™12C2 


-1 


(ID 


h  -{A‘C:  |yUL  -  iL  *  <‘*> 

z2  ‘{x*  C :  |AVl  -  h-i  -  «n  *  ™>2c2-‘ci  |  *  <15> 

The  structure  of  Eq.  (11)  shows  that  the  Z^  roots  (associated 

with  the  upper  left  partition  of  A^)  become  unobservable.  Also,  the  Z2 

roots  (associated  with  the  lower  right  partition  of  A  )  become 

o 

uncontrollable. 

The  final  set  of  roots,  Z^,  are  identified  with  the  "fast  mode" 
and  are  determined  as: 

Z3  “^0|AlL  ~  h  *  Wo  |  "  °j  (16> 


As  the  sampling  frequency  (f)  increases,  these  roots  provide  the 

dominant  transient  output  for  the  system.  Evaluation  of  Eq.  (11)  yields 

/■*/ 

the  result  that  P  (X)  “  0.  This  leaves  the  asymptotic  transfer 
function  as 

P  (X)  -  f  (X)  -  CUL  -  \  *  C2B2Ko)  lc2B2Ko  (17) 


This  transfer  function  is  the  result  in  the  limit  as  the  sampling 
time  goes  to  zero.  The  main  thrust  of  this  development  shows  that  as  the 
sampling  time  is  decreased,  only  the  fast  modes  are  present  in  the 
output.  The  matrix  is  chosen  such  that 

C2B2Ko  -  2!“  -  diag  £*-• ,  •  •  •  *1}  (18) 


Then  the  asymptotic  transfer  function  matrix  becomes 


T(X>  -  ♦  2.  > 


-l 


(19) 


..  „  c  _£j__  gj  <n 


Therefore  decoupling  is  accomplished. 

Equation  (19)  can  be  simplified  even  further  by  specifying  ^ 


1^.  Equation  (19)  now  becomes 


■tJ 


(20) 
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Equation  (18)  dictates  the  form  of  the  Kq  matrix: 

‘o-(c2‘2rlE  <21) 

For  the  purposes  of  compatibility  with  the  computer  program, 
MULTI,  and  the  other  design  methods,  and  and  given  in  Eqs.  (6)  and 
(7)  and  are  repeated  below: 


Ko  -  «c(C2B2)  1  2]  (6) 

*1  ‘  <C2"2rl  £  <7> 
limitatiohs 

Transmission  zeros  become  a  problem  if  they  are  located  at  the 
origin  or  in  the  right-half  plane  (Ref.  7).  When  the  state  and  output 
equations  have  the  form  of  Eqs.  (3)  and  (4),  Porter  shove  that  the 
transmission  zeros  are: 


ZT  “  A11  +  A12C2  C1 


(22) 


It  is  noted  that  this  expression  contains  terms  which  are  similar 
to  those  found  in  Eq.  (15).  This  indicates  that  the  second  set  of  finite 
roots,  Z2>  approach  the  transmission  zeros  as  the  sampling  frequency 
approaches  infinity.  If  these  transmission  zeros  are  at  the  origin  or  in 
the  right-half  plane,  the  Z2  roots  may  become  unstable  as  they  approach 
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the  transmission  zeros.  If  this  occurs,  it  may  be  necessary  to  define 
different  output  variables  to  achieve  a  more  desirable  transmission  zero 
location . 


(B)  IRREGULAR  PLAHTS 

Irregular  plants  are  those  whose  first  Markov  parameter  (C2B2^ 
does  not  have  full  rank.  Since  is  constrained  to  have  full  rank  by 
Eq.  (3),  C2  must  be  augmented  to  make  up  its  deficiency.  This  is 
accomplished  with  a  transducer  or  measurement  matrix  M.  The  feedback 
vector  is  now  defined  as  W. 

X.(t) 

w(t)  -  lF.,F_]  (23) 
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It  is  evident  from  Eqa.  (24)  through  (26)  thet  W  ( t  )  ^*> 

y(t)es* 

A  modified  block  diagram  for  the  irregular  plant  deaign  is  shown 
in  Fig.  2. 


Fig.  2  Biock  diagram  -  Irregular  Plant  and  Controller 

The  only  difference  between  Figs.  1  and  2  is  the  extra 

measurements  of  the  states  in  the  feed-back  loop.  The  new  error  vector 
is : 

*  -  v  -  w  (27) 

The  new  first  Markov  parameter  is  now  F2B2.  Tl1*8  requires  F^  to 

have  full  rank  in  order  for  the  first  Markov  parameter  to  be  invertible. 
F2  -  [C2  ♦  MA12]  (28) 

where 

F2  is  J?  x 
C2  ia 

M  ia  X  x  (*.-£) 

*l2  im 
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Some  insight  is  needed  in  order  to  pick  a  measurement  matrix. 


This  is  addressed  later  in  this  section. 


Once  M  is  chosen  properly  and  is  calculated,  the  controller 


matrices  are: 

K 

K1  "  <F2B2)_1  -2T 


-1 


(29) 

(30) 


Note  the  similarities  between  Eqs.  29  and  30,  and  Eqs.  6  and  7, 
respectively.  The  control  law,  Eq.  5,  is  unchanged. 

As  with  the  regular  plant  design,  using  f  *=  as  the 
perturbation  parameter,  the  asymptotic  transfer  function  as  f-*- at)  assumes 
the  form: 


T(A)  -  r q)  *  Pa) 


(31) 


where 


Fa>  -  ICjVW0!  -  c2f2‘1fi>  1  >  -  “o' 


-1. 


.-1 


TA12F2 


-1 


r  a)  -  [c2j  i  a  iL  -  iL .  h2v2i'X 


and 


V1') 

A12F2'V\ 


A11  "  A12F2  F1 


(32) 

(33) 

(34) 


Note  the  similarities  between  Eqs.  32  through  34  and  Eq6.  11 
through  13.  If  *  Cj  in  Eqs.  32  through  34,  the  result  is  Eqs.  11 
through  13,  respectively. 

As  in  the  regular  plant  design,  the  Aq  matrix  is  block  structured 
and  the  roots  can  be  easily  determined.  Equation  32  yields  two  6ets  of 


♦ 
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finite,  closed-loop  eigen  values  associated  with  the  "slow  aodes” 


-.3 


(35) 


Z2  ■  [X‘C-  I  AH-J.  •  *I-L  -  “a  *  “n»2'\|  •  «}  <“> 

Analogous  to  the  regular  plant,  inodes  associated  with  the  Zj 
roots  become  uncontrollable  as  £-*o*  and  effectively  dissappear  as 
sampling  time  T  is  decreased.  On  the  other  hand,  the  mode  associated 
with  the  Z^  roots  remain  observable  due  to  the  extra  measurements 
generated  by  the  M  matrix.  This  is  quite  different  from  the  regular 
plant  design  where  the  modes  associated  with  the  Z^  roots  became 
unobservable  as  T  is  decreased. 

The  third  set  of  roots  are  associated  with  the  "fast  modes" 
derived  from  Eq.  33. 

I3-£X*C:|xh.-H.*  W.|-°3  <37> 

Again,  the  roots  are  the  desired  outputs  since  tracking  becomes 
"tighter"  and  interaction  minimal  as  f  -♦  . 

The  asymptotic  transfer  function  matrix  for  the  irregular  plant 
contains  both  slow  and  fast  modes,  rather  than  just  the  fast  modes  for 
the  regular  plant  seen  in  Eq .  17. 

rft)  -  IC,  -  CjFj-'FjI  [>Vl  -  ^  -  U||  ♦  tT*uV’] 

(SLOW  MODES) 


+ 

IC2F2_1]  ^W"1  [W 

(FAST  MODES) 


(38) 
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It  can  be  shown  that  this  is  a  diagonal  matrix  and  therefore  interaction 
is  reduced  a6  f 00  . 

Equation  37  dictates  the  form  of  the  controller  matrices  by 
choosing  Kq  such  that 


Wo  "  (C2  +  MA12)B2Ko  dUg  -Or]  (39) 

This  leads  to  the  expression  for  K^: 

'o  •  <V2r‘2  <“> 

For  the  purpose  of  compatibility  of  the  computer  program,  MULTI, 


and  the  other  design  methods,  the  matrices  Kq  and 

i,  -  <P2,2)'1zr 


are  re-defined  as 

(29) 

(30) 


MKASmtEKKHT  MATRIX 


The  transmission  zeros  of  the  augmented  system  are: 

ZT  -£«C|Ahl  "  AH  +  A12F2  Fll“°]  (41J 

Note  that  this  equation  contains  similar  terms  ae  Eq.  36,  the 

expression  for  the  second  set  of  finite,  slow  modes,  Z^.  Since 

F2_lyi  “  [C2  +  MBI21~1  IC1  +  **11*  (42) 

It  is  evident  that  any  choice  of  M  influences  the  transmission 

zero  location  for  the  plant.  A  proper  choice  of  M  is  a  rather  easy 
choice  if  M  is  an  L  x  1  vector,  in  other  words,  if  the  number  of  states 
is  one  more  than  the  number  of  outputs.  This  involves  solving  Eq.  36  for 
one  root.  However,  this  becomes  rather  involved  when  Eq.  36  contains  two 
or  more  roots.  Ridgely  (Ref.  8)  has  a  possible  solution  to  this  problem 
by  suggesting  the  use  of  eigenstructure  assignment.  Although  his  work  is 
for  an  analog  system,  the  method  is  adapted  in  this  thesis  for  digital 


16 


system  designs. 


The  structure  roots  are  determined  by 

|  A  V-L  '  >H  -  “ll  ♦  “l  A  \  |  *  *  <«> 

This  is  similar  in  structure  to  the  state  feedback  eigenatructure 
approach  where 


(44) 

A  “  XH-L  +  **11 

(45) 

B  -TA12 

(46) 

l  -  P2_1,l 

(47) 

Equations  45  and  46  can  be  calculated  for  a  specific  sampling 
time,  T«  They  are  entered  into  the  interactive  program,  CESA  (Ref.  9), 
as  the  ASD  and  BSD  matrices,  respectively.  The  CESA  program  requests  the 
number  of  eigen  values  to  be  assigned.  The  program  then  generates  a  null 
space  for  each  eigenvalue  and  the  user  picks  one  or  more  linearly 
independant  eigenvectors  from  those  null  spaces.  The  program  then 
generates  the  control  law  matrix  K  .  Once  K  is  known,  the  measurement 
matrix  can  be  formed  as  follows : 

K  -  [C2  ♦  MA12I_1  [Cj  +  KAjj]  (48) 

M[A12K  -  An]  -  tCj  -  C^]  (49) 

M  -  ICj  -  C^)  [A12K  "  Au)_1  (50) 

There  is  a  solution  to  Eq.  50  if  the  inverse  of  the  matrix  [A^K  -  A^] 
exists . 

The  matrix  manipulations  shown  above  are  included  in  option  18  of 
MULTI,  the  computer  package  developed  as  part  of  this  thesis.  Chapter 
III  and  Appendices  B  and c  contain  a  complete  description  of  the 
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program.  Option  18  returns  an  M  matrix  once  a  K  matrix  from  CESA  is 
entered . 


The  one  draw  back  to  this  method  is  that  the  K  matrix,  and 
consequently  the  M  matrix,  is  not  a  unique  solution  to  the  eigenvalue 
assignment  problem.  K  is  assumed  to  be  an  independent  variable  by  the 
program  CESA.  However,  the  interpretation  of  K  for  the  output  feedback 
methods  of  this  thesis  is  that  it  is  a  function  of  A^  and  A^  which  are 
also  present  in  the  ASD  and  BSD  matrices. 

Porter  used  measurement  matrices  that  include  0.25  as  one  of  the 
non  zero  elements.  This  may  be  a  trial  and  error  figure  that  returned 
suitable  results.  (Refs.  3  and  10).  In  both  reference  papers  the  rank 
deficiency  of  the  matrix  is  eliminated  by  clever  placement  of  non-zero 
elements  in  the  M  matrix. 

A  final  word  of  caution  when  choosing  a  measurement  matrix  is  to 
observe  the  vector  space  of  A^2  and  C^.  Since  the  F2  matrix  is  a  linear 
combination  of  A^2  and  C2  (via  the  M  matrix),  full  rank  may  not  be 
possible  if  A^2  and  C2  share  one  or  more  row  vectors.  If  A^2  and  C2 
share  one  or  more  row  vectors,  the  dimension  of  the  vector  space  from 
which  F2  is  chosen  may  be  less  than  $  .  This  can  occur  if  one  of  the 
output  equations  is  identical  to  one  of  the  equations  and  can  lead  to 
a  non- invertible  first  Markov  parameter. 

PMtMOWH  PLATS 

Porter's  paper  (Ref.  1)  addresses  the  need  for  obtaining  an 
accurate  state-space  representation  of  a  dynamic  plan  and  then 
synthesizing  a  digital  controller  for  it.  Porter's  contention  is  that 
such  a  detailed  investigation  of  the  plant's  dynamic  qualities  is 
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unnecessary  and  that  a  digital  controller  can  be  synthesized  using  the 
plant's  steady-state  transfer  function  matrix,  G(0).  The  matrix  G(0)  is 
easily  obtained  by  offline  tests  and  serves  as  the  fir6t  Markov  parameter 
for  the  plant. 

The  control  law  is  stated  as: 

u  -  TaaMp  £  KQe(kT)  +  (51) 

where 

u  is  the  input  vector 

T _ is  the  sampling  time 

samp 

Kq  is  the  proportional  error  controller  matrix 
Kj  is  the  integral  error  controller  matrix 
e  is  the  error  vector,  v  -  y 

x  is  the  integral  error  vector 

The  only  difference  between  Eq.  51  and  the  control  law  equation 

for  known  plants,  Eq.  5,  is  the  scalar  multiplier  T  (instead  of  f). 

samp 

This  is  the  gain  parameter  6hown  in  Fig.  1.  The  controller  matrices  are 
given  by 

K0  -  ~[G<0)]-1  53  (52) 

Kj  -  1<J(0)]_12Z  (53) 

where 

is  the  proportion  of  direct  to  integral  control 
G(0)  is  the  steady-state  tranfer  function  matrix 
£  is  the  output  weighting  matrix 
Equations  51  through  53,  used  with  the  controller  structure  in 
Fig.  1,  lead  two  sets  of  closed  loop  characteristic  roots,  and 
Z1  "  -  I*  -  1A  +  0(T2)  |  -  0}  (54) 
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(55) 


*2  ■fr*tK  -  h  *  **£  ♦  0<t3)  I  ■  ®3 

where 

Zj  “  set  of  roots  doe  to  plant  states 

Z j  “  set  of  roots  due  to  integrator  states 

I  -  number  of  plant  states 

L  ■  number  of  inputs  “  number  of  outputs 
T  “  sampling  time 


A  “  open-loop  plant  matrix 

0(T2)  -  analytic  expression  whose  value  is  on  the  order  of  T2 


(XT3)  “  analytic  expression  whose  value  is  on  the  order  of  T3 


For  example,  a  plant  whose  open-loop  matrix  contains  only  the 
diagonal  elements  £  -1 ,-l ,-l^  ,  the  Z^  roots  would  be: 

Zj  -  f 1  -  T  ♦  0(T2),  1  -  2T  ♦  0(T2),  1  -  3T  +  0(T2)j  (56) 

If  a  £  matrix  is  chosen  that  contains  only  the  diagonal  elements 
£l20,60^  ,  the  Z2  roots  would  be: 

Z2  -  £l  -  120T2  +  Ott3),  1  -  60T2  ♦  (XT3)}  (57) 

It  is  readily  apparent  that  sampling  time  has  a  dominant  effect 


on  the  closed-loop  roots.  As  T  approaches  zero,  all  roots  converge 
inside  the  unit  circle  to  the  Z  -  1  point. 

While  it  is  possible  to  obtain  an  analytic  solution  for  the 
closed-loop  roots,  the  premise  of  this  design  method  is  that  it  is  better 


to  observe  the  output  time  responses  of  the  system  and  to  adjust  the 
sampling  time  to  obtain  suitable  responses.  This,  in  effect,  bypasses 
the  conventional  approach  of  assigning  closed-loop  eigen  values  and  lets 
the  designer  deal  strictly  with  the  system  output  in  the  time  domain. 


LIMITATIONS  AMP  PARAMETER  CHOICE 

The  unknown  plant  design  method  has  a  few  limitations  of 
particular  interest  to  the  flight  control  engineer.  The  main  limitation 
is  that  the  open-loop  system  must  be  stable,  i.e.  all  roots  of  the 
discrete-time  representation  of  the  system  must  be  inside  the  unit 
circle.  This  is  evident  from  the  following  example.  In  an  unstable 
plant  with  a  single  root  at  •  ■  +3  is  put  into  the  form  of  Eq.  53,  the 
Zj  root  is: 

Zj  -  £l  ♦  3T  -  0(T2)^  (58) 

No  value  of  T  can  bring  this  root  inside  the  unit  circle.  The 
limitation  of  requiring  open-loop  stability  precludes  this  design 
method's  use  for  aircraft  with  relaxed  static  stability  in  the  pitch  axis 
or  spiral  divergence  in  the  yaw  axis. 

The  other  limitation  is  that  this  design  can  not  be  used  if  the 

plant  has  zero-valued  transmission  zeros,  since  this  reduces  the  rank  of 

G(0).  This  leads  to  a  non-invertible  first  Markov  parameter  which  is  the 

basis  of  the  controller  matrices,  K  and  K,  .  This  limitation  occurs 

o  1 

primarily  when  the  outputs  are  angular  rates  such  as  P,  Q,  and  R.  This 
can  be  avoided  by  reforming  the  output  matrix  C  so  that  just  the  angles 
are  the  observed  output. 

While  this  method  is  identified  as  useful  for  an  unknown  plant, 
it  can  also  be  used  when  the  plant  model  is  known.  The  plant  equations 
are  set  up  in  the  normal  matrix  equation  format  using  the  A,  B,  C,  D 
matrices.  Servos  can  be  included  in  the  equations. 

Parameters  affecting  the  response  are  sampling  time  T,  the 
proportion  of  direct  and  integral  feedback  and  the  output  weighting 
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matrix^  •  There  is  an  additional  Bcalar  parameter,  c  that  multiplies 
the  sigma  matrix.  Since  the  design  is  an  iterative  process,  the  most 
logical  route  is  to  set  nominal  values  for  the  parameters  and  to  adjust 
accordingly.  Sampling  time  can  be  nominally  set  at  T  ■  0.01.  Alpha  can 
be  initially  1  for  equal  direct  and  integral  feedback.  Setting  J  *  I 

L 

and  •  ■  y  provides  an  iterative  baseline  from  which  adjustments  can  be 
made.  If  one  output  is  too  slow,  an  increase  in  its  corresponding 
element  in  the^^  matrix  will  correct  this.  If  the  input  energies  and 
outputs  are  all  too  small,  an  increase  in  the  sigma  matrix  multiplier 
corrects  this.  Decreasing  sampling  time  causes  the  outputs  to  track  the 
command  vector  faster.  Finally,  adjustment  of  -*<  affects  the  amount  of 
undershoot  or  overshoot  of  the  outputs. 

COHCLPSIOH 

Three  distinct  design  methods  are  presented:  known  plants  - 
regular,  known  plants  -  irregular,  and  unknown  plants.  Each  has  its  own 
benefits  and  limitations.  Figure  3  represents  a  designer  with  a  flow 
chart  that  outlines  the  design  process. 
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PLANT 

EQUATIONS 
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CHAPTER  III 


"MPLTI"  COMPUTES  PROGRAM  STRUCTURE  AHD  DEVKI-OPMET 

IHTRODOCTIOH 

The  computer  program  "MULTI"  is  an  interactive,  user  oriented 
computer  program  that  develops  digital  controllers  for  error  actuated 
trackers  using  singular  perturbation  methods.  It  is  designed  tc  fulfill 
the  following  objectives. 

1.  The  computer  program  must  be  interactive,  user  oriented,  and 
have  the  ability  to  recover  from  user  input  errors. 

2.  The  program  must  be  flexible  enough  to  allow  for  a  meaningful 
interactive  design  process. 

3.  It  must  be  able  to  test  various  plant  parameter  variations 
(such  hs  different  flight  conditions)  with  a  single  control  law  in  order 
to  demonstrate  the  robustness  of  a  design. 

4.  It  must  contain  a  flexible  plotting  package  in  order  to 
graphically  interpret  design  parameter  changes. 

5.  The  program  must  have  the  ability  to  write  information  to  and 
read  information  from  local  data  files.  This  provides  fast  design 
iterations  and  removes  the  laborious  necessity  of  entering  plant  data 
every  time  the  program  is  accessed. 

6.  It  should  incorporate  all  three  design  methodologies  :  Known 
-  regular  plants,  known  -  irregular  plants,  and  unknown  plants.  It 
should  also  direct  the  user  to  an  alternate  design  procedure  if  the  first 
design  has  limitations. 

7.  It  should  incorporate  such  features  as  inclusion  of  servo  and 


sensor  dynamics  in  the  simulation. 


8.  It  should  also  feature  a  provision  to  include  limits  for  the 
input  values  (such  as  rudder  and  aileron  limits)  in  the  simulation. 

9.  The  program  must  be  written  in  a  top-down  approach  so  that 
future  improvements  can  be  made  with  minimal  problems. 

The  program  "MULTI"  is  based  on  a  core  of  algorithms  written  by 
A.  Hemani,  entitled  "Package  200"  (Ref.  4).  An  analysis  of  this  program 
reveals  that  it  does  not  satisfy  all  of  the  previously  cited  objectives. 
Capt.  Porter  (Ref.  11)  undertook  the  major  effort  of  developing  and 
refining  the  program  MULTI  by  modifying  the  algorithms  in  Hemani's 
program  in  order  to  satisfy  some  of  the  desired  objectives.  The  author 
of  this  thesis  made  a  large  contribution  to  that  effort.  MULTI  is  fully 
interactive,  compatible  with  operation  on  the  CDC  6600's  INTERCOM,  and  is 
easily  accessed.  Lt.  Paschall  provided  the  matrix  inversion  subroutine, 
Capt.  Porter  wrote  the  plotting  package,  and  the  author  provided  the 
routine  to  solve  the  differential  equations  with  the  control  surface 
constraints . 

The  inclusion  of  these  subroutines  into  "Package  200”  required 
too  much  core  memory  for  use  on  INTERCOM.  It  also  lacked  the 
known-regular  and  known-irregular  design  packages,  and  did  not  have  a 
provision  to  save  information  on  data  files.  The  program  also  lacked  the 
inherent  flexibility  required  in  a  user-oriented,  iterative  design 
package. 

The  problem  of  flexibility  is  solved  by  the  use  of  "options"  and 
is  similar  to  the  approach  used  in  the  interactive  program  "TOTAL"  (Ref. 
13).  Plant,  sensor  and  servo  information  is  entered  in  options  1-9. 
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Design  parameter  information  is  entered  and  calculation  of  the  control 
law  matrices  is  performed  in  options  11-19.  Simulation  parameters  are 
entered  and  the  simulation  is  performed  in  options  21-29.  Plotting 
options,  both  CALCOMP  and  terminal  plots,  are  included  in  options  31-39. 
Options  0,10,20  and  30  provide  the  lists  of  input,  design,  simulation  and 
plotting  options,  respectively.  Options  101  to  139  provide  a  print  out 
of  the  data  values  entered  in  options  1  to  39,  respectively. 

QVKKLAY  STRUCTURE 

One  main  problem  in  this  stage  of  the  development  of  MULTI  is  the 
resultant  size  of  the  program.  It  can  not  be  compiled  on  INTERCOM,  which 
has  a  65,000  limit.  In  order  to  adapt  the  program  to  fit  the  computer 
limits  there  are  two  software  structures  to  chose  from:  (1)  segmented 
load  using  loader  commands  to  direct  program  execution  or  (2)  th»  use  of 
overlays.  The  use  of  overlays  is  chosen  as  the  approach  for  this  thesis 
in  order  to  reduce  the  amount  of  storage  required  and  to  make  efficient 
use  of  the  available  field  length.  The  program  is  divided  into  overlays. 
Each  overlay  is  an  executable  program  and  the  overlays  are  a  collection 
of  programs  combined  into  an  overlay  structure.  Since  the  program  is 
divided  into  options  and  a  top-down  approach  is  used,  the  overlay 
structure  is  easily  implemented. 

Overlays  are  classified  as  main  overlays,  primary  overlays  and 
secondary  overlays.  There  is  a  single  main  (or  zero)  overlay  that  is  the 
executive  program  and  directs  the  loading  of  other  overlays.  The  main 
overlay  is  loaded  first  and  remains  in  memory  at  all  times.  A  primary 
overlay  can  be  loaded  following  the  main  overlay,  and  a  secondary  overlay 
can  be  loaded  immediately  following  the  primary  overlay.  This  allows  the 
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program  to  load  into  memory  only  the  primary  overlay  that  ia  required, 
and  a  secondary  overlay  if  necessary.  Since  primary  overlays  can  be  as 
large  or  small  as  necessary,  a  large  amount  of  code  can  be  partitioned 
among  the  primary  overlays  in  order  to  drastically  reduce  the  size  of  the 
resident  memory  core.  The  program  "MULTI"  uses  12  overlays.  These  are 
discussed  in  detail  in  Appendix  B.  A  simplified  block  diagram  is  shown 
in  Fig.  4.  Additional  information  concerning  overlays  can  be  found  in 
the  CDC  FORTRAN  5  manual  (Ref.  14). 

DATA  STORAGE  AND  ERROR  PROTECTION 

Data  is  passed  between  overlays  by  the  use  of  labeled  common 
blocks.  If  an  overlay  uses  any  variable  or  array  that  is  used  in  the 
main  overlay  or  any  other  primary  overlay,  it  must  be  included  in  a 
labeled  common  block  in  the  specific  overlays.  Otherwise,  any  arrays 
that  are  used  or  initialized  are  destroyed  when  another  overlay  is 
loaded.  There  are  24  labeled  common  blocks  in  "MULTI". 

A  useful  feature  of  "MULTI"  is  the  ability  to  open  data  files  and 
write  information  to  them.  Similarly,  the  program  can  also  open 
specified  data  files  to  read  information  from  them.  Overlay  (14,0)  is 
called  when  the  program  is  halted.  This  automatically  writes  the  plant, 
servo  and  sensor  state  equation  information  to  a  file  called  "MEMO".  It 
also  writes  the  design  parameters  and  control  law  matrices  to  a  file 
called  MEM10".  These  local  files  may  be  cataloged  in  the  normal  fashion 
to  save  their  contents. 

In  order  to  read  plant,  servo  and  sensor  state  equation 
information,  the  user  calls  Option  #9  and  types  in  the  name  of  the  local 
file  to  be  accessed.  If  the  user  wishes  to  read  the  G(0)  matrix 
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xlcutive  program 


information  (used  for  unknown  plant  design),  the  user  calls  Option  #8. 
In  order  to  read  in  the  design  parameters  and  control  lav  matrices,  the 
user  calls  Option  #19.  This  gives  the  designer  the  ability  to  enter 
different  flight  condition  for  an  aircraft  and  to  test  them  against  a 
single  control  lav  design. 

Error  protection  is  provided  on  2  levels.  The  first  level  is  the 
ability  of  the  designer  to  correct  a  mistake  when  entering  matrices 
without  having  to  enter  all  the  information  once  again.  The  second  level 
of  error  protection  informs  the  user  if  he  has  bypassed  an  entry  before 
progressing  further.  For  example,  if  all  the  design  parameters  have  not 
been  entered  (Options  11-13),  the  control  laws  are  not  calculated.  In 
this  case  the  program  instructs  the  user  as  to  which  parameter  is 
missing. 

As  a  final  safeguard.  Options  101  through  139  print  out  the 
current  values  of  data  corresponding  to  Options  1  through  39, 
respectively.  This  permits  a  double-check  of  all  parameters. 

With  the  modifications  noted  in  this  Chapter  all  nine  objectives 
listed  previously  are  met. 
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The  aircraft  model  used  in  this  thesis  is  a  "paper"  aircraft.  It 
is  the  result  of  a  joint  effort  between  the  Flight  Dynamics  Lab  and  three 
companies,  Lockheed,  Honeywell  and  Pratt  &  Whitney,  to  explore  the 
effects  of  coupling  the  propulsion  and  flight  control  systems  of  an 
aircraft.  The  resulting  study  (Ref.  15)  and  a  follow-up  study  (Ref.  16) 
demonstrated  the  possiblities  of  flight  propulsion  control  coupling 
(FPCC)  and  its  performance  advantage. 

The  study  set  up  a  methodology  to  choose  between  a  number  of  new 
technologies  in  order  to  satisfy  the  following  guidelines: 

Greatly  expanded  flight  envelope 
Increased  agility 
Improved  stability  and  control 
Improved  handling  qualities 
Reduced  sixe/veight/cost 
Extended  opterational  lifetime 

The  technologies  chosen  to  be  incorporated  into  the  aircraft  are 
presented  in  Table  I  along  with  the  relative  merits  of  a  canard 
configuration  or  an  aft-tail  configuration  for  the  aircraft. 

The  canard  configuration  is  chosen,  primarily  because  the 
aft-tail  configuration  proves  to  be  a  limitation  when  using  jet 
flaps/vectored  thrust. 
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TABLE  I  ADVANCED  AIRCRAFT  TECHNOLOGY  SUITABILITY 


ADVANCED  TECHNOLOGY 

CANARD  CONFIGURATION 

AFT-TAIL  CONFIGURATION 

CANDIDATfe 

RATIONALE 

CANDIDATE 

RATIONALE 

Fly-By-Wire 

(Digital) 

Yes 

RFQ  Reqmt, 

Adv  Control 
Mechanization 

Yes 

RFP  Reqmt 

Adv  Control 
Mechanization 

Relaxed  Inherent 
Static  Stability 

Yes 

Improved 
Controllabil ity 

Yes 

Trim  Drag; 
Reduction 

Maneuver  Flaps 

Yes 

Reduced  Drag, 
Delayed  Buffet 
Onset 

Yes 

Reduced  Drag, 
Delayed 

Buffet  Onset 

Jet  Flap/ 

Vectored  Thrust 

Yes 

Reduced  Drag, 
Delayed  Buffet 
Onset 

No 

Wing-Tail 

Interference, 

Control 

Limitation 

Direct  Side  Force/ 
Lift  Control 

Yes 

Improved  Aiming 
Accuracy,  Ride 
Quality 

Yes 

Improved 

Aiming 

Accuracy, 

Ride  Quality 

Thrust/ 

Drag 

Control 

Thrust 

Reversers 

Tradeoff 

Integration/ 

Redundancy 

With  Jet  Flap 

Tradeoff 

Integration/ 

Redundancy 

With  Jet  Flap 

Modulating 

Dive 

E rakes 

Yes 

Permits  Higher 
Engine  RPM 

Hold 

Yes 

Permits  Higher 
Engine  RPM 

Hold 

CCV 

•  Flutter 

Suppression 

•  Maneuver  Load 
Control 

•  Gust  Load 
Alleviation 

Yes 

Reduced  Weight/ 
Drag  Potential 

Yes 

Reduced 

Weight/Drag 

Potential 

Composite  Primary 
and  Secondary 
Structure 

Yes 

Reduced  Weight 

Yes 

Reduced  Weight 
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The  bulk  of  the  analysis  of  the  use  of  jet  flaps  was  done  by 
Honeywell  with  engine,  inlet  and  aerodynamic  models  provided  by  Lockheed 
and  Pratt  &  Whitney. 

GBHMAT.  MtfiCttIPTKM 

The  aircraft  studied  is  shown  in  Figure  5.  It  is  a  single  seat, 
high-performance,  supersonic  aircraft  designed  for  air  combat  with  an 
air-to-aurf ace  secondary  role.  It  is  powered  by  two  Pratt  &  Whitney 
F-100  engines  which  have  modulating  jet  flap  nozzles  along  the  trailing 
edges  of  the  wings.  A  horizontal  canard  provides  pitch  control  and 
counter  balancing  when  the  jet  flaps  are  deployed. 

The  aircraft  is  statically  unstable  at  M-0.9  and  L-30,000  ft. 

The  flight  control  functions  and  their  respective  limits  are 
indicated  in  Table  II.  In  addition  to  the  conventional  pitch,  roll  and 
yaw  controls,  vertical  cbin  canards  are  available  for  direct  side  force. 
The  function  and  deflection  limits  for  the  surfaces  are  also  shown  in  the 
table.  The  control  deflection  sign  convention  is  chosen  so  that  a 
positive  deflection  tends  to  produce  a  positive  lift,  side  force,  or 
positive  rolling  movement. 
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TABLE  II 


COrrtOL  SPEFACE  FD1CTIMS  ARP  LIMITS 


contOL 

FUCTIOB 

CORROL 

SURFACE 

DEFLECTIOW  LIMITS 

M<1.0  1.0<M 

LIFT 

MAHEUVER  FLAP 
L.E.  FLAP 

JET  FLAP 

CAWABn 

♦15° 

0°  TO  20° 

0°  TO  90° 

(SEE  BELOW) 

LIMITED 

SIDE  FORCE 

CHIU  CANARDS 
LOWER  RUDDER 

+35° 

(SEE  BELOW) 

±15° 

FITCH 

P-AWABT) 

caiabk  FLAP 

+20° 

♦30° 

♦5°  TO  -10° 
±15° 

TAW 

UPPER  RUDDER 
LOWER  RUDDER 

♦30° 

♦30° 

LOCKED 

±15° 

ROLL 

AILEROW  (PER 
SIDE} 

♦20° 

- 

The  horizontal  canard  is  a  variable  camber  surface  having  geared 
leading  edge  and  trailing  edge  flaps.  For  the  purposes  of  this  study, 
they  are  assumed  to  be  fixed.  The  two  vertical  chin  canards  are  slab 
surfaces  which  are  deflected  in  parallel.  The  ailerons  and  rudder  are 
conventional  surfaces. 

FPCCSIM  PROgRAM 

Appendix  A  contains  the  longitudinal  equations  of  motion  for  the 
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aircraft.  The  scope  of  this  study  did  not  include  the  lateral  dynasties. 
As  can  be  seen  from  the  equations  in  the  appendix,  there  are  several 
non-linear  equations  dealing  with  the  supercirculation  effect  and  its 
corresponding  lift  and  drag.  Obviously,  these  equations  do  not  lend 
themselves  to  the  usual  methods  of  determining  the  equations  of  motion. 
In  order  to  analyze  the  aircraft,  the  computer  program  FPCCSIM  (Ref.  5) 
is  used.  FPCCSIM  is  a  full  nonlinear  6  DOF  simulation.  It  transfers  the 
block  diagram  into  state  equations  for  each  of  10  flight  conditions 
listed  in  Table  III.  This  study  uses  the  longitudinal  data  from  flight 
conditions  7,  8,  and  10.  These  correspond  to  a  low  altitude  subsonic 
condition,  high  altitude  transonic  condition  and  a  high  altitude 
supersonic  condition,  respectively.  The  engine  model  used  by  the  program 
had  been  transformed  to  a  state  equation  format  by  Pratt  &  Whitney. 

The  FPCCSIM  program  is  large;  the  source  listing  fills  eight 
boxes  of  punched  cards.  Much  of  the  program  deals  with  the  engine 
simulation.  The  program  was  modified  by  John  E.  Houtz  of  the  Flight 
Dynamics  Lab  so  that  it  would  provide  just  the  state  equations  and  eigen 
values  of  the  aircraft  itself  without  the  flight-control  system  (Ref.  6). 
This  was  necessary  because  the  original  program  would  provide  a 
simulation  complete  with  the  engine  dynamics  and  the  flight  control 
system  inputs  which  would  be  useless  for  the  purpose  of  this  thesis.  The 
revised  program  is  much  shorter  and  less  time-consuming.  The  matrices 
for  the  three  flight  conditions  are  generated  using  a  data  deck  similar 
to  the  one  listed  in  Ref.  4.  Different  flight  conditions  are  analyzed  by 
switching  a  single  parameter  change  card  (PCC)  in  the  deck. 

The  output  of  the  program  is  in  state  equation  format. 
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x  -  Ax  +  Bu  (59) 

y  ■  Cx  ♦  Do  (60) 

For  flight  condition  7,  Mach  0.6,  zero  ft.,  the  state  equations 


are: 
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•=  canard  deflection;  degrees 

■  maneuver  flap  deflection;  degrees 

-  jet  flap  deflection;  degrees 

-  net  thrust,  engine  1,  2;  lb. 

-  inlet  drag  coefficient;  engine  1,  2 
«  angle  of  attack;  radians 

»  total  airspeed;  ft. /sec. 

*  accel  of  C.G  along  x-body  axis;  ft^/seCj 

■  accel  of  C.G  along  z-body  axis;  ft/sec2 
*=  body  axis  pitch  rate;  rad/sec 

“  z-body  axis  velocity;  ft/sec 
«*  x-body  axis  velocity;  ft/sec 

■  pitch  angle;  radians 

-  altitude;  ft. 

■  down  range  position  referenced  to 
initial  body  axis;  ft. 


All  variables  are  perturbation  quantities  from  the  equilibrium 

point . 

A6  discussed  in  Chapter  II,  the  design  methodology  for  known 
plants  for  synthesizing  the  flight  control  system  (Ref.  2)  has  no 
provision  for  a  feed  forward  matrix,  D.  This  can  be  dealt  with  by  either 
transforming  the  equations  into  a  new  6tate  space  that  includes  actuator 
dynamics,  or  by  reconfiguring  the  output  matrix  to  include  terms  that  do 
not  contain  feed-forward  terms.  The  latter  route  is  chosen  after  the 
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former  method  required  a  large  measurement  matrix  that  was  3X6.  There 
is  little  guidance  in  choosing  the  elements  of  a  large  measurement 
matrix.  The  two  engine  thrusts  are  combined  into  one  input,  F  total. 
The  CDI  terms  are  dropped  on  the  assumption  that  the  inlet  drag  does  not 
change  significantly  at  each  flight  condition.  The  jet  flaps  and 
maneuver  flaps  are  geared  together  so  that  full  deflection  occurs 
simultaneously  -15°  for  maneuver  flaps  and  75°  for  the  jet  flaps.  This 
is  done  because  the  matrix,  with  four  inputs,  does  not  have  full  rank. 
These  two  surfaces  provided  the  same  input,  large  lift  with  small 
negative  pitching  moment,  and  were  therefore  redundant.  The  output 
matrix  consists  of  flight  path  angle,  pitch  angle  and  forward  velocity. 
The  resulting  simplified  6  state,  3  input,  3  output  euations  for  flight 
condition  #7  are: 
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■=  flight  path  angle;  radians 

«=  lift  surface  deflection,  geared  cosibination  of 
jet  flap  and  aaneuver  flap;  degrees 
F  total  =  total  net  thrust;  lb. 
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The  same  approach  is  used  for  flight  condition  #9,  high  altitude 
transonic  cruise,  Mach.  0.9,  30,000  ft.  The  equations  are: 
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The  equations  for  flight  condition  #10,  high  altitude  super8onic 
cruise,  Mach  2.3,  40,000  ft.,  are: 
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CHAPTER  ? 


The  longitudinal  tracker  developed  for  the  FPCC  aircraft  iB  baaed 
on  a  4-state  model.  This  model  reduction  results  from  an  inability  to 
choose  a  measurment  matrix  for  the  6-state  models  developed  in  Chapter 
IV.  The  original  design  is  accomplished  using  flight  condition  #9  -  Mach 
0.9,  30,000  ft.  (which  is  considered  to  be)  the  most  critical  of  the 
three  flight  conditions  under  study.  Thus,  by  applying  the  control  law 
for  this  transonic  flight  condition  to  the  subsonic  (flight  condition  #7) 
and  the  supersonic  (flight  condition  #10)  conditions,  a  rigorous  test  of 
Porter's  claim  of  robustness  can  be  made.  His  claim  that  fast  sampling, 
error-actuated  control  is  robust  enough  to  handle  wide  parameter  changes 
in  the  plant  is  supported  by  the  results  obtained  in  this  chapter.  In 
hindsight,  however,  it  appears  that  the  most  critical  flight  condition 
may  have  been  the  subsonic  condition,  based  on  the  velocity  tracking 
results . 

Three  command  vectors  are  chosen:  positive  pitch  pointing, 
vertical  translation,  and  straight  climb.  The  positive  pitch  pointing 
command  vector  is 

VT  -  10.1,0,0]T  (69) 

This  is  a  0.1  radian  positive  pitch  angle  with  no  change  in  flight  path 
angle  and  no  change  in  velocty. 

The  vertical  translation  command  vector  consisted  of: 

▼T  -  [0,0.1,0]T  (70) 
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This  represente  a  0.1  radian  positive  flight  path  angle  with  no  change  in 
pitch  angle  nor  velocity. 

The  straight  climb  command  vector  consisted  of: 

▼T  -  10.1,0.1 ,0]T  (71) 

This  represents  a  combination  of  a  0.1  radian  pitch  angle  coupled  with  a 
0.1  radian  flight  path  angle.  This  effectively  regulates  angle  of  attack 
to  aero.  Velocity  change  is  regulated  to  zero. 

TEAK SONIC  FLIGHT  CONDITION 

The  four  state  plant  and  output  equations  for  the  transonic 
flight  condition,  Mach  0.9,  30,000  ft.  is  as  follows: 
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The  servos  are  chosen  as  first-order  lags,  based  on  the  design 
specifications  from  Lockheed  for  the  FPCC  aircraft.  The  engine  was 
modeled  as  a  first-order  lag  based  on  the  Pratt  &  Whitney  test  data  (Ref. 
16).  They  are  as  follows: 

cammed 


S+100 

The  design  parameters  are  chosen  as: 

SAMPLING  TIME  -  .01  sec 

ALPHA  -  2.5 

SIGMA  MATRIX  -  [1..5..03] 

DIAGONAL  ELEMENTS 

EPSILOM  -  0.1 

MEASUREMENT  -  [.2,0,0]T 

MATRIX 

These  design  parameters  are  the  result  of  a  lengthy,  iterative 
design  study.  The  rationale  behind  the  choice  of  these  parameters  is 
included  in  Chapter  VI.  These  parameters  result  in  an  excellent 
combination  of  input  and  output  responses  for  all  three  command  vectors 
with  respect  to  input  and  output  figures  of  merit.  Figures  of  merit  for 
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the  output**  include  overshoot,  rise  time  and  settling  time.  The  figures 
of  merit  for  the  inputs  were  the  maximum  deflections  and  number  of 
oscillations  (#  of  peaks). 

The  controller  matrices  formed  from  the  design  parameters  are: 


*T 


2.281 

-0.706 

-391.2 

0.9125 

-0.2825 


-156.5 


6.604  0. 554X1 O"5 
53.45  -0.218X10"3 
1208  3.433 

2.642  0. 221X1 O-3 
21.38  -O.843X10-4 
483  1.373 


(74) 


(75) 


The  input  and  output  responses  are  shown  in  Fig.  6a  to  Fig.  6e 

for  the  positive  pitch  pointing  command.  Theta  (pitch  angle)  and  gamma 

(flight  path  angle)  are  combined  on  a  single  graph  since  these  two 

outputs  are  of  prime  importance.  Theta  reaches  a  steady  state  value  of 

0.1  radian  while  gamma  approaches  zero.  Velocity  perturbation  has  a 

minimum  of  -2.5  ft/sec  and  a  maximum  of  1.3  ft/sec.  This  perturbation  is 

from  a  nominal  velocity  of  830  ft/sec,  and  is  hardly  noticeable.  The 

steady-state  values  for  the  inputs  are: 

G4MABD:  -14° 

LIFT  FLAP:  -70° 
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Figure  6a.  Canard  Deflection,  Positive  Pitch  Pointing  Command  (Transonic  Flight  Condition) 


POSITIVE  PITCH  POINTING  COMMAND  (TRANSONIC  FLIGHT  CONDITION) 


POSITIVE  PITCH  POINTING  COMMAND  ( TRANSONIC  FLIGHT  CONDITION) 


POSITIVE  PITCH  POINTING  COMMAND  (TRANSONIC  FLIGHT  CONDITION) 


The  engine  is  still  oscillatory  at  T»10  sec. 

As  seen  in  the  lift  flap  response,  a  negative  steady  state  angle 
for  the  lift  flap  is  needed  to  achieve  positive  pitch  pointing.  This  is 
needed  to  counteract  the  additional  lift  generated  by  maintaining  an 
abnormal  angle  of  attack.  The  resultant  positive  pitch  moment  is 
countered  by  downloading  the  canard.  This  is  a  feasible  method  to 
achieve  positive  pointing  except  that  the  jet  flap  (geared  together  with 
the  maneuver  flap  to  form  a  lift  flap)  can  not  have  a  negative 
deflection.  If  the  point-of-viev  is  taken  that  this  aircraft  is  still  in 
the  design  stage,  these  results  would  dictate  the  need  for  a  jet  flap 
capable  of  positive  and  negative  deflections. 

The  vertical  translation  command  vector  results  in  the  inputs  and 
outputs  shown  in  Fig.  7a  through  Fig.  7e.  In  this  mode,  gamma  is 
commanded  to  0.1  radian  with  theta  regulated  to  zero.  Velocity 
perturbation  has  a  minimum  of  -2.4  ft/sec  and  a  maximum  of  1.4  ft/sec. 
The  steady  state  values  of  the  inputs  are: 

CAHAKD:  14° 

LIFT  FLAP:  70° 

The  engine  is  still  oscillatory  at  T*=10  secs. 

.  The  straight  climb  vector  results  in  the  input  an  outputs  shown 
in  Fig.  8a  through  Fig.  8e.  Both  theta  and  gamma  track  to  0.1  radian 
velocity  perturbation  has  a  minimum  of  -4.1  ft/sec  and  a  maximum  of  2 
ft/sec.  The  steady  values  for  the  input  are: 

CAHAKD:  0° 
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Figure  7a.  Canard  Deflection,  Vertical  Translation  Command  (Transonic  Flight  Condition) 
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VERTICAL  TRANSLATION  COMMAND  (TRANSONIC  FLIGHT  CONDITION) 
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STRAIGHT  CLIMB  COMMAND  (TRANSONIC  FLIGHT  CONDITION) 
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STRAIGHT  CLIMB  C 


LIFT  FLAP:  0° 

BHGIHK:  1300  lb. 

All  pertinent  figures  of  merit  for  the  three  command  vectors  are 
presented  in  Table  IV  for  this  flight  condition. 

SUBSONIC  FLIGHT  CONDITION 

The  four  state  plant  and  output  equations  for  flight  condition  #7 
-  Mach  0.6,  zero  ft.  are  as  follows: 
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TABLE  IV  FFCC  CONTROL  LAW-THREE  FUCSTT  CONDITIONS 


COWAND 

VECTOR 

INPUT  OR 
OUTPUT 

FIGURES  OF 
MERIT 

TRANSONIC 

MH).9,hr=3Qk 

SUBSONIC 

M=O.6,h=0 

POSITIVE 

PITCH 

POINTING 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

17% 

5.9  sec 

1  pk 

8% 

5  sec 

1  pk 

GAM4A 

INTERACTION 
#  OF  PEAKS 

12% 

1  pk 

17% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.5  to  1.3 

-13  to  7.8 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

-3*  to  -16° 

3  pks 

-7°  to  -17° 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

-12*  to  -75* 

1  pk 

-7* to  -117° 

1  pk 

ftotal 

MAX  VALUES 
#  OF  PEAKS 

1400  to  -750 

3  pks 

-1800  to  6900 

2  pks 

VERTICAL 

TRANS¬ 

LATION 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

3% 

3.9  sec 

2  pks 

0 

5.8  sec 

2  pks 

THETA 

INTERACTION 
#  OF  PEAKS 

-20% 

1  pk 

+22%  to  -14% 

4  pks 

VELOCITY 

MAX  VALUES 

-2.4  to  1.4 

9.2  to  -5.5 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

39* to  10° 

2  pks 

21  to  6* 

5  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

300*  to  15* 

3  pks 

300*  to  3 f 

3  pks 

ftotal 

MAX  VALUES 
#  OF  PEAKS 

6700  to  -1400 
3  pks 

4300  to  -1070 
3  pks 

STRAIGHT 

CLIMB 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

0 

0.7  sec 

1  pk 

1% 

1  sec 

3  pks 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

9% 

1.7  sec 

2  pks 

10% 

0.8  sec 

3  pks 

VELOCITY 

MAX  VALUES 

-4.1  to  2 

-4.4  to  2.5 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

36*to  -9° 

2  pks 

14"  to  -10* 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

298*1X3  -25° 

4  pks 

O  o 

302  to  -38 

2  pks 

ftotal 

MAX  VALUES 
#  CF  PEAKS 

7000  to  -570 

4  pks 

10040  to  -2900 
4  pks 

SUPERSONIC 

M=2.3,h=40k 

16% 

5.6  sec 
1  pk 

8% 

1  pk 

-54  to  37 

-t  to  -15* 

4  pks 

-9*  to  -50° 

1  pk 

-600  to  ;34000 

2  pks 

12% 

2.6  sec 

2  pks 

-39% 

3  pks 

51  to  -36 

64* to  -11° 

2  pks 

251°  to  -42* 

4  pks 

-31800  to  124 

2  pks 

0 

0.7  sec 

1  pk 

17% 

0.5  sec 

3  pks 

18  to  -17 

49°  to  -13° 

2  pks 

239*' to  -68* 

2  pks 

-7000  to  7600 
2  pk 


The  servos  end  actuators  are  the  same  as  stated  previously. 

In  order  to  check  the  robustness  of  the  transonic  design,  the 

same  controller  matrices,  K  and  K.  ,  and  the  same  measurement  matrix,  M, 

o  1 

are  used  for  the  subsonic  flight  condition. 

The  positive  pitch  pointing  command  vector  produces  the  input  and 

output  results  shown  in  Figs.  9a  through  Fig.  9e.  The  theta  overshoot  iB 

less  than  the  transonic  condition,  but  the  gamma  interaction  is  slightly 

higher.  The  velocity  perturbation  is  larger  with  a  -13  ft/sec  minimum 

and  a  7.8  ft/sec  maximum.  This  is  the  perturbation  from  a  nominal 

velocity  of  670  ft/sec  and  may  or  may  not  be  readily  apparent  to  the 

pilot.  The  steady  state  values  for  the  inputs  are: 

CHUBS:  -17° 

LIFT  FLAP:  -117° 

The  engine  is  still  oscillating  at  T*10  sec.  The  steady  state 
value  for  the  lift  flap  exceeds  the  physical  maximum  of  75°  which 
indicates  that  the  command  vector  is  too  large  or  that  the  lift  flap  does 
not  have  the  "authority"  to  perform  the  needed  task.  This  should  be 
taken  into  account  if  the  FFCC  aircraft  evolves  any  further  than  its 
current  status  as  a  "paper”  aircraft. 

The  inputs  and  outputs  for  the  vertical  translation  command 
vector  are  shown  in  Fig.  10a  through  Fig.  lOe.  Both  the  gamma  overshoot 
and  the  theta  interactior  are  less  than  the  for  transonic  condition  but 
the  velocity  perturbation  is  greater.  The  maximum  perturbation  is  9.2 
ft/sec  and  the  minimum  is  -5.5  ft/sec.  Although  this  is  less  than  a  22 
deviation  from  the  nominal  value  of  670  ft/sec,  it  may  be  noticeable. 

The  steady  state  values  for  the  inputs  are: 
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Figure  9a.  Canard  Deflection,  Positive  Pitch  Pointing  Command  (Subsonic  Flight  Condition) 
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POSITIVE  PITCH 
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POSITIVE  PITCH  POINTING  COMMAND  (SUBSONIC  FLIGHT  CONDITION) 


Figure  10a.  Canard  Deflection,  Vertical  Translation  Command  (subsonic  Flight  Condition) 
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VERT  I C  RL  TRANSLATION  COMMAND  (SUBSONIC  FLIGHT  CONDITION) 
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CA1AED:  17° 

LIFT  FLAP:  117° 

The  engine  is  still  oscillatory  at  T«10  seconds.  The  steady  state  value 

for  the  lift  flap  again  indicates  that  either  the  command  vector  is  too 

large  or  the  lift  flap  lacks  the  "authority"  to  perform  a  vertical 

translation  with  a  0.1  radian  flight  path  angle. 

The  straight  climb  command  vector  result  are  shown  in  FigB.  11a 

through  lie.  Theta  and  gamma  track  a  0.1  radian  command  with  slightly 

more  oscillatory  behavior  than  the  transonic  condition  results  of  Figs. 

6a  through  Figs.  8e.  The  velocity  perturbation  is  also  quite  similar  to 

the  transonic  results ;  the  maximum  perturbation  is  2.5  ft/sec,  the 

minimum  is  -4.4  ft/sec.  These  would  hardly  be  noticed  by  the  pilot. 

The  steady  state  values  for  the  inputs  are: 

CAHAED:  0° 

LIFT  FLAP:  0° 

The  engine  is  still  oscillatory  at  T-10  sec. 

The  pertinent  figures  of  merit  for  all  three  command  vectors  are 
presented  in  Table  IV  under  the  subsonic  flight  condition  heading. 
gmnasorrc  flight  common 

Tfcr  four  state  plant  and  output  equations  for  the  supersonic 
flight  condition  #10  -  Mach.  2.3,  40,000  ft  -  are  as  follows: 


o 

o 

<D 
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Figure  11a.  Canard  Deflection,  Straight  Climb  Command  (Subsonic  Flight  Condition) 
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_ TIME,  SECONDS _ 

STRAIGHT  CLIMB  COMMAND  (SUBSONIC  FLIGHT  CONDITION) 
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_ HUE.  SECONDS _ 

COMMAND  (SUBSONIC  FLIGHT  CONDITION) 
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0 

0.497x10 

1.0 
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The  servos  and  sensors  are  the  same  as  stated  previously. 

This  flight  condition  is  the  final  check  for  robustness.  The 
controller  matrices,  Kq  and  and  the  measurement  matrix  are  the  same 

ones  developed  for  the  transonic  flight  condition. 

The  positive  pitch  pointing  results  are  shown  in  Figs.  12a 
through  12e.  The  theta  overshoot  and  gamma  interaction  are  reduced  from 
the  values  for  the  transonic  flight  condition.  Velocity  perturbation  has 
a  minimum  of  -54  ft/sec  and  a  maximum  of  37  ft/sec  from  the  nominal 
velocity  of  2230  ft/sec.  This  corresponds  to  less  than  a  2  1/2Z 
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POSITIVE  PITCH  POINTING 
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POSITIVE  PITCH  POINTING  COMMAND! SUPERSONIC  FLIGHT  CONDITION) 


POSITIVE  PITCH 


variation  in  apeed  which  may  be  detected  by  the  pilot. 

Steady  atate  input  valuea  are: 

CAIAU):  >10° 

LIFT  FLAP:  -45° 

The  engine  ia  still  oscillatory  at  T*10  sec.  The  maximum 
variations  from  the  nominal  thrust  ia  +34,000  and  -600  lb.  Since  the 
maximum  augmented  thrust  from  both  engines  is  31,500  lb.,  it  is  apparent 
that  this  command  vector  is  too  large  for  the  aircraft  to  maintain 
nominal  velocity. 

The  vertical  translation  command  vector  results  are  shown  in 
Figs.  13a  through  13e.  The  gamma  overshoot  is  12Z  but  the  theta 
interaction  peaks  at  -0.039,  a  391  negative  interaction.  This  is  the 
only  glaring  problem  in  the  robustness  test.  Velocity  perturbation  has  a 
maximum  of  51  ft/sec  and  a  minimum  of  -36  ft/sec  from  the  nominal 
velocity  of  2230  ft/sec.  This  is  slightly  more  than  a  2Z  variation  in 
velocity  which  may  be  noticeable. 

The  input  steady  state  values  are: 

CAHA1D:  10° 

LIFT  FLAP:  45° 

The  engine  is  still  oscillatory  at  T“10  sec.  The  total  thrust  varies 
from  -31,800  lb  to  +124  lb  from  the  nominal  thrust.  This  again  is 
clearly  outside  the  normal  operation  of  the  engine.  The  command  vector 
is  too  large  for  the  aircraft  at  this  supersonic  condition. 

The  straight  climb  vector  results  are  shovn  in  Figs.  14a  through 
14e.  The  overshoot  is  17Z,  larger  than  the  transonic  case.  Theta  still 
exhibits  no  overshoot.  Velocity  perturbation  is  quite  small;  maximum 
perturbation  is  18ft/sec,  minimum  is  -17  ft/sec.  This  is  less  than  a  II 
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Figure  13a.  Canard  Deflection,  Vertical  Translation  Command  (Supersonic  Flight  Condition) 
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VERTICAL  TRANSLATION  COMMAND  (SUPERSONIC  FLIGHT  CONDITION) 
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variation  from  the  nominal  value  and  would  be  imperceptable  to  the  pilot. 


The  steady  state  values  for  the  inputs  are: 

CAMA1D:  0° 

Lin  FLAP:  0° 

The  engine  is  still  oscillatory  at  T«10  secs.  The  thrust  perturbation; 
-7,000  lb  to  +7,600  lb  from  the  nominal,  is  well  within  the  limits  of  the 
engine.  This  would  indicate  that  a  0.1  radian  straight  climb  is  quite 
feasible  at  this  flight  condition. 

COMCLOSIOM 

The  results  of  this  chapter  confirm  Porter's  claim  of  robustness 
of  the  design  method.  The  only  limitation  on  performance  is  the  physical 
limits  of  the  aircraft.  The  use  of  a  single  controller  could  prove  to  be 
a  radical  departure  from  the  current  practice  of  gain- scheduling  for  the 
controller  elements  based  on  flight  conditions. 

A  few  words  are  in  order  concerning  the  size  of  the  transients  of 
the  inputs,  particularly  in  the  vertical  translation  and  straight  climb 
modes.  The  lift  flap  initial  deflections  may  reach  300°  in  some 
instances,  obviously  out  of  the  normal  operating  range.  This  is  caused 
primarily  by  the  command  vector  step  inputs.  These  transients  could  be 
reduced  or  even  eliminated  by  shaping  orramping  the  inputs  to  the  final 
value.  This  could  be  a  final  design  parameter  to  tailor  the  input 
response  to  an  acceptable  form  while  maintaining  output  fidelity. 
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CHAPTER  TI 


■<4  4\*V. 


or  htoividhal  DESIOI  pm«mg 


lETRODPCTIOH 

The  design  of  the  longitudinal  tracker  control  lav  for  the 
transonic  flight  condition  presented  in  Chapter  V  is  the  result  of  a 
lengthy  iterative  design  procedure.  This  chapter  is  vritten  to  assist 
any  designer  who  wishes  to  develop  a  digital  control  lav  using  the 
singular  perturbation  methods  discussed  in  Chapter  II.  It  is  assumed 
that  the  designer  has  determined  vhich  type  of  plant  and  design  he  will 
use.  It  is  further  assumed  that  a  first-cut  design  may  have 
objectionable  characteristics,  such  as  excessive  overshoot,  and  that  the 
designer  needs  guidance  in  selectiting  the  various  design  parameters. 
The  parameter  changes  presented  here  are  compared  to  the  base  line 
figures  of  merit  of  the  design  presented  in  Chapter  V.  Since  the  design 
is  based  on  an  irregular  plant,  a  designer  UBing  a  regular  or  an  unknown 
plant  design  does  not  need  the  measurement  matrix  information  section. 
The  conclusions  concerning  the  other  parameters  is  the  same  for  all  three 
design  methods,  except  as  noted. 

One  of  the  biggest  factors  in  the  design  procedure  is  the 
presence  of  the  servoB.  The  servo  dynamics  are  not  included  in  the 
design  due  to  the  resulting  complexity  and  difficulty  of  a  proper  choice 
of  a  measurement  matrix.  By  excluding  the  servo  dynamics  in  the  design, 
the  response  is  made  slow  enough  bo  that  the  servo  effects  can  not  be 
seen  (in  the  form  of  rapid  transients)  in  the  inputs  and  outputs.  This 
is  quite  difficult  with  the  engine  model  "r.ervo"  whose  time  constant  is 
1.77  sec. 


101 


HUBS  or  Dirmxrr  samplimg  times 

In  all  three  design  methods,  e  smaller  sampling  time  gives 
tighter  control  -  less  overshoot  and  shorter  settling  time.  Three  sample 
times  are  examined: 

BASELUE  (T-O.Ol) 

SMALL  (TO.005) 

LARGE  (T-0.05) 

In  this  design,  however,  smaller  sample  times  resulted  in  an  increase  in 
the  frequency  and  magnitude  of  the  transients.  This  is  unacceptable  in 
terms  of  engine  fatigue,  since  the  number  of  peak  thrust  transients 
correlate  to  the  decrease  in  engine  life.  The  large  sampling  time  of 
0.05  sec  causes  the  aircraft  to  go  unstable  for  all  three  command 
vectors.  Since  the  aircraft  is  unstable  in  the  pitch  axis,  this  sampling 
time  appears  to  be  too  large  to  counter-act  the  instability.  The  figures 
of  merit  for  all  3  sampling  times  and  for  all  3  command  vectors  are 
presented  in  Table  V.  Figure  15  is  a  plot  of  pitch  and  flight  path 
angles  for  a  positive  pitch  pointing  command  vector,  T  ■  0.005  sec.  Fig. 
16  is  the  same  plot  but  for  T  “0.05  sec.  These  can  be  compared  to  the 
baseline  in  Fig.  6d. 

EFFECT  OF  DIFFERENT  MKARITBKKKWT  MATRICES 

In  the  irregular  design,  a  proper  choice  of  a  measurement  matrix 
yields  satisfactory  roots  as  discussed  in  Chapter  II.  For  the  4 
state,  3  input,  3  output  model  developed  for  the  transonic  flight 
condition,  a  measurement  matrix  of  the  form: 

M  -  [Mj  ,0,0]*  (80) 

yields  a  single  Z^  root: 


\ 
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TABLE  V  EFFECTS  CF  DIFFERENT  SAILING  TIMES 


CttMJWD 

INPUT  OR 

FIGURES  OF 

VECTOR 

OUTPUT 

MERIT 

POSITIVE 

THETA 

OVERSHOOT 

PITCH 

SETTLING  TIME 

POINTING 

#  CF  PEAKS 

GAtf-lA 

INTERACTION 
#  OF  PEAKS 

VELOCITY 

MAX  VALJJiS 

CANARD 

MAX  VALUES 

CF  PEAKS 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

FTCTAL 

MAX  VALUL'S 
#  OF  PEAKS 

’VERTICAL 

GAM>JA 

OVERSHOOT 

TRANS¬ 

SETTLING  TIME 

LATION 

#  OF  PEAKS 

THETA 

INTERACTION 
#  CF  PEAKS 

VELOCITY 

MAX  VALUES 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

IJFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

F 

MAX  VALUES 

TOTAL 

#  OF  PEAKS 

STRAIGHT 

THETA 

OVERSHOOT 

CLJMH 

SETTLING  TIME 
#  CF  PEALS 

GAWi 

OVERSHOOT 
SETH  ENG  TIME 
#  OF  PEAKS 

VELOCITY 

MAX  VALUES 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

ftotal 

MAX  VALUES 
#  CF  PEAKS 

BASELINE 

17% 

8% 

Unstable 

5.9sec 

4.5  sec 

Unstable 

1  pk 

1  pk 

Unstable 

12% 

6% 

Unstable 

1  pk 

1  pk 

Unstable 

-2.5  to  1.3 

-1.5  to  0.7 

Unstable 

-3° to  -16° 

-3° to  -28° 

Unstable 

3  pks 

6  pks 

Unstable 

1" to  -75° 

-2. 5" to  -73° 

Unstable 

1  pk 

3  pks 

Unstable 

1400  to  -750 

4080  to  -1900 

Unstable 

3  pks 

7  pks 

Unstable 

3% 

9% 

Unstable 

3.9  sec 

2.4  sec 

Unstable 

2  pks 

2  pks 

Unstable 

-20% 

-9% 

Unstable 

1  pk 

1  pk 

Unstable 

-2.4  to  1.4 

-2.3  to  1.7 

Unstable 

39* to  10" 

26° to  -6.7° 

Unstable 

2  pks 

4  pks 

Unstable 

300° to  15° 

180°  to  -105° 

Unstable 

2  pks 

4  pks 

Unstable 

6700  to  -1400 

5000  to  -1500 

Unstable 

3  pks 

6  pks 

Unstable 

0 

0 

Unstable 

Q.7sec 

0.7  sec 

Unstable 

1  pk 

1  pk 

Unstable 

9% 

11% 

Unstable 

1.7  sec 

0.5  sec 

Unstable 

2  pks 

3  pks 

Unstable 

-4.1  to  2 

-3  to  1.5 

Unstable 

36' to  -9° 

-101*  to  100° 

Unstable 

2  pks 

3  pks 

Unstable 

298' to  -25" 

12 4" to  -168° 

Unstable 

4  pks 

4  pks 

Unstable 

7000  to  -570 

5400  to  -3500 

Unstable 

4  pks 

6  pks 

Unstable 
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Figure  15.  Pitch  and  Flight  Path  Angles,  Positive  Pitch  Pointing  Command  (T=0.005) 
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Figure  16.  Pitch  and  Flight  Path  Angles,  Positive  Pitch  Pointing  Command  (T=0.05) 


(81) 


The  baseline  measurement  matrix  is: 

M  -  [0. 2,0,0]”* 

for  T  *0.01  sec,  Z 2  ■  0.95. 

Two  other  measurement  matrices  were  examined.  Measurement  matrix 


A  is: 

Ma  -  10.25.0,0]* 

This  yields  Z 2  ■  0.96  (for  Tsamp  ■  0.1) 

Measurement  matrix  B  is: 

M,  -  [0.15,0,0]* 

which  yields  Z^  “  0.933. 

This  tends  to  indicate  that  M  will  gives  a  slower  response  than 

A 

the  baseline.  This  is  confirmed  in  the  results  of  Table  VI.  Measurement 
matrix  A  produces  a  slower  response  with  larger  overshoots  and  settling 
times  than  the  baseline  response.  Measurement  matrix  B  gives  marginally 
better  responses  than  the  baseline  but  at  the  expense  of  increased 
transients . 

When  choosing  a  measurement  matrix,  the  designer  may  wish  to 
place  his  Z 2  roots  farther  into  the  unit  circle  than  the  values  used 
here.  This  can  lead  to  stability  problems  in  some  cases.  Due  to  the 
reciprocal  nature  of  Eq .  81,  a  small  value  for  M^  places  the  Z2  root 
close  to  the  origin  in  the  unit  circle.  This  may  be  desirable,  but  the 
problem  remains  that  the  extra  measurements  generated  by  the  measurement 
matrix  may  also  be  too  small  to  be  significant.  Thus,  the  designer  must 
optimize  his  choice  of  measurement  matrix  based  on  the  "region  of 
activity"  in  the  Z-plane,  usually  close  to  the  Z*1  point  for  small 
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TABLE  VI  EFFECTS  OF  DIFFERENT  MEASUREMENT  MATRICES 


COMAND 

VECTOR 

INPUT  OR 
OUTPUT 

FIGURES  OF 
MERIT 

BASELINE 

M=(.2,0,0)T 

ma 

M=(.25,0,0)T 

M=(.15,0,0)T 

POSITIVE 

PITCH 

POINTING 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

17? 

5.9  sec 

1  pk 

- 77i - 

6.3  sec 

1  pk 

- 351 - 

5.4  sec 

1  pk 

GAM4A 

INTERACTION 
#  CF  PEAKS 

12% 

1  pk 

12% 

1  pk 

15% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.5  to  1.3 

-2.6  to  1.4 

-2.5  to  1.2 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

-3*  to  - 16  ° 

3  pks 

-2°  to  -15* 

3  pks 

-5°to  -25° 

5  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

l'to  -75° 

1  pk 

-f  to  -80* 

1  pk 

— l*to  -80~ 

3  pks 

FTOTAL 

MAX  VALUES 
#  OF  PEAKS 

1400  to  -750 

3  pks 

980  to  -560 

2  pks 

2380  to  -1100 
5  pks 

VERTICAL 

TRANS¬ 

LATION 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

3% 

3.9  sec 

2  pks 

3% 

3.8  sec 

2  pks 

3% 

4.0  sec 

2  pks 

• 

THETA 

INTERACTION 
#  OF  PEAKS 

-20% 

1  pk 

-25% 

1  pk 

-18% 

4  pks 

VELOCITY 

MAX  VALUES 

-2.4  to  1.4 

-2.3  to  1.4 

-2.5  to  1.5 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

39*  to  10T 

2  pks 

39*  to  9° 

2  pks 

40°  to  10* 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

300 "to  15° 

2  pks 

300° to  15° 

2  pks 

300* to  15° 

2  pks 

F 

TOTAL 

MAX  VALUES 
#  OF  PEAKS 

6700  to  -1400 
3  pks 

6700  to  -100 

3  pks 

6700  to  -210 

4  pks 

STRAIGHT 

CLIMB 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

0 

0.7  sec 

1  pk 

0 

0.9  sec 

1  pk 

0 

0.5  sec 

2  pks 

GAMBIA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

9% 

1.7  sec 

2  pks 

8% 

2.1  sec 

2  pks 

11% 

1.4  sec 

1  pk 

VELOCITY 

MAX  VALUES 

-4.1  to  2 

-4  to  1.9 

-4.1  to  2.0 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

36*to  -9° 

2  pks 

O  O 

37  to  -3 

2  pks 

34* to  -13w 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

298* to  -25* 

4  pks 

297’to  -15° 

2  pks 

300° to  -43° 

2  pks 

FTOTAL 

MAX  VALUES 
#  CF  PEAKS 

7000  to  -570 

4  pks 

6900  to  160 

4  pks 

7400  to  -570 

4  pks 

107 


sampling  times.  Fig.  17  shows  the  pitch  and  flight  path  angles  for  a 

vertical  translation  command  vector  and  measurement  matrix  M.  .  Fig.  18 

A 

is  the  same  response  but  for  measurement  matrix  M  .  These  can  be 

D 

compared  with  Fig.  7d  for  the  baseline  response.  These  figures 
illustrate  the  effects  of  different  measurement  matrices. 

EFFECTS  OF  DIFFEREHT  AT.PHA  PARAMETERS 

Alpha  is  the  proportion  of  direct  to  integral  feedback.  In  the 
continuous  domain,  SISO  case,  alpha  dictates  the  position  of  the  zero  in 
the  lag-lead  network  formed  by  the  proportional-plus-integral  controller. 

In  the  digital  domain,  alpha  has  a  primary  influence  on  settling  time 
with  a  secondary  influence  on  overshoots.  As  ex  i8  increased,  integral 
action  is  diminished  and  settling  time  increases.  This  is  coupled  with 
slightly  lower  overshoots  in  some  cases.  The  reverse  is  also  true.  As  or 
is  decreased,  settling  time  decreases  but  the  lag  contributed  by 
increased  integral  action  contributes  to  larger  overshootB. 

The  choice  of  alpha  in  this  design  must  be  constrained.  If«*c 
is  chosen  less  than  2,  the  velocity  grows  as  an  exponential  sinusoid. 
This  may  be  due  to  the  lag  introduced  by  the  engine  model.  Choosing 
alpha  as  2.5  maintains  stability  and  still  achieves  adequate  performance. 
Two  test  cases  are  run,  one  with  or*  1  and  the  second  with  or*  4.  The 
results  are  presented  in  Table  VII. 

For  or*  1,  the  overshoots  are  much  larger  than  for  the  baseline 
results  (  <x  *  2.5).  The  settling  times  are  smaller  and  the  velocity 
perturbation  is  a  growing  exponential  sinusoid. 

For  of*  4,  the  settling  times  are  larger  and  most  overshoots  are 
smaller  than  the  baseline  results.  Thus,  the  choice  of  «*f  *2.5  seems  to 
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Figure  18.  Pitch  and  Flight  Path  Angles,  Vertical  Translation  Command  (  M= [0.15, 0,0] 


TABLE  VII  EFFECTS  CF  ALPHA  PARAMETER  CHANGES 


COMAND 

VECTOR 

INPUT  OR 
OUTPUT 

FIGURES  OF 
MERIT 

BASELINE 

-^=2.5 

LOWER 

^"=1.0 

HI®ER 
^  =4.0 

POSITIVE 

PITCH 

POINTING 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

17% 

5.9  sec 

1  pk 

68% 

3.2  sec 

1  pk 

11% 

7.8  sec 

1 

GAMMA 

INTERACTION 
#  CF  PEAKS 

12% 

1  pk 

36% 

1  pk 

8% 

6  pks 

VELOCITY 

MAX  VALUES 

-2.5  to  1.3 

Unstable 

-2  to  0.3 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

-3* to  -16* 

3  pks 

9. 6* to  -19* 

1  pk 

18*to  -10* 

1  pk 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

1°  to  -75* 

1  pk 

-2 ’to  -102° 

1  pk 

0.3* to  -74* 

1  pk 

FTOTAL 

MAX  VALUES 
#  OF  PEAKS 

1400  to  -750 

3  pks 

1260  to  -1180 

3  pks 

2800  to  -670 

5  pks 

VERTICAL 

TRANS¬ 

LATION 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

3% 

3.9  sec 

2  pks 

1  % 

1.7  sec 

1  pk 

15% 

5  sec 

2  pks 

£ 

t 

* 

THETA 

INTERACTION 
#  CF  PEAKS 

-20% 

1  pk 

-58% 

1  pk 

-13% 

1  pk 

V 

t 

t 

f 

VELOCITY 

MAX  VALUES 

-2.4  to  1.4 

Unstable 

-2.6  to  1.1 

i 

* 

* 

c 

j  \ 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

39'  to  10* 

2  pks 

24°  to  14* 

2  pks 

37° to  -Id* 

4  pks 

1  i  l. 

'  '  Y 

’ 

I  .TFT  FLAP 

MAX  VALUES 
#  CF  PEAKS 

300 ’to  15° 

2  pks 

140"  to  70° 

2  pks 

252* to  -15/ 

3  pks 

f 

i 

F 

TOTAL 

MAX  VALUES 
#  OF  PEAKS 

6700  tc  -1400 

3  pks 

4700  to  300 

4  pks 

5700  to  -3400 
7  pks 

►  * 

STRAIGHT 

CLLMB 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

0 

0.7  sec 

1  pk 

12% 

1.4  sec 

2  pks 

0 

0.7  sec 

1  pk 

l 

GAMMA 

CVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

9  % 

1.7  sec 

2  pks 

15% 

3  sec 

1  pk 

21% 

0.7  sec 

3  pks 

'  '  ' 

:  ! 

VELOCITY 

MAX  VALUES 

-4.1  to  2 

Unstable 

-3.4  to  0.6 

*  * 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

36°  to  -9° 

2  pks 

28*to  0' 

2  pks 

18' to  -20° 

4  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

29 8' to  -25" 

4  pks 

188'to  -0.5** 

2  pks 

255° to  -204° 

3  pks 

ftotal 

MAX  VALUES 
#  CF  PEAKS 

7000  to  -570 

4  pks 

3450  to  -150 

2  pks 

8500  to  -3900 
5  pks 
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be  the  best  trade  off  between  settling  times,  overshoots  and  stability 
for  the  system.  Figs.  19  and  20  show  flight  path  and  pitch  angles,  for  a 
positive  pitch  pointing  command  vector  (for  or  ■  1  and  o<  -  4).  These 
can  be  compared  to  the  baseline  performance  seen  in  Fig.  6d. 

EFFECT  OF  DIFFERENT  KPRTI-OW  PABAOTCTOfi 

Epsilon  is  a  scalar  parameter  that  multiplies  each  element  in  the 
output  weighting  marix.  By  increasing  a,  the  entire  2  weighting  matrix 
is  increased,  leading  to  quicker  and  tighter  control:  faster  setting 
times,  lower  overshoots.  If  e  is  decreased,  the  control  responses  are 
characterized  by  slower  settling  times  and  larger  overshoots. 

Two  values  of  epsilon  are  used,  a  larger  value,  (e  ■  0.2),  and  a 
smaller  value,  (c  -  0.05),  than  the  baseline  value,  (e  -  0.1).  The 
results  of  the  two  trials  are  compared  to  the  baseline  results  in  Table 
VIII.  For  c  *  0.2,  the  result  is  tighter  control  in  the  form  of  lower 
overshoots,  interaction  and  settling  times.  This  is  coupled  with  the 
previously  mentioned  problems  associated  with  tighter  control;  an 
increase  in  the  frequency  and  magnitude  of  initial  transients.  For 
epsilon  -  0.05,  looser  control  is  seen  in  the  form  of  large  overshoots 
and  larger  settling  times.  These  results  are  seen  in  Figs.  21  and  22. 
Fig.  21  shows  flight  path  and  pitch  angles,  for  a  positive  pitch  pointing 
command  vector  and  c  “  0.2.  Fig.  22  is  the  same,  but  for  c  ■  0.05. 
These  figures  can  be  compared  to  the  baseline  plot  of  Fig.  6d. 

EFFECTS  OF  SIGMA.  MATBTX  PJjjEWT  CHAHGES 

The  sigma  matrix  determines  the  position  of  the  closed  loop  fast 
roots  (Z^)  discussed  in  Chapter  II.  Each  element  on  the  diagonal  has  a 
direct  influence  on  the  response  of  its  corresponding  output.  In  a  fast 
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Figure  20.  Pitch  and  Flight  rath  Angles,  Positive  Pitch  Pointing  Coimands  (<x=4) 


TABLE  VIII  EFFECTS  OF  EPSILON  PARAMETER  CHANGES 


CGM4AND 

VECTOR 

INPUT  OR 
OUTPUT 

FIGURES  OF 
MERIT 

BASELINE 
*  =0.1 

LARGER 
£  =0.2 

SMALLER 

£■=0.05 

POSITIVE 

PITCH 

POINTING 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

17% 

5.9  sec 

1  pk 

8% 

4.5  sec 

1  pk 

49% 

7.1  sec 

1  pk 

GAtMA 

INTERACTION 
#  OF  PEAKS 

12% 

1  pk 

6.3% 

1  pk 

30% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.5  to  1.3 

-1.5  to  0.7 

-4.1  to  1.7 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

-3* to  -16* 

3  pks 

-2' to  -30° 

5  pks 

4° to  -17° 

1  pk 

LEFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

1°  to  -75° 

1  pk 

l°to  -75° 

5  pks 

-1.8° to  -88“ 
1  pk 

FTOTAL 

MAX  VALUES 
#  OF  PEAKS 

1400  to  -750 

3  pks 

4400  to  -2300 

5  pks 

900  to  -840 

3  pks 

VERTICAL 

TRANS- 

1ATICN 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

3% 

3.9  sec 

2  pks 

11% 

2.4  sec 

5  pks 

0 

4.9  sec 

2  pks 

THETA 

INTERACTION 
#  OF  PEAKS 

-20% 

1  pk 

-10% 

3  pks 

-48% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.4  to  1.4 

-2.5  to  1.7 

-2.2  to  0.7 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

39*to  10° 

2  pks 

34*to  -21° 

6  pks 

27° to  13“ 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

300* to  15° 

2  pks 

253'to  -209° 

5  pks 

212° to  70° 

2  pks 

F 

TOTAL 

MAX  VALUES 
#  OF  PEAKS 

6700  to  -1400 

3  pks 

6700  to  -3900 

5  pks 

4760  to  810 

2  pks 

STRAIGHT 

CLIMB 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

0 

0.7  sec 

1  pk 

0 

0.7  sec 

1  pk 

0 

1  sec 

1  pk 

GAMMA 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

9% 

1.7  sec 

2  pks 

18% 

0.7  sec 

5  pks 

5% 

3.6  sec 

1  pk 

VELOCITY 

MAX  VALUES 

-4.1  to  2 

-3.1  to  1.5 

-5.1  to  2.4 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

36*  to  -9° 

2  pks 

16*  to  -20* 

5  pks 

31  to  0° 

2  pks 

LEFT  FLAP 

MAX  VALUES 
#  CF  PEAKS 

298“to  -25° 

4  pks 

196*to  -212 

6  pks 

210*  to  -.02“ 
2  pks 

ftotal 

MAX  VALUES 
#  CF  PEAKS 

7000  to  -570 

4  pks 

7900  to  -6200 

5  pks 

4025  to  530 

2  pks 
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sampling  environment  there  is  s  minimum  effect  on  the  other  outputs  due 
to  the  decoupling  st  high  sample  rates. 


If  a  single  element  in  the  £  matrix  is  changed,  the 
corresponding  output  responds  in  a  similar  manner  discussed  previously  in 
the  epsilon  parameter  change  section.  If  the  element,  dp,  is  increased, 
tighter  control  and  faster  response  occurs.  If  the  element  is  decreased, 
a  slower  response  will  occur.  Changing  an  individual  element  also 
changes  the  input  that  has  the  most  direct  influence  on  the  specific 
output.  This  is  illustrated  by  changing  the  third  element,  o^ ,  in  the 
matrix.  Two  values  are  used  to  determine  the  effects  on  the  system  - 
higher  value,  (o^  -  0.06),  and  a  lower  value,  ■  0.015),  than  the 

baseline  value,  (<y^  *  0.03).  The  results  are  presented  in  Table  IX. 

It  is  readily  apparent  that  the  element  affects  only  its 
corresponding  output  (velocity).  The  only  input  affected  is  the  thrust. 
This  indicates  a  high  degree  of  decoupling  for  the  syBtem.  For  ■ 
0.06,  the  velocity  perturbation  is  smaller  in  all  cases,  coupled  with  an 
increase  in  initial  transients  in  the  Ftotal  input.  This  agrees  with  the 
previous  statements  concerning  tighter  control  and  fast  response.  For 
<3^  ■  0.015,  the  velocity  perturbations  are  larger  and  thrust  transients 
are  smaller.  This  also  agrees  with  the  previous  conclusions  concerning 
looser  control  and  slower  response.  These  results  are  seen  graphically 
in  Figs.  23  and  24.  Fig.  23  is  the  velocity  response  for  a  straight 
climb  vector  and  ■  0.06.  Figure  24  is  the  same  response  but  for  0^  - 
0.015.  These  can  be  compared  to  the  baseline  in  Fig.  8e. 

C0MCUJ8I0M 

The  design  process  to  develop  a  digital  control  design  using 


TABLE  IX  EFFECTS  CF  SICMA  MATRIX  ELEMJ7T  CHANGES 


COMAND 

VECTOR 

POSITIVE 

PITCH 

POINTING 

INPUT  OR 

OUTPUT 

THETA 

FIGURES  OF 
MERIT 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

BASELINE 
ay=o.  03 

17% 

5.9  sec 

1  pk 

HIGHER 

<3V=0.06 

18% 

6  sec 

1  pk 

LOWER 

<J?f=0.015 

18% 

6  sec 

1  pk 

GAIWA 

INTERACTION 
#  CF  PEAKS 

12% 

1  pk 

12% 

1  pk 

12% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.5  to  1.3 

-1.6  to  0.8 

-3.8  to  1.1 

CANARD 

MAX  VALUES 
#  OF  PEAKS 

-3° to  -16° 

3  pks 

-3°to  -16° 

3  jks 

-3 "to  -16° 

3  pks 

LIFT  FLAP 

MAX  VALUES 
#  OF  PEAKS 

l*to  -75° 

1  pk 

— 1*  to  -77* 

2  pks 

-l°to  -77° 

2  pks 

FTOTAL 

MAX  VALUES 
#  OF  PEAKS 

1400  to  -750 

3  pks 

1400  to  -590 

6  pks 

1500  to  -800 
5  pks 

VERTICAL 

TRANS- 

IATICN 

GAM4A 

OVERSHOOT 
SETTLING  TIME 
#  OF  PEAKS 

3% 

3.9  sec 

2  pks 

3% 

3.9  sec 

2  pks 

3% 

3.9  sec 

2  pks 

THETA 

INTERACTION 
#  OF  PEAKS 

-20% 

1  pk 

-20% 

1  pk 

-20% 

1  pk 

VELOCITY 

MAX  VALUES 

-2.4  to  1.4 

-2.3  to  1.6 

-2.8  to  1.4 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

39* to  10° 

2  pks 

39*  to  10* 

2  pks 

39* to  10* 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  CF  PEAKS 

300® to  15* 

2  pks 

300*  to  14° 

2  pks 

30(f  to  15° 

2  pks 

ftctal 

MAX  VALUES 
#  CF  PEAKS 

6700  to  -1400 
3  pks 

7000  to  160 

4  pks 

6500  to  -540 
4  pks 

STRAIGHT 

CLIMB 

THETA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

0 

0.7  sec 

1  pk 

0 

0.7  sec 

1  pk 

0 

0.7  sec 

1  pk 

GAf-WA 

OVERSHOOT 
SETTLING  TIME 
#  CF  PEAKS 

9% 

1.7  sec 

2  pks 

9% 

1.8  sec 

2  pks 

9% 

1.8  sec 

2  pks 

VELOCITY 

MAX  VALUES 

-4.1  to  2 

-3.1  to  1.6 

-5.7  to  2.4 

CANARD 

MAX  VALUES 
#  CF  PEAKS 

36  *to  -9a 

2  pks 

36* to  -9' 

2  pks 

-36° to  -8° 

2  pks 

LIFT  FLAP 

MAX  VALUES 
#  CF  PEAKS 

298* to  -25° 

4  pks 

298° to  -27° 

2  pks 

298^ to  -27° 

2  pks 

ftotal 

MAX  VALUES 
#  CF  PEAKS 

7000  to  -570 

4  pks 

7300  to  300 

4  pks 

6900  to  -525 
4  pks 
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Figure  23.  Velocity  Perturbation,  Straight  Climb  Command  (  O' =0.06) 
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Figure  24.  Velocity  Perutrbation ,  Straight  Climb  Command  (  <TL  =0.015) 


singular  perturbation  methods  is  inherently  iterative.  The  designer 
needs  a  grasp  of  the  influence  of  parameter  changes  on  the  design  with 
respect  to  figures  of  merit  for  outputs  and  limitations  on  the  inputs. 
Table  X  provides  the  designer  with  this  information.  Experience  has 
shown  that  the  designer  must  look  at  all  the  required  commands  to  the 
system,  such  as  pitch  pointing,  vertical  translation,  and  climb,  rather 
than  just  optimizing  around  a  single  input  vector. 
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TABI£  X  SUMMARY  CF  PARAMETER  CHANGE  INFLUENCE 


PARAMETER 

CHANCE 

OVERSHOOT 

SETTLING  TIME 

SAMPLING  TIME 

INCREASE 

INCREASE (P) 

INCREASE  (P) 

DECREASE 

DECREASE  (P) 

DECREASE (P) 

MEASUREMENT 

MA33RIX 

INCREASE 

ELEMENT 

INCREASE (P) 

INCREASE  (S) 

DECREASE 

ELEMENT 

DECREASE (P) 

DECREASE (S) 

AIPHA 

INCREASE 

DECREASE (S) 

INCREASE  (P) 

DECREASE 

INCREASE (S) 

DECREASE  (P) 

EPSILON 

INCREASE 

DECREASE (P) 

DECREASE (P) 

DECREASE 

INCREASE  (P) 

INCREASE  (P) 

SIGMA 

MATRIX 

INCREASE 

HI 

DECREASE  (P-I) 

DECREASE (P-I) 

DECREASE 

ELEMENT 

INCREASE (P-I) 

INCREASE  (P-I) 

(P)=  PRIMARY  INFLUENCE 

(S)=  SECONDARY  INFLUENCE 

(P-I)  =  PRIMARY  INFLUENCE-INDWIDUAL  OUTPUT 

TRANSIENTS 
DECREASE  (P) 
INCREASE  (P) 

DECREASE  (P) 
INCREASE (P) 

INCREASE  (P) 
DECREASE (P) 

INCREASE (P) 
DECREASE  (P) 

INCREASE (P-I) 
DECREASE (P-I) 
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CHAPTER  VII 


COECLPSIOES  AID  tECOMMEEDATIOHS 


COMCLOSIOES 

An  interactive  computer  program,  MULTI,  was  written  and  used  to 
develop  and  test  a  longitudinal  tracker  for  the  FPCC  aircraft.  Chapter 
II  provides  a  brief  summary  of  three  papers  presented  by  Porter  that  form 
the  basis  for  the  design  methods  (Ref.  1  to  3).  A  flow  chart  is  given  at 
the  end  of  that  chapter  to  assist  the  reader  in  choosing  the  proper 
design  method  based  upon  the  constraints  imposed  by  each  of  the  three 
design  methods:  known-regular  plants,  known-irregular  plants  and  unknown 
plants . 

The  computer  program  MULTI,  an  interactive  user-oriented  computer 
aided  design  tool,  was  developed  to  satisfy  the  requirements  listed  in 
Chapter  III.  The  chapter  also  contains  a  brief  outline  of  the  program's 
structure  and  flow.  A  much  more  detailed  description  iB  found  in 
Appendix  B.  Appendix  C  contains  a  user's  manual  that  assists  the  MULTI 
user  in  applying  the  design  methods.  These  two  appendices  are  from  Capt. 
Porter's  thesis  and  are  included  for  completeness  (Ref.  11).  The  program 
developement  was  a  joint  effort  between  Capt.  Porter  and  the  author. 

A  four-s tate,  3  input  and  3  output  model  for  the  FPCC  aircraft 
was  developed  for  3  flight  conditions,  subsonic  (Mach  0.6,  0  ft), 
transonic  (Mach  0.9,  30,000  ft)  and  supersonic  (Mach  2.3,  40,000  ft). 
Chapter  IV  describes  the  model  evolution  from  the  original  6  state,  6 
input  and  6  output  model  to  its  current  version.  Appendix  A  contains  the 
non-linear  equations  that  are  used  by  the  FPCC  SIM  program  to  derive  the 
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state  equations. 


Chapter  V  contains  the  results  from  a  design  based  upon  the 
transonic  flight-condition.  Three  command  vectors  are  applied  to  the 
system:  positive  pitch  pointing,  vertical  translation  and  straight  climb. 
Very  good  results  are  obtained.  The  physical  limits  on  the  lift  flap  are 
exceeded  in  some  cases,  but  this  can  be  remedied  by  applying  a  ramped-up 
step  input  in  order  to  avoid  large  initial  transients.  The  design  is 
applied  to  both  the  subsonic  and  supersonic  models  to  determine 
robustness.  Very  good  results  are  obtained,  including  the  finding  that 
either  a  0.1  radian  vertical  translation  or  a  0.1  radian  positive  pitch 
pointing  maneuver  i6  too  severe  at  Mach  2.3  for  the  engine  to  sustain 
zero  velocity  change. 

Chapter  VI  is  essentially  a  sensitivity  study.  Each  of  the 
design  parameters  is  raised  or  lowered  individually  to  determine  the 
overall  effect  on  the  system.  This  also  has  the  secondary  purpose  of 
validating  the  original  design.  A  table  is  presented  at  the  end  of  the 
chapter  to  provide  the  designer  with  a  quick  overview  of  design  parameter 
influence . 

RBCOMMEHDATI OH S 

A  limitation  exists  in  the  method  for  choosing  a  measurement 
matrix  for  the  irregular  design.  Once  a  systematic  approach  i6 
developed,  this  design  method  should  find  universal  acceptance. 

Any  future  follow-on  effort  in  designing  a  flight  control  system 
for  the  FPCC  aircraft  should  include  the  lateral  dynamics  as  well  as  the 
longitudinal.  This  will  contain  4  additional  states,  3  inputs  and  3 
outputs.  The  feasibility  of  developing  separate  longitudinal  and  lateral 
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control  systems  can  also  be  studied. 

The  computer  program  MULTI  provides  the  designer  with  a  powerful 
tool.  The  usefulness  may  be  improved  by  inclusion  of  the  following 
points : 

(1)  The  use  of  random  access  files; 

(2)  The  ability  to  identify  transmission  seros; 

(3)  The  ability  to  calculate  the  closed-loop  asymptotic  roots; 

(4)  The  ability  to  determine  open-loop  atability; 

(5)  The  use  of  a  more  rigorous  error  protection  schcse; 

(6)  The  ability  to  derive  figures  of  merit. 
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APPENDIX  As  DEVELOPMENT  OP  EQUATIONS  OF  MOTION  FOR  THE  FFCC  AIRCRAFT 


The  longitudinal  equations  of  motion  for  the  aircraft  in  body  axes 
are  presented. 

TRANSLATIONAL  EQUATIONS 

u.  -  F  +  rV  -  qw  -  g  sin  ©  (A-l) 

"  -XB  1  1 

m 

w  «  F  +  pV  +  qu  +  g  cos  ©  cos  4  (A-2) 

1  -ZB  1  1 

m 

ROTATIONAL  components 

«  ■  "b  *  azz  -  rxx,rp  *  ;XZ  '  r2  -  V2)  <a-3) 


EULER  RATES 


©  =  q  cos  |  -  r  sin  j)  (A-4) 


EXTERNAL  FORCES  IN  BODY  AXES 


F  “X  cos  OK  -  Z  Bin oT  +  x  +  X  (A-5) 

*B  s  s  eng^  eng2 
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F  ■  Z  cosac  +  X  sin^C  +  z  +  Z 

zB  s  s 


engx  eng2 


(A-6) 


MOMENTS  IN  BODY  AXES 

and  Az  represent  c.g.  offset  from  m.r.p. 

M^-M+FAz  +  F^x  +  X  z  +X  z 

"B  s  xb“  zb"  eng1  engx  eng2  eng2 


-  Z  x  -  Z  x 

eng1  engx  eng2  eng2 


(A-7) 


AERODYNAMIC  FORCES  AND  MOMENTS  IN  STABILITY  AXES 


XS  =  '^VT  SC 


2  *  "  k=l  ^ 


GA-8) 


D  -  w  *  V 
\  C1 


(A-9) 


'  '/?  ^  S  CI 


(A-10) 


M  =  /0  Vm2  S  c  C 
S  /f  T  a 


(A-ll) 


DIRECT  PROPULSION  FORCES  IN  BODY  AXES 


eng 


=  m  _  F,  cos  9.  cosy 
1  v  f  1  j  engL 


-  C, 


r,  1  PV m 

d  jr  t 
AI 


(A-12) 


eng. 


X  =  <n  _  F,  cos  9 .  cos  'j  -  c  1  0  V  “  s 

-  t  F  2  j  eng  D  -y  T 


*2  I  2 ' 


(A-13) 


S  -  -01  _  F.  sin  9 , 

engx  F  1  3 


(A-14) 


eng. 


+L  F  F2  Sin  ®j 


(A-15) 
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AERODYNAMIC  FORCES 


(CT )  +  CT  jf  +  AC.  +  CT 

L  %d°  **  hfy  ' 


C-if  +  _£_  <  cr  <1  +  c  0C ) 

0  Jfi  ™  L~  L« 


Si  f2"J'2  2VT  '« 


(A-27) 


CD  “  CD, 


MIN 


* k  (cc”  V2  +  %  W  *  ^l4 


CT  )' 

lmf 


*  CD,  fjf  *  Ct  Sit 

hi,  ijf. 


(A-28) 


AERODYNAMIC  MOMENTS 


s,  •  =*_ « '♦  v  S„  *  %.  Sc  *  \  Jjf  ♦ 

fMF  ic  ijf.  1 


H,  (jf  *  J-  <  C„ 

*jf„  2  V_  q 


(A-29) 


MASS  PROPERTIES 

Airplane  mass  at  combat  weight 

M  *  902  slugs 

Moments  of  Inertia 

I  ■  57,500  slug-ft2 
xx 

I  -  130,000  slug- ft2 

yy 

I  ■  175,000  slug-ft2 
zz 

2 

I  ■  5,400  slug-ft 

xz 


(A- 30) 

(A-31) 

(A-32) 

(A-33) 

(A-34) 
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MOMENT  REFERENCE  POINT  (0.25  MAC) 


XMRP  “  3.5  ft 

C.G.  OFFSET  FROM  MRP 
Ax  -Ay  =Az  =  0 


(A-35) 


(A-36) 


THRUST  GEOMETRY 


x  =  x  =  -12.5  ft 

xeng^  eng2 


y  =  -5.16  ft 

eng1 


y  =  5.16  ft 

eng2 


2  =2  =0 
eng1  eng2 


THRUST  LINE  TOW- IN  ANGLE 


T  =  °c 


thrust  adjustment  factor 


-  0-g8 


(A-37) 

(A- 38) 

(A-39) 

(A-40) 


(A-41) 


(A-42) 
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5.  Option  Pre-requisites 


PROGRAMMER'S  MANUAL  FOR  PROGRAM  MULTI 


1 .  Introduction 

This  guide  provides  the  documentation  needed  for  future 
modifications  to  MULTI.  It  is  also  intended  to  help  a  programmer  analyze 
the  flow  structure  of  the  program  so  that  MULTI  can  be  modified  if 
unexpected  errors  occur.  This  manual  describes  the  overlay  structure, 
local  and  attached  subroutines  and  the  names  of  all  elements  of  the 
program. 

MULTI  is  written  in  FORTRAN  V  and  is  fully  documented  by  COMMENT 
statements  throughout  the  program.  The  program  structure  is  such  that 
programmers  familiar  with  FORTRAN  V  can  understand  MULTI 's  operation. 

The  programmer  should  have  a  full,  working  knowledge  of  the  theory  behind 
the  design  of  discrete-time,  error-actuated  controllers  for  linear, 
multivariable  plants  as  developed  by  Professor  Brian  Porter,  University 
of  Salford,  England  (Ref.  1,  2,  3).  The  programmer  should  also  have  used 
MULTI  interactively. 

MULTI  retains  the  code  from  the  University  of  Salford  which  deals 
with  measurement  matrices  for  the  C  and  G(0)  matrices.  In  order  to  avoid 
confusion,  these  matrices  are  not  utilized  in  the  current  design  code, 
and  thus  throughout  MULTI,  the  following  matrix  identities  exist: 

CM  (I,J)  -  C  (I.J) 

GMO  (I.J)  -  GO  (I.J) 

The  reader  should  obtain  a  MULTI  source  listing  before  continuing 
with  this  guide  and  have  a  CDC  FORTRAN  V  REFERENCE  MANUAL  (Ref.  14)  and  a 
CDC  LOADER  REFERENCE  MANUAL  (Ref.  17)  available.  Other  programming 
guides  which  may  be  helpful  to  a  potential  MULTI  programmer  include: 
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IMSL  Library  Reference  Manual  (Ref.  18) 

ASD  Computer  Center  GALCOMP  Plotting  Guide  (Ref.  19) 

Solution  of  Ordinary  Differential  Equation*  of  the  GDC 
6600/CYBER  74  Processor*  (Ref.  20) 

Computer  Printout  for  Subroutine  FLOTIT  (Ref.  21) 

2.  Description  of  Overall  Structure 

MULTI  is  written  to  provide  a  user  with  an  interactive  design 
tool  for  the  design  of  control  laws  needed  to  attain  tracking  and 
disturbance  rejection  in  a  multi-variable  plant.  It  is  necessary  for  the 
program  to  have  the  capability  to  evaluate  the  input  and  output  responses 
after  the  control  law  is  simulated.  In  addition,  it  is  necessary  for  an 
iterative  design  method  that  the  program  retain  all  parameter  values 
between  designs.  Finally,  the  program  must  be  fully  interactive. 

MULTl'must  fulfill  all  of  these  requirements  and  be  able  to 

t 

operate  in  the  limited  65,000g  words  of  memory  core  that  is  available  on 

AFIT's  INTERCOM.  In  its  original  form  MULTI  required  in  excess  of 

110,000Q  words  of  memory.  Thus  it  was  necessary  to  redesign  MULTI  using 
o 

en  overlay  structure  and  labeled  common  block. 

2 . 1  Overlay*  , 

Overlays  are  used  to  reduce  the  storage  requirements  of  large 
programs  by  dividing  the  program  into  modules.  All  modules  are  separate 
programs  in  their  own  right  and  are  linked  together  by  the  use  of  a  main 
executive  module. 

The  main  executive  module  and  all  common  variables  of  the  program 
constantly  remain  in  operational  core.  The  main  executive  directs  the 
program  flow  by  calling  the  primary  overlays  into  operational  core  as 
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they  are  needed. 

The  data  from  common  variables  is  passed  between  the  main  and 
primary  overlays  by  declaring  the  variable  in  COMMON  statements  in  the 
main  executive  overlay. 

The  CDC  FORTRAN  V  REFERENCE  MANUAL  and  the  CDC  LOADER  REFERENCE 
MANUAL  for  NOS  and  NOS/BE  can  provide  a  programmer  with  more  detailed 
information  on  overlays  and  COMMON  statements. 

2.2  MULTI'a  Overlay  Structure 

MULTI  is  composed  of  one  main  overlay  and  13  primary  overlays. 

The  main  overlay  provides  the  executive  directing  function  of  the 
program.  As  an  executive,  the  overlay  initializes  and  stores  data  in 
specified  common  blocks.  When  an  option  number  is  entered  by  the  user, 
the  main  overlay  checks  an  IF/ELSEIF  structure  to  attach  the  primary 
overlay  needed  to  satisfy  the  option  request.  The  main  overlay  also 
contains  subroutines  which  are  used  by  more  than  one  primary  overlay. 

A  main  overlay  description  is  given  by  Fig.  B-l . 


MAIN  OVERLAY 

Overlay  # 

Program  Marne 

Options  Included 

(0,0) 

MULTI 

AT-T. 

Overlay  # 

Subroutines 

Purpose 

(0,0) 

MATPR 

Prints  matrices 

QPRIHT 

Asks  if  data  is 
to  be  printed 

ANSWER 

Asks  if  data 

is  correct 

INVERT 

Inverts  matrices 

Fig.  B-l  Description  of  Main  Overlay 
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Initially  the-  *?ere  three  primary  overlays;  each  vas  chosen  to 
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conform  to  three  developmental  aspect*  of  controller  design:  data  input, 
control  law  systhesis  and  simulation.  However,  to  further  reduce  memory 
core  requirements,  the  overlays  and  related  subroutines  are  now  organised 
and  defined  as  shown  in  Fig.  B-2. 


PRIMARY  OVERLAYS 

Overlav  # 

Program  Name 

Options  Included 

(  1,0) 

OPTO 

#0  -  #9 

(  2,0) 

OFTIO 

#10  -  #13,  #15  -  #17,  #19 

(  3,0) 

0PT140 

#14  -  Unknown  Plants 

(  4,0) 

0PT14R 

#14  -  Regular  Plants 

(  5,0) 

OPT14I 

#14  -  Irregular  Planta 

(  6,0) 

OPT  18 

#18 

(  7,0) 

OPT20 

#20  -  #29 

(10,0) 

OFTPLT 

#30  -  #39 

(11,0) 

OPI31 

#31  -  #33 

(12.0) 

0PT34 

#34  -  #36 

(13,0) 

ERROR 

ALL 

(14,0) 

MEMORY 

#99 

(15,0) 

PRINT 

#100  -  #130 

Overlav  # 

Subroutines 

Purpose 

(  7,0) 

CLPASS 

Forms  different  equations 

TOUT 

Calculates  output  values 

Pig.  B-2 

Description  of  Primary  Overlays 

2.3  MULTI 

's  Data  Elements 

Labeled  common  blocks  are  used  in  MULTI  to  transfer  data  values 
between  overlays.  These  data  values  are  also  retained  in  memory  and 
determine  the  basic  need  for  memory  core  space.  Each  COMMON  block  is 
selected  so  that  only  the  data  required  by  the  primary  attached  to  the 
main  overlay  is  transferred  into  the  primary  overlay.  In  addition,  all 
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character  variables  oust  be  in  separate  blocks  and  not  in  numerical 
COMMON  blocks. 

The  numerical  sequence  of  the  blocks  corresponds  to  the  overlay 
structure  of  the  program.  In  COMMON  block  7,  the  postscripts  "A"  and  "S" 
refer  to  actuator  and  sensor.  Figure  B-3  shows  which  common  blocks  are 
used  in  each  overlay. 

Arrays  and  matrix  elements  which  are  defined  in  the  labeled 
common  blocks  are  dimensioned  in  the  same  statements  rather  than  with  an 
additional  set  of  DIMENSION  statements.  Arrays  and  matrices  which  are 
not  common  to  more  than  one  overlay  are  dimensioned  in  the  individual 
overlay. 

Initializations  for  the  common  blocks  is  accomplished  via  DATA 
statements  in  the  executive  overlay.  All  values  are  set  equal  to  zero 
except  as  follows: 

Actuators  and  sensors  are  set  to  be  the  transfer  function... 


9999 

a  ♦  9999 

Control  limits  are  set  equal  to..ft 
value  —  -1.0  x  1010 

Maximum  value  —  1.0  x  10*® 

Initial  values  for  all  other  arrays  and  matrices  are  entered  in 
the  individual  overlays  by  the  use  of  DO  loops  or  DATA  statements. 

2.4  Proxam  Labels 

The  program  labels  are  selected  to  reflect  certain  operations  in 
the  program  flow.  Each  label  quickly  identifies  the  code  associated  with 
the  label  and  provides  a  method  for  easy  editing  and  error  checking.  The 
MULTI  program  labels  are  defined  as  follows: 
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4 

0000  -  1999 

2000  -  2999 

3000  -  3999 
4000  -  4999 
8000  -  8999 
9000  -  9999 


1WQ SB 

Assigned  sequentially  in  the  progran  to  label 
general  operations. 

Assigned  to  correspond  to  individual  option 
transfers  (i.e.  OPTICS  #1  is  labeled  2001). 

Assigned  to  error  statenents  in  Overlay  (13,0). 

Assigned  to  FODfAT  statements. 

Assigned  to  error-free  exit  from  overlays. 

Assigned  to  CONTINUE  statenents  for  erroneous 
exit  from  overlays. 


2.5  Overlay  Structure 

Each  of  HULTI's  overlays  are  structured  in  the  same  basic  format. 


The  format  has  several  elements  beginning  with  overlay  identification  and 


ending  with  an  END  statement.  Each  overlay,  however,  may  not  contain  all 
elements.  The  format  used  in  MULTI' s  structure  is  as  follows: 


A.  Overlay  Beginning,  including... 

Overlay  Identification 
Program  lame 

CHARACTER,  INTEGER,  UAL  Declarations 
COMMON  Blocks 
DIMENSION  Statements 
DATA  Statements 

B.  Option  Flag  Checking 

C.  Option  Routing 

D.  Option/Overlay  Code 

Option  Flag  Initialisations 

E.  FORMAT  Statements 

F.  Overlay  Ending,  including... 

Error  Flag  Initialisations 
END  Statement 


Comment  cards  are  used  throughout  the  program  to  help  a 
programmer  interpret  the  overlays.  In  general,  four  lines  of  asterisks 
box  an  overlay-heading  comment  card,  two  rows  of  asterisks  lead  all 
subroutines  and  format  blocks  and  a  single  line  of  asterisks  separate  the 
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options  in  each  overlay.  Finally,  consent  carda  vitb  dashed  lines 
explain  various  operations  or  the  next  lines  of  code,  or  the  dashed  lines 
are  used  to  separate  sections  of  the  program. 

2.6  Overlay  Calls 

The  main  overlay  is  identified  vith  a  name  and  numerical 
designator: 

OVERLAY  (MULTI, 0,0) 

The  primary  overlays  are  identified  only  by  numerical 
descriptions : 

OVERLAY  (n,0) 

vhere  "n"  is  an  octal  number . 

Each  overlay  is  called  by  designating  the  main  overlay  name, 
MULTI,  and  the  numerical  identification  of  the  primary  overlay  vith  "nM 
in  decimal.  For  example,  to  attach  Overlay  (11,0)  the  program  call 
statement  is: 

CALL  OVERLAY  (MULTI, 9,0) 

vhere  9  is  the  decimal  equivalent  of  the  octal  number  11. 

This  aspect  of  overlay  usage  is  critical  and  is  not  defined  veil 
in  the  CDC  FORTRAN  V  REFERENCE  MANUAL. 

2.7  Option  Flaga  and  Krror  Flags 

MULTI  uses  a  simple,  but  effective,  method  to  assure  that  options 
are  accomplished  in  the  correct  order.  This  is  necessary  since  program 
flov  must  begin  vith  data  input  and  proceed  to  simulation.  As  each 
option  is  accomplished,  a  corresponding  option  flag  is  set  before  program 
flov  leaves  the  option. 

To  check  program  flov,  the  operational  code  beginning  each 
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The  lines  of  code  described  above  are  only  executed  once;  this 
occurs  vhen  the  program  is  first  begun.  After  MULTI  is  opened,  the 
program  flov  remains  in  the  option  request-loop  which  begins  at  statement 
#9000. 

3.2  Error  Checking 

All  program-determined  errors  set  an  error  flag  denoting  the  type 
of  error.  When  the  program  variable,  IERR,  is  not  equal  to  zero,  Overlay 
(13,0)  is  called  to  warn  the  user  of  the  error  condition. 

This  error  checking  is  the  first  operational  line  of  the  option 
request-loop  code. 

3.3  Option  Request-Loon  Code 

The  option  request-loop  begins  with  the  statement: 

FRUIT  '(//A)',  'OPT I<M,  PLEASE  #' 

and  is  used  to  obtain  the  option  number  request  from  the  user  and  direct 
the  program  flow  accordingly.  The  routing  is  accomplished  by  checking  a 
series  of  IF/ELSEIF  statements. 

There  are  two  additional  checks  accomplished  vhen  OPTION  #14  is 
selected  or  if  any  plotting  options  are  chosen. 

If  OPTION  #14  is  requested,  the  program  asks  the  user  to  select 
the  type  of  design  that  is  to  be  accomplished.  The  single-character 
variable,  METHOD,  is  used  to  further  route  the  program  flow  to  obtain  the 
overlay  needed  for  the  computations. 

When  OPTIONS  #30-39  are  chosen.  Overlay  (10,0)  is  called  to  form 
a  plotting  matrix  for  terminal  plots  or  CALCOMP  plots.  Upon  completion 
of  OVERLAY  (10,0),  the  program  returns  to  the  main  executive  for  more 
option  checks.  If  a  terminal  plot  is  requested,  Overlay  (11,0)  is  called 
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to  accomplish  the  terminal  plot.  If  OPTIONS  #34-36  are  requested  to 
obtain  a  CALCOKP  plot.  Overlay  (12,0)  is  brought  into  operational  core  to 
produce  the  PLOT  file. 

3.4  Main  Executive  Overlay  Subroutines 

There  are  four  subroutines  attached  to  the  main  overlay.  These 
subroutines  are  used  by  more  than  one  overlay  and  therefore  remain  in 
operational  core  at  all  times. 

SUBROUTINE  MATPR  (TE.IR.IC)  Subroutine  MATPR  is  used  for  printing  all 

matrices.  The  subroutine  has  three  parameters  as  described  below: 

TR  -  k  real  matrix  of  maximum  dimension  10x10. 

IE  -  An  integer  value  denoting  the  #  of  rows  of  the  matrix. 

IC  -  An  integer  value  denoting  the  #  of  columns  of  the  utrix. 

RtlRBonTIMB  OPRJWT  fmn,*)  This  subroutine  is  used  to  ask  if  data  should 

be  printed.  There  are  two  parameters  used  in  the  subroutine  call: 

CHAR  -  Character  string  with  maximum  of  30  characters. 

*  -  Line  number  denoting  where  program  flow  should  go  if  data 

is  not  to  be  printed. 

If  the  dat-a  is  to  be  printed,  the  code  directing  the  printing 
should  follow  immediately  after  the  subroutine  call. 

SUBROUTINE  ANSWER  (*.*)  ANSWER  is  used  after  the  program  echoes  input 

data  to  the  user  and  askB  if  the  data  is  correct.  If  the  answer  from  the 

user  is  affirmative,  the  program  continues  with  the  statements  following 

the  subroutine  call.  If  the  data  is  not  correct  or  the  option  is  to  be 

aborted,  the  subroutine  parameters  direct  the  program  flow  as  follows: 

1st  *  -  Line  number  denoting  incorrect  data  directing  return  to 
data  input  point. 

2nd  *  -  Lina  number  denoting  option  abort  directing  return  to  end 
of  option. 

For  the  second  asterisk,  the  line  number  should  be  selected  so  as 
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to  set  I ERR  equal  to  zero  and  so  that  the  option  flag  is  not  set 


MimonTlMg  INVERT  (A.AIMV.l.IA.*)  This  subroutine  is  used  to  invert  a 
matrix.  The  calling  parameters  are: 

A  -  A  real  matrix  to  be  inverted  of  maximum  dimension  10x10. 

AHV  -  The  resulting  real,  inversed  matrix. 

>  -  The  #  of  rows  in  the  A  matrix. 

IA  -  The  —  »in«  rov  dimension  of  A  matrix  as  described  by 
the  external  program. 

*  -  Line  number  denoting  where  program  flow  should  be  directed 
if  the  A  matrix  cannot  be  inverted. 

The  line  number  parameter  is  selected  to  route  program  control  to 
a  point  where  an  error  flag  is  set  corresponding  to  a  program-directed 
error  statement.  This  error  statement  indicates  which  matrix  is  not 
inver table. 

Subroutine  INVERT  accomplishes  the  matrix  inversion  by  use  of  the 
1MSL  surbroutine,  L1NV2F  (Ref.  18).  The  error  flag  from  Subroutine 
LINV2F  indicates  if  the  inversion  cannot  be  obtained.  Print  statements 
relating  the  problem  are  contained  in  Subroutine  INVERT  so  that  the  user 
does  not  have  to  ref or  to  the  IMSL  directives. 

4.  Description  of  MDLTI's  Primary  Overlays 

Program  MULTI  contains  13  primary  overlays.  The  program  code  is 
written  in  FORTAN  V  code  and  is  easy  to  understand.  Comments  are 
included  in  some  portions  of  the  code  to  help  explain  the  function  of 
certain  sections  of  the  program. 

Each  overlay  follows  the  structure  described  in  Section  2.5.  An 
additional  structural  aspect  at  the  beginning  of  Overlays  (1,0),  (2,0), 
(7,0),  (10,0)  and  (15,0)  is  the  routing  code  to  the  individual  options. 
Each  overlay  contains  lines  similar  to  the  following  line  taken  from 
Overlay  (7,0): 


CO  TO  (2021,2022,2023,2024,2025,2026,2027,2028,2029)  SOFT 
This  line  is  used  to  route  OPTIONS  #20  to  #29  to  the  lines  beginning  each 
option.  N0PT  is  an  integer  variable  which  is  equal  to  the  option  number 
the  user  requests.  Upon  entry  in  to  Overlay  (7,0),  NOPT's  value  is 
change  to  NOPT-20.  For  example,  if  N0PT  originally  was  Z3,  it  is  changed 
to  equal  3 ,  and  the  progrem  flow  is  directed  to  label  2023  which  is  the 
third  label  listed  in  the  GO  TO  statement.  Option  routing  in  the  other 
overlays  is  similar.  If  NOPT’s  modified  value  becomes  zero,  program  flow 
is  to  the  statements  directly  following  the  GO  TO  statement.  N0PT  must 
be  returned  to  its  original  value  before  leaving  an  overlay  in  case  an 
error  flag  has  been  set  during  the  overlay's  operation. 

The  reader  can  review  the  source  listing  for  a  description  of 
each  primary  overlay.  An  overview  of  each  of  the  overlays  follows  in 
this  section.  Section  5.  contains  a  listing  of  pre-requisites  for  each 
option. 


This  overlay  is  used  for  entering  data  values  which  describe  the 
plants.  The  codefor  the  option  in  this  overlay  consists  mainly  of  READ 
and  PRINT  statements. 

OPTION  #3  is  used  to  enter  the  plant  A,  B,  C  and  D  matrices.  It 
is  noteworthy  that  the  code  is  currently  written  so  that  if  a  matrix  data 
entry  is  made  incorrectly,  the  user  must  re-enter  the  entire  matrix 
rather  than  change  the  individual  element  value.  For  this  reason  each 
matrix  is  echoed  back  to  the  user  for  checking  by  the  calls  to 
Subroutines  MATPR  and  ANSWER. 


OPTIONS  #4  and  #5  include  code  to  reset  actuator  and  sensor 
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values  to  the  transfer  function: 


9999 

a  ♦  9999 

in  order  to  eliminate  their  effect  in  the  simulation. 

OPTIONS  #8  and  #9  read  plant  data  from  a  local  file.  The  data  is 
read  by  opening  a  file,  LFN,  specified  by  the  user.  LFN  is  a  variable 
that  can  be  up  to  30  characters  in  length.  The  file  is  first  opened  on  a 
program-selected  unit  number,  the  data  is  read,  and  the  unit  is  then 
closed.  It  is  essential  that  the  unit  be  closed  at  the  end  of  the  option 
so  that  multiple  entries  to  the  option  can  occur.  This  might  be  desired 
if  the  user  wishes  to  introduce  new  plant  data  values  for  a  different 
flight  condition. 

4.2  Overlay  (2.0)  -  Options  #10  -  #19 

This  overlay  contains  the  code  for  options  #10  through  #19, 
except  that  the  code  for  OPTIONS  #14  and  #18  are  contained  in  subsequent 
overlays.  OPTIONS  #15  thru  #17  are  not  used. 

The  code  for  the  options  in  this  overlay  is  mostly  composed  of 
PRINT  and  READ  statements  pertaining  to  design  parameters. 

OPTION  #19  contains  the  FORTRAN  lines  to  read  design  parameters 
from  a  local  file  which  must  be  specified  by  the  user.  The  operation  of 
data  file  opening  and  closing  is  given  in  the  discussion  of  Overlay 
(1,0). 

4.3  Overlay  (3.0)  -  Option  #14  -  Unknown  Plants 

The  controller  matrices,  K0  and  K1 ,  for  unknown  plants  are 
calculated  in  Overlay  (3,0).  In  addition  to  the  COMMON  blocks  of  the 
overlay,  there  are  four  additional  scratch  matrices  (CM,  W,  WW  and  ZZ) 
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dimensioned  end  initialized.  The  matrix  values  for  these  matrices  are 
not  retained  after  the  overlay  is  terminated. 


The  unknown  design  can  be  accomplished  by  use  of  the  plant  G(0) 
matrix  or  by  forming  this  matrix  from  the  plant's  A,  B,  C  and  D  matrices. 
In  the  latter  case  the  calculations  are: 


Prrtflr—  Calculation 

W-A-1 

ZZ-W*B 


GO-D  -  C*ZZ 
VW-D  -  C*ZZ 


Actual  Calculation 

W-A_1 

ZZ-A_1*B 

G(0)-D-*A  “  1  *  B 
GM(0)-D  -  C*A_1  *  B 


At  this  point  either  the  calculated  G(0)  matrix  or  the 
user-originated  matrix  is  used  to  create  the  controller  matrices.  The 
calsulations  are: 


Program  Calculation 

wmxT1 

S«W*GAMA 


Kl-S*KPSLOB 

KO-ALPHA*S*KPSLOM 


Actual  Calculation 
WfG(O)  _1 


K1-G(0)_1*GAMA*KPSIL0B 

rO-ALPHA*G(0)_1*CAMA*KPSILOH 


The  single  character  variable,  METHOD,  is  set  at  the  end  of  the 
overlay  such  that : 

METHOD  -  'O' 

denoting  unknown  plant  design  for  normal  overlay  termination  or: 

METHOD  -  'X* 


denoting  abnormal  overlay  termination. 

4.4  Overlay  (4.0)  -  Option  #14  -  tegular  Plants 

This  overlay  calculates  controller  matrices,  K0  and  K1 ,  for 
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regular  plants.  There  are  three  scratch  matrices  (VV,  W  and  ZZ)  which 
are  dimensioned  and  initialized  at  the  beginning  of  the  overlay.  The 
matrices  are  destroyed  when  the  overlay  is  completed. 

After  pre-requisite  checks,  the  controller  matrices  are  formed  by 
the  following  calculations: 

Program  Calculation  Actual  Calculation 

W«C*B  V?-C*B 

zz-(w)-1  zz-(c*b)”1-b~1*c~1 

VW-ZZ*GAMA  WW-B_1*C-1*GAKA 

KO-VW*ALFHA*KPSLOB  KO-B~1*C~1GAMA*ALPBA*EPSILOH 

Xl-WV*EPSLOB  K1 "B-1 *C-1 *GAMA*EPSILOB 

For  a  normal  termination  _o  the  overlay,  the  single  character  variable, 
METHOD,  is  set  equal  to  "R"  denoting  a  regular  plant  design.  If  abnormal 
termination  occurs,  METHOD  is  set  equal  to  "X". 

4.5  Overlay  (5.0)  -  Option  #14  -  Irregular  Plants 

An  irregular  plant  design  of  the  controller  matrices,  K0  and  K1 , 
is  accomplished  in  Overlay  (5,0).  Thre  are  four  temporary  matrices  (VV, 
WW,  ZZ  and  VWZ)  which  are  dimensioned  and  initialized  for  use  only  while 
this  overlay  is  attached  to  operational  core. 

For  this  type  of  design  the  program  requires  a  measurement 
matrix,  MM,  to  augment  the  rank  deficiencies  of  the  plant  matrices. 
After  pre-requisite  checks,  the  next  block  of  program  code  determines  if 
the  measurement  matrix  already  exists,  if  the  user  must  initilize  MM,  or 
if  the  user  desires  to  re-initialize  the  values  of  the  matrix.  The 
measurement  matrix  values  are  originally  a1. 1  set  equal  to  zero  in  the 
main  overlay.  If  the  user  has  not  initialized  the  values  by  selecting 


150 


OPTIONS  #9  or  #18,  the  matrix  values  must  be  entered  throught  this 
overlay. 

The  following  equations  then  occur: 


Prosram  Calculation 

Actual  Calculation 

W-C  ♦  MK*A 

F-C  +  M*A 

w-vw2 

ZZ-W*2 

f2-c2  +  m*a12 

-1 

ww-zz 

WW-F, 

-1 

ZZ-W*GAMA 

ZZ-F2  *gaka 

vw-b2 

WW“B2 

-1 

_  -1 

vwz-ww 

WZ-B  1 

-1  -1 

WW-VWZ*ZZ 

VW-B2  *F2  x*gana 

controller  matrices  are 

then  calculated  as  follows: 

Program  Calculation 

Actual  Calculation 

KO-WW*ALPHA*KPSLO> 

ko-b2_1*f2_1*gaka*alpha*kpsilo« 

Xl-WW*EPSLOH 

Kl-B2_1*F2_1*GAKA*KPSILOB 

The  overlay  terminates  by  setting  by  setting  the  single  character 
variable,  METHOD,  equal  to  "I",  denoting  an  irregular  plant  design,  if 
the  overlay  terminates  normally.  It  is  set  equal  to  "X"  if  the  overlay 
terminates  abnormally. 

VMULFF,  a  high  precision,  IMSL  Library  subroutine,  is  used  in 
this  overlay  to  multiply  two  matrices  together.  When  using  this 
subroutine,  the  second  matrix  which  is  entered  into  the  subroutine  is 
overwritten  during  the  matrix  multiplication  process.  Thus,  the  contents 
of  this  matrix  are  destroyed. 
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4.6 


OverltT  (6.0)  -  Option  #18 

Overlay  (6,0)  is  provided  for  two  purposes.  One  purpose  is  to 
create  a  measurement  matrix  for  use  in  an  irregular  plant  controller 
design,  and  the  other  purpose  is  to  provide  a  tool  for  row  and  column 
operations  on  any  matrix. 

This  overlay  uses  four  scratch  matrices  (VV,  WW,  ZZ  and  VWZ) 
during  its  computations.  These  matrices  are  discarded  when  the  option 
ends. 

The  integer  variable,  ISKIP,  is  set  equal  to  1 ,  2  or  3  depending 
on  the  user's  choice  of  measurement  matrix  creation,  row  and  column 
operations  on  [  C  *  B  ] ,  or  row  and  column  operations  on  any  matrix  of 
interest.  The  routing  for  the  different  choices  is  determined  by  an 
IF/ELSEIF  block  of  code. 

Before  creating  a  measurement  matrix,  the  user  must  obtain  the 
feedback  gain  matrix,  KI ,  from  CESA  OPTION  #38  (Ref.  9).  After  this 
matrix  is  read  into  the  program  as  matrix,  VWZ,  and  is  verified  by  the 
user,  the  measurement  matrix  is  calculated. 


Proem  Calculations 

Actual  Calculations 

MM-Cj 

MM-Cj 

ww-c2 

ww-c2 

ZZ-MM  -  WW*VWZ 

ZZ-Cj  -  Cj*Kl 

■**12 

"**11 

W-A12 

WV-Aiz 

W^W*V¥Z  -  MM 

"■*12“  -  A11 
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wv-w"1 

HKZ«W 

The  program  code  for  performing  row  and  column  operations  on  [  C 
*  B  ]  and  performing  row  and  column  operations  on  any  other  matrix  is 
Identical.  The  only  difference  between  the  two  choices  is  the  matrix 
that  is  used  in  the  operations.  In  the  first  case,  the  [C  *  B  ]  matrix 
is  automatically  formed  by  the  program  without  any  input  from  the  user. 
In  the  second  case,  the  user  muBt  supply  a  matrix  of  maximum  dimension 
10x10. 

The  program  provides  for  four  choices  of  operations  on  a  matrix. 
The  user's  choice  of  operations,  which  is  identified  by  the  integer 
variable,  ISKIP,  directs  program  flow  through  an  IF/ELSEIF  block.  When  a 
new  matrix  is  formed,  it  is  echoed  back  to  the  user.  The  program  then 
returns  to  the  operation-selection  point. 

4.7  Overlay  (7.0)  -  Options  #20  -  #29 

All  aspects  of  the  system  simulation  are  contained  in  this 
overlay.  The  overlay  also  contains  two  subroutines,  CLPASS  and  YOUT, 
which  are  used  in  OPTION  #26.  Since  subroutine  CLPASS  is  used  as  a 
parameter  in  a  CALL  statement,  it  must  be  declared  as  an  EXTERNAL  in  the 
beginning  of  the  overlay.  There  are  two  scratch  matrices  (CM  and  F)  and 
eight  scratch  arrays  (IWORK,  WORK,  X,  Z,  Y,  E,  AX  and  MMAX)  that  are 
dimensioned,  initialized,  used  and  discarded  as  the  overlay  operates. 
The  F  matrix  of  this  option  should  not  be  confused  with  the  F  matrix 
generated  by  OPTION  #14  for  irregular  plants.  There  is  also  a  character 
variable  called  STRING  hich  is  declared  at  the  beginning  of  the  overlay. 


-  AU> 

■Me,  -  -  Ajj) 
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Before  the  variable  NOPT  routes  the  program  flow  to  the  different 
options  of  the  overlay,  the  CM  and  F  matrix  values  must  be  set.  Although 
these  matrices  are  only  superficial,  it  is  necessary  that  they  be  equal 
to  these  values: 

CM  (I,J)  -  C  <1,J> 

F  (I,J)  “  Identity  Matrix 

These  two  matrices  are  found  in  the  original  program  code  obtained  from 
the  University  of  Salford  and  are  used  in  OPTION  #26. 

OPTION  #20  to  #25  and  OPTION  #27  are  mainly  FORTRAN  READ  and 
PRINT  statements  which  are  used  in  obtaining  state  and  integrator  initial 
values,  the  output  command  vector,  simulation  time  parameters  and  the 
control  input  limitations.  0PIT0N  #27  also  contains  the  code  to  remove 
control  limits  by  setting  the  control  input  minimum  value  to  -1.0  x  10*® 
and  the  control  input  maximum  to  1 .0  x  10*®. 

OPTION  #29  opens  unit  number  40  to  read  the  design  parameters 
from  a  user-specified  local  file.  The  operation  of  this  option  is  the 
same  as  the  operation  of  OPTION  #9  which  is  discussed  in  Section  4.1. 

The  heart  of  the  simulation  is  contained  in  the  OPTION  #26  code. 
As  previously  stated,  this  code  is  nearly  identical  to  the  code  received 
from  the  University  of  Salford.  The  option  begins  by  checking  for  option 
flags  which  may  not  be  set  and  by  checking  for  the  presence  of  actuators, 
sensors  and  contol  limits.  Next,  the  total  number  of  states  is 
determined  by  adding  the  number  of  states  from  the  sensors,  actuators  and 
plant  matrices. 

There  are  three  major  loops  in  this  option.  The  outer  loop  for 
each  sampling  time  is  blocked  by  label  #1285;  the  middle  loop  based  on 


sampling  period  is  blocked  by  label  #1250;  the  inner  loop  based  on  step 
time  is  blocked  by  label  #1240. 

QgTfR  LOOP  The  outer  loop  occurs  since  the  program  allows  the  simulation 
to  be  run  with  more  than  one  sampling  time.  Each  time  a  diferent 
sampling  time  is  run,  the  program  must  initialize  several  variable  and 
array  values.  A  calculation  of  the  number  of  incremental  points  is  also 
made  based  upon  the  total  simulation  time,  TT,  the  sampling  time,  SAMT, 
and  the  step  time,  ST.  The  following  formulas  determine  the  number  of 
total  time  increments,  NT: 

■II  ■  Integer  Value  of  (  TT  /  SAMI  ♦  0.5  )  (H) 

■TT  ”  Integer  Value  of  (  SAMI  /  ST  ■  0.5  )  (B-2) 

■T  -  *1*  *  BTT  (B-3) 

If  the  step  time  is  chosen  to  be  greater  than  the  sampling  time,  NTT  is 

set  to  a  value  of  one  prior  to  the  calculation  of  total  number  ot  time 

increments  in  Eq.  (B-3).  The  diagram  in  Fig.  B-4  shows  the  relationship 

between  TT,  SAMP  and  ST. 

If  the  total  number  of  time  increments  is  greater  than  100,  the 
simulation  exceeds  the  dimension  of  the  matrices  which  hold  the  input  and 
output  values.  Thus,  an  integer  variable,  IPACK,  is  calculated  to 
determine  the  rate  at  which  data  is  packed  into  a  matrix  with  100  rows. 
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Fig.  B-4  Relationship  Between  Total  Tune,  Sampling  Tune  and  8tep  Time 
If  there  are  125,  or  less,  time  increments,  the  data  from  the  first  100 
increments  is  retained.  If  there  are  more  than  125  increments,  IPACK  is 
calculated  as: 

IPACK  ■  [  Integer  Value  of  (  HT  /  100  )  ]  ♦  1  (B-4) 

As  an  example  of  the  packing  process,  if  there  are  300  time  increments, 
IPACK  is  calculated  as  three,  and  every  third  input  and  output  value  is 
retained. 

The  last  part  of  the  code  for  the  outer  loop  sets  the  integer 
variable,  ISKIP,  which  is  used  to  suppress  the  time-sequential  printout 
of  input  and  output  values  as  the  solutions  to  the  state  differential 
equations  are  found. 

irnmi-K  LOOP  The  middle  loop  is  passed  through  KIM  times  as  calculated  in 
Eq.  (B-l).  First,  the  ouput  measurement  vector,  YM,  is  formed.  In  MULTI 
this  vector  is  equivalent  to  the  output  vector,  Y,  and  is  calculated 
from 


I  -  fl*  *  I 


The  error  vectors  are  then  formed  as  follows: 

For  Irregular  Plants... 

Program  Calculation  Actual  Calculation 

AX-Aj*X 
MMA1“MM1  *Aj *1 
S-Q 

E-V  -  Q  -  Mf1*A1*X 

Actual  Calculations 
S-F*Q  8-0 

E-V  -  S  E-V  -  Q 

For  unknown  plants  the  controller  is  then  formed  from: 
D-(rg*K  +  lQ*Z)  /  SAMT  (B-6) 

and  for  all  othe  plants  the  controller  iB  formed  from: 

U-(K0*E+K1*Z)*  SKAT  (B-7 ) 

However,  before  these  inputs  are  used  in  the  simulation  they  are  compared 
to  the  current  input  control  limitations  and  modified  if  necessary. 

A  call  to  Subroutine  YOUT  is  the  final  line  of  code  in  the  middle 
loop.  This  subroutine  call  determines  the  controlled  output  vector  from 
the  states  of  the  system  at  any  sspecified  time  and  is  described  at  the 
end  of  this  section. 

IHML  LOOP  The  inner  loop  is  passed  through  NTT  times  as  calculated  by 
Eq.  (B-2).  The  first  section  of  the  loop  is  concerned  with  the  data 
packing  procedure  described  earlier  and  sequentially  sets  the  row  indices 
for  the  input  and  output  matrices,  UP  and  YP.  These  two  matrices  are 
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| _  NT-Tolal  Number  of  Time  Increments 

Fig.  B-5  Representation  of  the  3-D  Input  Matrix,  DP 

retained  throughout  MULTI  and  contain  the  input  and  output  data  as 
calculated  by  OPTION  #26.  The  matrices  are  3-dimensional  where  the  row 
dimension  is  equivalent  to  an  input  or  output  vector  at  each  time 
increment  The  first  column  holds  the  time  increment  values,  T,  while  the 
other  columns  bold  the  sequential  data  for  each  input  or  output.  The 
third  dimension  is  equivalent  to  the  sampling  time.  See  the  diagram 
presented  in  Fig.  B-5. 

MULTI  uses  Subroutine  ODE  (Ref.  20)  from  the  ASD  CC6600  Library 

to  solve  the  set  of  differential  equations  formed  by  Subroutine  CLPASS. 

The  CLPASS  subroutine  is  discussed  at  the  end  of  this  section.  Currently 

-4 

the  precision  limits  for  ODE's  computations  are  set  at  1.0  x  10  .  If  an 

error  occurs  during  the  solution  process  an  error  message  is  printed  to 
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alert  the  user  of  the  condition.  The  program  provides  the  error  flag 
number  to  the  user  and  states  that  the  ODE  Manual  is  to  be  referenced. 


The  inner  loop  is  ended  by  a  call  to  Subroutine  YODT  to  form  a 
new  output  vector,  Y,  from  the  values  returned  from  the  ODE  subroutine. 
When  the  inner  loop  is  complete,  a  new  integrator  vector  is  calculated 
and  program  flow  returns  to  the  middle  loop. 

When  all  time  increments  for  the  total  simulation  time  have  been 
processed  through  the  inner  and  middle  loops,  the  program  code  for 
printing  control  limit  information  is  entered.  A  "simulation  complete" 
message  is  then  provided  and  the  outer  loop  is  run  again  until  each 
sampling  time  is  simulated. 

Before  the  option  is  completed,  the  integer  variable  NT  is  set 
equal  to  the  actual  number  of  data  points  in  the  input  and  output 
matrices  of  the  last  simulation.  This  is  necessary  so  that  the  plotting 
routines  receive  the  correct  number  of  data  values  in  the  matrices. 
However,  since  NT  does  not  have  a  value  for  each  sampling  time  and  an 
input/output  matrix  can  be  formed  where  the  first  sampling  time  dimension 
has  less  data  values  than  the  second  sampling  time  dimension,  a 
termination  error  can  result  unless  MULTI  operates  under  one  of  the 
following  procedures: 

(a)  MULTI  should  be  run  with  all  sampling  times  entered  in 

order  of  increasing  value,  and  step  time  should  be  set 
so  as  to  be  equal,  or  less  than,  the  smallest  sampling 
time. 

or  (b)  MULTI  should  be  run  using  only  one  sampling  time. 
gUBlOUTIME  CLP  ASS  (T.X.ZDOT) 

This  subroutine  is  nearly  identical  to  the  CLFASS  subroutine 
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attached  to  the  PAK200  program  (Ref.  4).  Minor  modifications  are  made  in 
some  parts  of  the  code  so  that  it  may  be  used  as  an  external  program  in 
Subroutine  ODE. 

This  subroutine  forms  the  differential  equations  from  the  plant, 

actuator,  sensor  and  error  integral  matrices.  The  parameters  are: 

T  -  Time 
X  -  State  Vector 

XDOT  -  Derivative  of  State  Vector,  X 

CLPASS  defines  the  values  of  the  array,  XDOT,  at  time  T.  Array  X 
includes  all  the  states  of  the  composite  system.  This  array  contains  the 
plant  states  first,  followed  by  the  actuator  states  and  then  by  the 
sensor  states. 

The  system  input,  U,  is  applied  to  the  actuators.  The  output  of 
the  actuators,  W,  provides  the  input  to  the  plant.  Finally,  the  plant 
output,  YM,  is  the  input  to  the  sensors,  and  the  sensor  output,  Q,  is  the 
actual  system  output.  Fig.  B-6  helps  the  reader  visualize  this 
arrangement . 

SUBROUTINE  TOOT  (X.Y.W.C.D) 

Subroutine  YOUT  is  a  modified  version  of  the  YOUT  subroutine 
provided  with  PAK200.  MULTI's  YOUT  subroutine  uses  a  COMMON  block  to 
pass  data  into  the  subroutine  rather  than  passing  the  data  through  the 
calling  parameters. 

This  subroutine  determines  the  controller  output  vector  from  the 

states  of  the  system  at  any  time.  The  parameters  in  the  subroutine  call 

are  defined  as: 

X  -  System  State  Vector 
T  -  Controlled  Output  Vector 
W  -  Actuator  Output  Vector 
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C  -  System  Output  Matrix 
D  -  System  Direct  Output  Matrix 

4.8  Overlay  (10.0)  -  Options  #30  -  #39 

When  any  type  of  plot  is  to  be  generated,  MULTI  uses  Overlay 
(10,0)  to  generate  a  plotting  matrix,  PLMAT,  which  is  filled  with  data 
values  from  the  input  and/or  output  value  matrices,  UP  and  YP.  After  the 
PLMAT  is  formed  the  actual  terminal  plot  is  generated  by  Overlay  (12,0). 
These  two  overlays  are  discussed  in  the  next  two  sections. 

It  is  pointed  out  to  the  reader  that  COMMON  blocks  B13  and  B13A 
change  the  character  representations  of  the  integer  variable,  NT,  and  the 
matrices,  YP  and  UP.  In  this  overlay  and  Overlays  (11,0)  and  (12,0), 
these  variable  names  are  referred  to  as  N,  Y  and  P,  respectively. 

The  FORTRAN  coding  in  this  overlay  is  very  condensed  and  the 
program  flow  is  very  complex.  The  coding  is  condensed  so  that  the 
overlay  can  operate  successfully  while  using  minimum  core  space.  The 
core  space  requirement  is  critical  since  this  overlay  uses  the  three 
largest  matrices  of  the  program,  YP,  UP  and  PLMAT.  The  program  flow  is 
complex  since  MULTI  provides  (1)  terminal  plots  and  CALCOMP  plots,  (2) 
four  types  of  each  kind  of  plot,  (3)  a  long  version  and  two  shortened 
versions  of  the  plot  request.  All  these  provisions  are  completed  in  this 
single  overlay.  In  addition,  if  PLMAT  for  different  flight  conditions  is 
generated,  the  program  flow  loops  between  this  overlay  and  the  other 
options  of  the  program  until  PLMAT  is  completed.  In  this  case,  the 
overlay  is  required  to  keep  a  count  of  the  number  of  times  that  the  plot 
sequence  has  been  entered. 

To  correctly  complete  a  plotting  matrix,  PLMAT,  this  overlay 
extracts  columns  of  data  from  the  U  and  Y  matrices.  The  columns  are 


transferred  to  PLMAT  in  a  pre-designated  order  after  the  user  enter*  the 
value*  for  the  variable*  ICODE,  IFLTCN ,  LINES,  NDEPTH,  IDEPTH,  CHOICE, 
NDM,  ICOLMN  and  ICLM.  The  values  for  these  variables  may  be  set  by  the 
program,  rather  than  by  the  user,  depending  on  the  type  of  plot 
requested.  The  definitions  for  these  variables  are  found  in  the  last 
section  of  this  guide. 

After  the  PLMAT  matrix  is  complete,  the  integer  variables,  NDM 
and  ICLM,  are  reset  to  the  values  required  for  correct  operation  in 
Overlays  (11,0)  and  (12,0).  This  ends  the  overlay. 

4.9  Overlay  (11.0)  -  Terminal  Plot 

MULTI  has  the  capability  of  producing  terminal  plots  by  calling 
upon  Subroutine  PLOTIT.  PLOTIT  is  a  subroutine  designed  by  Major  Michael 
R.  Stamm  of  the  Department  of  Physics  at  the  Air  Force  Institute  of 
Technology  (Ref.  21). 

This  overlay  begins  with  the  addition  of  a  new  COMMON  block 
called  SCALIT  which  is  common  only  to  the  PLOTIT  subroutine.  This  COMMON 
block  provides  for  the  initialization  of  the  terminal  plot's  calling 
parameters.  As  discussed  in  Section  4.8,  this  overlay  also  renames  the 
variables,  NT,  UP  and  YF  as  N,  U  and  Y,  respectively. 

Overlay  (11,0)  must  also  provide  the  PLOTIT  subroutine  with  the 
independant  and  dependant  data  arrays,  XAXIS  and  YAXIS,  plus  the  minimum 
and  maximum  values  of  there  arrays.  The  XAXIS  array  is  always  filled 
with  values  from  the  first  column  of  the  output  matrix,  Y.  The  YAXIS  is 
sequentially  set  equivalent  to  the  columns  of  the  PLMAT  matrix  as  the 
PLOTIT  subroutine  is  called.  PLOTIT  is  called  until  all  the  data  in  the 
PLMAT  matrix  is  scanned  and  entered  into  the  terminal  plot  routine.  The 
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terminal  plot  is  then  generated. 

The  overlay  is  ended  vith  a  series  of  PRINT  statements  which 
correspond  to  symbols  on  the  terminal  plot  curves  to  the  input/output 
numbers  of  the  plant  states. 

4.10  Overlay  (12.0)  -  CALCOMP  Plot 

MULTI  has  the  capability  of  producing  CALCOMP  plots  that  can  be 
routed  to  the  CALCOMP  plotter  upon  program  termination.  The  various 
CALCOMP  subrout inee s  used  are  well  documented  in  the  CALCOMP  user  guide 
(Ref.  19).  Only  five  plots  may  be  put  on  a  plot  file  at  one  time.  Once 
this  maximum  is  reached,  the  user  must  exit  MULTI  and  route  the  plot  file 
to  the  plotter.  MULTI  keeps  count  on  the  number  of  plots  on  file  and 
prompts  the  user  when  the  maximum  of  five  is  reached. 

The  overlay  begins  by  incrementing  the  plot  number  counter.  This 
is  followed  by  setting  the  XAXIS  equal  to  the  first  column  in  the  YP 
matrix  which  contains  the  time  increment  listing.  The  PLMAT  matrix  iB 
then  scanned  for  its  minimum  and  maximum  values  in  order  to  accurately 
set  the  YAXIS  range.  A  title  with  a  20  character  maximum  is  then 
entered. 

The  CALCOMP  subroutines  are  then  begun.  The  XAXIS  and  YAXIS 
arrays  are  scaled  for  their  minimum  and  maximum  values.  The  proper 
columns  from  the  PLMAT  array  are  chosen  and  plotted.  Finally,  a  box  is 
drawn  around  the  entire  plot. 

The  last  section  of  Overlay  (12,0)  is  the  warning  message 
concerning  the  current  number  of  plots  on  file. 

4.11  Overlay  (13.0)  -  Error  Statements 

This  overlay  is  used  to  print  error  messages  for  all  operational 
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failures  detected  by  the  program.  The  operational  code  of  the  overlay 
consists  entirely  of  PRINT  statements. 


I 


k 


The  first  30  error  statements  refor  to  option  errors.  The  number 
of  the  error  statement  relates  directly  to  the  number  of  an  error  flag 
which,  in  turn,  corresponds  directly  to  an  option  number.  Error 
statements  #31  thru  #39  provide  information  about  matrix  errors.  Error 
statements  #40  thru  #44  give  miscellaneous  error  information.  When  an 
option  does  not  exist,  a  PRINT  statement  is  still  included  to  provide  an 
operational  check  of  MULTI .  These  error  statements  should  never  be 
printed. 

As  a  working  example,  consider  a  prerequisite  check  at  the 
beginning  of  an  option.  OPTION  #3  requires  the  option  flag  from  OPTION 
#2  to  be  set  before  the  plant  matrices  can  be  entered.  If  IFLAG(2)  is 
not  set  to  one  prior  to  selecting  OPTION  #3,  I ERR  is  set  equal  to  two  and 
Overlay  (13,0)  is  called.  A  value  of  two  for  I  ERR  routes  the  program 
flow  to  Statement  #3002  which  prints : 

#  OF  STATUS,  INPUTS  &  OUTPUTS  KISSING. . .SEE  OPTION  #2 
4.12  Overlay  (14.0)  -  Option  #99  -  Memory  Files 

When  the  user  selects  OPTION  #99  to  end  program  MULTI,  Overlay 
(14,0)  is  called  to  create  the  data  memory  files.  There  are  three  memory 
files  created. 

MEMO  is  created  to  hold  plant  data  as  directed  by  the  integer 
variable,  IPLANT.  If  IPLANT  equals  "1",  all  state,  sensor  and  actuator 
matrices  are  saved;  if  IPLANT  equals  "2",  the  plant  G(0)  matrix  is  copied 
to  the  file.  If  IPLANT  equals  "0",  plant  data  from  OPTIONS  #2  to  #5  are 
still  transferred  to  the  memory  file,  and,  when  the  file  is  read  back 
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into  MULTI,  the  IPLANT  value  alerts  the  user  that  the  file's  data  may  be 
incomplete. 


MEM10  ia  used  to  save  all  design  data,  and  MEM20  is  opened  to 
contain  the  simulation  data. 

The  overlay  ends  with  several  PRINT  statements  relating  the  file 
names  that  have  been  created  and  the  data  that  each  contains. 

4.13  Overlay  (15.0)  -  Options  #100  -  #139 

MULTI's  last  overlay  provides  data  value  printing  for  most  of  the 
options.  If  an  option  is  not  concerned  with  data  input  or  data  creation, 
the  program  uses  error  flag  #43  to  print: 

THERE  ARE  HO  VALUES  11  MEMORY  CORE  FOR  0PTI0B  # 

If  data  is  available,  the  values  are  printed  by  routing  the 
program  flow  to  the  proper  section  of  the  overlay.  All  printing  options 
are  numbered  by  adding  a  value  of  100  to  the  related  option  number  of  the 
main  program.  For  example,  data  values  entered  or  created  by  OPTION  #12 
are  printed  by  OPTION  #112. 

The  data  values  from  some  options  are  combined  to  form  a  single 
block  of  data  information.  The  data  from  OPTIONS  #11  and  #13  is 
generated  when  either  OPTION  #111  or  #113  is  selected.  The  same  occurs 
with  data  from  OPTIONS  #21  and  #22;  the  data  is  printed  by  using  either 
OPTION  #121  or  #122.  And  finally,  all  data  values  from  OPTIONS  #23  to 
#25  are  provided  by  selecting  any  option  number  from  123  to  125. 

The  code  for  the  overlay  is  simple  and  self-explanatory.  The 
reader  is  directed  to  the  MULTI  source  listing  for  further  detail  on  this 
overlay's  operation. 


166 


5 


MULTI* 


The  following  list  provides  the  pre-requisite  requirements  for 
options. 

OPTIOa  »  PRB-RIODISITK  OPTIC  * 


3 

2 

4 

2 

5 

2 

12 

2 

14 

2,  3,  11  TO  13 

18  (entries  1  &  2) 

5 

21 

2 

22 

2 

26 

14,  21  to  25 

29 

2 

31 

26 

32 

26 

33 

31  or  32 

34 

26 

35 

26 

36 

34  or  25 
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1.  Introduction  to  MULTI 

1.1  MULTI'e  Input  Modes 
OPTION  Mode 

DATA  Mode 
QUERY  Mode 

1.2  MULTI's  Options 

2.  Complete  Description  of  MULTI's  Options 

2.1  BLOCK  1  -  Plant  Input  Options 

2.2  BLOCK  2  -  Design  Input  Options 

2.3  BLOCK  3  -  Simulation  Options 

2.4  BLOCK  4  -  Plotting  Options 

2.5  BLOCK  5  -  Printing  Options 

3.  Summary  of  MULTI's  Options 

4.  Example  for  Unknown  Plants 
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USER'S  MANUAL  FOR  PROGRAM  MULTI 


Overview  of  MULTI 

MULTI  is  an  interactive  program  which  enables  a  computer  user  to  design 
and  simulate  digital,  multivariable  control  laws  for  discrete,  linear, 
multiple-output  systems.  The  three  design  methods  incorporated  in  this 
program  were  developed  by  Professor  Brian  Porter  of  the  University  of 
Salford,  England.  To  help  the  user  become  quickly  familiar  with  the 
program,  the  following  overview  is  provided. 

MULTI  is  able  to  design  and  simulate  control  laws  for  three  types 
of  plants: 

(1)  Regular  -  The  linear,  multivariable  plant  dynamics  are 
described  by  the  usual  state  and  output  equations,  and  the  first  Markov 
parameter,  CB,  has  full  rank  (Ref.  2). 

(2)  Irregular  -  The  plant  dynamics  are  described  by  state  and 
output  equations,  and  the  first  Markov  parameter  does  not  have  full  rank. 
Thus,  the  system  must  be  augmented  by  a  measurement  matrix  so  that  a 
control  law  may  be  developed  (Ref.  3). 

(3)  Unknown  -  The  plant,  Btate  and  output  equations  and  unknown, 
but  the  steady  state  transfer  function  matrix,  G(0),  is  obtainable  from 
off-line  tests  if  the  open-loop  plant  is  stable.  This  method  is  also 
applicable  to  known  plants  (Ref.  1). 

All  control  law  designs  may  be  evaluated  by  a  discrete-time 
simulation  of  the  system,  actuator,  and  sensor  equations.  The  user  may 
elect  to  obtain  a  terminal  plot  or  a  Calcomp  plot  of  the  system  input 
and/or  output  responses. 

MULTI  focuses  on  Porter ' s  digital  design  as  opposed  to  a 
continuous  design  method.  At  various  points  during  the  design  and 
simulation  process,  the  program  can  provide  intermediate  information 
about  system  matrices  or  input/output  values. 
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Overview  of  MULTI 

MULTI  is  an  interactive  program  which  enables  a  computer  user  to  design 
and  simulate  digital,  multivariable  control  laws  for  discrete,  linear, 
multiple-output  systems.  The  three  design  methods  incorporated  in  this 
program  were  developed  by  Professor  Brian  Porter  of  the  University  of 
Salford,  England.  To  help  the  user  become  quickly  familiar  with  the 
program,  the  following  overview  is  provided. 

MULTI  is  able  to  design  and  simulate  control  laws  for  three  types 
of  plants: 

(1)  Regular  -  The  linear,  multivariable  plant  dynamics  are 
described  by  the  usual  state  and  output  equations,  and  the  first  Markov 
parameter,  CB,  has  full  rank  (Ref.  2). 

(2)  Irregular  -  The  plant  dynamics  are  described  by  state  and 
output  equations,  and  the  first  Markov  parameter  does  not  have  full  rank. 
Thus,  the  system  must  be  augmented  by  a  measurement  matrix  so  that  a 
control  law  may  be  developed  (Ref.  3). 

(3)  Unknown  -  The  plant,  state  and  output  equations  and  unknown, 
but  the  steady  state  transfer  function  matrix,  G(0),  is  obtainable  from 
off-line  tests  if  the  open-loop  plant  is  stable.  This  method  is  also 
applicable  to  known  plants  (Ref.  1). 

All  control  law  designs  may  be  evaluated  by  a  discrete-time 
simulation  of  the  system,  actuator,  and  sensor  equations.  The  user  may 
elect  to  obtain  a  terminal  plot  or  a  Calcomp  plot  of  the  system  input 
and/or  output  responses. 

MULTI  focuses  on  Porter's  digital  design  as  opposed  to  a 
continuous  design  method.  At  various  points  during  the  design  and 
simulation  process,  the  program  can  provide  intermediate  information 
about  system  matrices  or  input/output  values. 


Input,  design  and  simulation  data  is  automatically  stored  on 
separate  memory  files  for  future  use.  This  prevents  the  time  consuming 
task  of  re-entering  data  when  the  same  system  is  under  study. 

The  program  performs  full  error  detection  and  diagnostics  when 
there  are  deficiencies  in  input  data,  design  data,  plotting  data,  or  if 
the  plant  equations  are  not  sufficient  for  the  type  of  design  being 
tried.  Input  error  detection  and  recovery  is  limited,  relying  mainly  on 
the  automatic  CYBER  error  detection  capability.  The  user  should  be  wary 
of  entering  values  inconsistent  with  the  limitations  described  in  this 
guide,  since  exceeding  these  limits  may  cause  unexpected  termination  of 
the  program.  If  this  occurs  the  memory  files  are  not  generated  and  all 
data  is  lost. 

It  must  be  stated  that  MULTI  is  not  as  dynamic  in  its  ability  to 
receive  input  data  as  other  programs  with  which  the  user  may  be  familiar 
(i.e.  TOTAL  [Ref.  13],  CESA  [Ref.  9]).  At  present,  character  inputs  to 
display  current  matrix  values  or  system  design  values  are  not  available 
except  by  requesting  the  proper  option  number.  Also  at  present,  entering 
a  symbol,  rather  than  actual  data,  to  abort  an  option,  is  not 
recognized  unless  specifically  noted  in  the  program  instructions  below. 

1 .  Introduction  to  MULTI 

The  MULTI  computer  package  contains  approximately  40  ordered 
options  which  give  the  user  an  interactive,  iterative  approach  in  the 
design  and  simulation  of  control  laws  for  linear,  multivariable  plants. 
The  control  law  assures  that  for  constant  commands  the  output  tracks  the 
input  and  that  disturbance  rejection  is  accomplished.  The  method  can  be 
used  only  if  (1)  the  introduction  of  integral  action  preserves 
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stabilizability ,  (2)  the  number  of  outputa  is  equal  to  or  leas  that  the 
number  of  inputs,  and  (3)  the  forward  transfer  function  has  no 
zero-valued  transmission  zeros. 

Once  the  model  is  entered  and  the  parameters  for  the  control  lav 
are  selected,  a  simulation  can  be  run  which  provides  a  full 
time-sequential  listing  of  input/output  values  and  a  tabular  listing  of 
the  input/output  values,  if  the  user  desires.  The  following  discussion 
should  give  the  user  all  the  information  for  optimal  use  of  the  computer 
program  package  called  MULTI. 

To  attach  and  run  MULTI,  after  LOGIN,  the  user  enters: 

COMMAND  -  CORRECT, INPUT, OUTPUT 
COMMAND  -  ATTACH , MULTI , ID-AFIT , SH-AFIT 
COMMAND  -  SCREEN, TULL 
COMMAND  -  (ATTACH  DATA  FILES) 

COMMAND  -  MULTI 

1.1  MULTI' a  Innut  Modes 

MULTI  has  three  input  modes  in  which  it  requests  input 
information  from  the  user.  These  modes  are  called  the  OPTION  Mode,  the 
DATA  Mode  and  the  QUERY  Mode.  Each  mode  has  its  own  restrictions  on 
allowable  input  and  its  own  method  of  requesting  information. 

OPTION  MODE  The  OPTION  Mode  is  the  executive  command  mode  for 
MULTI.  It  has  the  following  prompt  message: 

OPTIOM,  PLEASE > 

Once  the  program  is  in  this  mode,  the  user  is  allowed  to  Belect 
any  option  referring  to  input,  design,  simulation,  plotting,  or  program 
termination.  This  mode  is  the  main  input  mode  of  MULTI  since  it  allows 
the  user  complete  access  to  the  program  and  data  files. 
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PATA  Modt  MULTI  eater*  the  date  mode  vhen  information  i*  needed 
to  perform  an  option.  Each  data  input  request  has  a  specific  statement 
associated  with  it,  and  each  request  is  terminated  by  the  symbol, 
Normally,  data  input  must  be  numerical  and  separated  by  either  a  cornua,  a 
space,  or  the  return  key  of  the  terminal.  There  are  other  miscellaneous 
data  mode  options  which  require  character  entries.  All  replies  can  be 
accomplished  by  a  single  character  except  when  the  program  asks  the  user 
to  enter  the  choice  of  "INPUT"  of  "OUTPUT".  In  this  case  the  program 
requires  the  entire  word,  spelled  correctly,  to  continue  properly. 

QUERY  Mode  MULTI  uses  the  QUERY  Mode  an  an  input  checking  tool 
and  to  suppress  date  printing.  The  name  QUERY  Mode  suggests  that  the 
program  is  asking  the  user  if  the  data  is  correct  or  is  to  be  printed. 

The  messages  is  this  mode  are  either  of  the  form: 

KRTKR  "O"  TO  SKIP  DATA  PRINTOUT 

BRER  "1"  TO  OBTAIN  DATA  PRINTOUT... 

or 

IS  THIS  CORRECT... YES tN0,$... 

In  the  latter  case,  a  "NO"  reply  will  return  the  user  to  the  data 
input  point  while  a  "$"  reply  will  terminate  the  option.  Although  the 
option  is  terminated,  the  values  just  entered  are  placed  in  memory.  The 
user  can  also  enter  a  simple  "Y"  of  "N". 

1.2  MULTI 's  Options 

MULTI  contains  40  main  options  that  allow  the  user  versatility  in 
the  input,  design,  simulation  and  analysis  of  control  laws  for  unknown, 
regular  or  irregular  plants.  The  options  are  grouped  together  into  four 
major  option  blocks  as  follows: 
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BLOCS  1  OFTIOIS  #  0  thru  #  9  Plant  Input  Options 

BLOCS  2  OFTIOIS  #10  thru  #19  Design  Parameter  Options 

BLOCS  3  OFTIOIS  #20  thru  #29  Simulation  Options 

BLOCS  4  OFTIOIS  #30  thru  #39  Plotting  Options 

Option  BLOCK  5  is  composed  of  printing  commands  vhich  print  the  data 
entered  into  the  program  from  the  related  main  option.  That  is,  to  print 
the  data  which  was  entered  into  the  program  by  using  Option  #5,  the  user 
selects  Option  #103. 

Although  the  selectiion  of  a  sampling  t:me  is  actually  a  design 
parameter,  BLOCK  2  is  used  only  to  form  the  control  law  matrices  which 
are  independant  of  sampling  time.  Thus,  the  sampling  time  selection  is 
accomplished  in  the  simulation  option  block  where  the  actual  control  is 
used. 


2.  Complete  Description  of  MULTI 'a  Options 

In  order  to  utilize  MOLTI's  assets  fully,  it  is  necessary  to  have 
a  complete  understanding  of  the  40  main  options  available  in  MULTI.  This 
manual  is  intended  to  provide  all  of  the  information  needed  to  accomplish 
this  requirement.  Each  of  the  following  options  described  may  be 
selected  by  simply  typing  in  the  option  number  while  the  program  is  in 
:*he  option  mode. 

2.1  BLOCK  1  -  Flant  Input  Options 

Since  a  control  law  design  cannot  be  started  without  the 
knowledge  of  the  plant  model  or  steady  state  transfer  function,  it  is 
necesary  to  enter  some  representation  of  the  system  for  which  the  control 
law  is  to  be  designed.  Thus  BLOCK  1  is  an  integral  part  of  the  MULTI 
program.  Flant  data  may  be  entered  via  G(0)  matrix  (OPTION  #1)  or  by 
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entering  the  individual  matrices  of  the  state-space  representation  of  the 
plant  (OPTIONS  #2  and  #3). 


Although  the  steady  state  transfer  function,  G(0),  can  be  used  to 
obtain  control  law  matrices,  it  does  not  provide  sufficient  information 
to  run  a  simulation  and  obtain  plots  of  the  responses  by  using  option 
BLOCKS  3  and  4. 

The  state  and  output  equations  of  the  plant  follow  the  form: 
x  (kl  +  T)  -  Ax  (kl)  ♦  B  u  (kl)  (C-l) 

X  (kT)  -  C  x  (kl)  +  D  u  (kT)  (C-2) 

where 

A  ■  sampled -data  plant  matrix 
B  “  samp led -data  control  input  matrix 
C  “  syatem  output  matrix 
D  ■  feed-forward  matrix 

Note  that  there  is  not  provision  for  a  disturbance  matrix  in  Eqn. 

C-l. 

The  plant  input  option  block  also  allows  the  user  to  input 

actuator  and  sensor  state  equations  into  the  program.  These  equations 

are  of  the  same  format  as  shown  for  the  syBtem  state  equations,  however 

there  is  no  feed-forward  matrix,  D.  If  the  actuator  and  sensor  values 

are  not  entered  by  the  user,  they  are  set  equal  to  single-order  servos 
-4 

with  1  X  10  time  constants. 

The  program  accepts  a  maximum  of  ten  states,  ten  inputs  and  ten 
outputs.  The  program  only  allows  for  2nd-order  actuators  and  sensors, 
thus  they  may  have  to  be  approximated  by  a  lst-order  or  2nd-order 
equivalent  servo. 
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OPTION  »Q  -  List  Option*  0  thru  9 

Thia  option  liata  the  plant  input  optiona  from  0  to  9. 

OPTlOj  #1  -  Enter  G(0)  Matrix 

Thia  option  entera  the  steady  state  transfer  function  matrix, 
j(0).  The  user  ia  asked  to  supply  the  number  of  inputs,  M,  and  the 
number  of  outputs,  P,  thus  setting  the  dimension  of  the  matrix.  The  data 
is  then  entered  by  row  following  each  prompt  message. 

OPTIQM  #2  -  Enter  lumber  of  States.  Inputs,  and  Outputs 

This  option  asks  the  user  to  enter  the  number  of  states,  N,  the 
number  of  inputs,  M,  and  the  number  of  outputs,  P.  These  values  must  be 
entered  sequentially  as  N,  M,  and  P. 

OPTIOK  #3  -  Enter  A.  B.  C  4  D  Plant  Matrices 

This  option  enters  the  plant  A,  B,  C  and  D  matrices.  Each  input 
has  the  same  format  and  the  user  enters  the  data  by  rows  after  the 
prompt.  After  each  matrix  is  entered,  its  values  are  automatically 
echoed  back  to  the  user  for  checking.  After  the  A,  B  and  C  matrices  are 
entered,  the  program  asks  if  there  is  a  D  matrix.  If  the  reply  is 
negative,  the  option  is  terminated,  and  the  program  sets  the  D  matrix 
values  to  zero. 

OPTIOBS  #4  and  fS  -  Enter  actuator/Sensor  State  Equation  Matrix  Data 

These  options  enter  or  eliminate  the  actuator  and  sensor  state 
equation  A,  B,  and  C  matrices.  When  these  options  are  selected  the 
following  prompt  will  appear: 

OTTER  "0"  TO  ELIKL1ATE  ACTUATORS / SE1S0RS 
OTTER  "1"  TO  SET  ACTUATOR/ SKRSOR  VALUES... 

If  "1"  is  entered,  the  user  is  first  asked  to  enter  the  number  of 
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states  of  each  actuator/sensor.  The  input  sequence  must  correspond  to 
the  sequence  of  the  actuators/sensors.  The  user  then  inputs  the  A  matrix 
values  by  row s,  the  B  column  values  and  the  C  row  values  after  each 
prompt . 

OPTIOBS  #6  and  #7  -  Reserved  Options 
OPTIOB  #8  -  Copt  G(0)  Data  from  Local  File 

An  alternate  data  input  mode  from  OPTION  #1  is  to  use  OPTION  #8. 
The  data  may  be  contained  in  any  local  file  and  need  not  contain  *EOR  or 
♦EOF  statements.  The  file  should  not  be  a  permanent  file  called  MEMO  and 
if  it  is  a  local  file  by  that  name,  its  contents  will  be  overwritten  upon 
normal  program  termination.  The  file  must  be  in  the  proper  format  as 
shown  below: 

100  *  2  Indicates  G(0)  information 

110  “33  M,P  values 

120  “  -1 .  5.  0.  G(0)  matrix  values  by  row 

130  -  1.  2.  7. 

140  -  1.2  4.  .123 E01 

The  matrix  data  values  can  be  entered  in  any  format.  The  line 
numbers  above  indicate  that  the  data  file  may  be  created  in  the  CYBER 
"EDITOR",  however  the  line  numbers  must  be  suppressed  when  the  file  is 
saved  into  a  local  file. 

When  using  this  option,  the  program  asks  the  user  to  specify  and 
verify  the  name  of  a  local  file  which  holds  the  data.  The  data  is  then 
read  into  memory  and  the  user  can  verify  the  data  entries  by  using 
OPTIONS  #101  and  #102. 

QFTIOW  #9  ~  Copt  Plant.  Actuator  and  Sensor  Info  from  Local  File 

This  option  copies  plant  Btate  equations,  actuator  and  sensor 
information  from  a  local  data  file  into  computer  storage  location.  The 
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file  name  restrictions  mentioned  in  the  OPTION  #8  discussion  also  apply 
to  this  option.  This  option  is  the  alternate  to  OPTIONS  #2  through  #5. 
The  following  lines  show  the  exact  format  for  the  data  file: 


100-1 

Indicates  state  equation  information 

110  -  3  2  2 

>,M,P  values 

120  -  1.  5.5  7B-02 

A  matrix  values  by  row 

130  -  .2  -1.2  0 

140  -  .5  .005  99.1 

150  -  4.  0. 

B  matrix  values  by  row 

160  -  1.  1. 

170  -  0.  1. 

180  -  0.  1  1. 

C  matrix  values  by  row 

190  -  1.  1.  1. 

200  -1 

Indicates  no  D  matrix 

210  -  K 

Indicates  no  actuators 

220  -  Y 

Indicates  no  sensor  data  follows 

230  -  1  2 

States  sensor  #1  and  sensor  #2 

240  -  -9999. 

Sensor  #1  A  matrix  value 

250  -  9999. 

Sensor  #1  B  matrix  value 

260  -  1. 

Sensor  #1  C  matrix  value 

270  -  0.  1. 

Sensor  #2  A  matrix  valuess  by  row 

280  -  -888.063  -147 

.2 

290  -  0.  293.141 

Sensor  #2  B  column  values 

300  -  1.  0. 

Sensor  #2  C  row  values 

The  matrix  data  values  can  be  entered  in  any  format.  The  line 
number  above  indicate  that  the  data  file  may  be  created  in  the  CYBER 
"EDITOR",  however  the  line  numbers  must  by  suppressed  when  the  file  is 
saved  into  a  local  file.  The  *EOR  and  *EOF  entries  at  the  end  of  the 
listing  are  not  required 

When  using  OPTION  #9,  the  program  asks  the  user  to  specify  and 
verify  the  name  of  a  local  file  which  holds  the  data.  The  data  is  then 
read  into  memory  and  the  user  can  verify  the  data  entries,  if  desired,  by 
using  OPTIONS  #102  thru  #105. 

2.2  BLOCK  2  ~  Design  Input  Options 

MULTI's  control  law  development  is  based  upon  forming  the 
controller  matrix,  K,  from  proper  partitions  of  the  G(0)  or  A,  B,  C  and  D 
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matrices  of  the  plant.  In  turn,  KO  and  Kl  matrices  are  created  using  the 
scalar  variables  ALPHA  and  EPSILON.  The  equations  are: 
for  unknown  plants... 

u  (kT)  -  T  [  KO  *  e  (kT)  -  Kl  *  a  (kT)]  (C-3) 

for  regular  or  irregular  plants... 

u  (kT)  -  f  IKO  *  e  (kT)  ♦  Kl  *  *  (kT)]  (C-4) 

where 

KO  -  ALPHA  *  EPSILO*  *  K  (C-5) 

Kl  -  KPSILOB  *  K  (C-A) 

f  -  1  /  T  (C-7) 

In  the  above  formulas, 
u  ”  input  vector 
e  -  proportional  error  vector 
z  -  integral  error  vector 
T  »  sampling  period 

The  K  matrix  above  is  equivalent  to  G(0)  *  *  SIGMA,  where  SIGMA 
is  a  diagonal  output  weighting  matrix,  and 
for  unknown  plants... 

G(0)  •  B  -  [  C  *  A'1  *  B  ]  (C-8) 

for  regular  plants... 

G(0)  -  I  C  *  E  ]  (C-9) 

for  irregular  plants... 

G(0)  -  [P2  *  B2  ]  (C-10) 

In  the  irregular  case,  the  F 2  matrix  is  a  partition  of  the 
measurement  equation  matrix  which  is  formed  as  follows: 

w  -  I  Fj  ,  P2  ]  -  ICj  ♦  M  *  A  n  .  C2  ♦  M  *  A12  ]  (C-ll) 
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The  program  refers  to  the  matrix,  M,  ia  Eqn.  (C-ll)  aa  the 


measurement  matrix. 

OPTIC  #10  -  Liat  Options  10  thru  19 

This  option  lists  the  design  parameter  input  options  10  to  19. 
OPTIC  #11  -  Enter  ALPHA 

This  option  is  used  to  set  the  proportion  of  direct  to  integral 
feedback,  such  that: 

KO  -  ALPHA  *  K1  (C-12) 

OPTIC  #12  ~  Enter  SIGMA  Weighting  Matrix 

In  this  option  the  SIGMA  weighting  matrix  diagonal  elements  are 
entered.  This  matrix  is  used  in  forming  the  controller  matrix,  K,  and 
has  the  effect  of  weighting  the  different  inputs  and  outputs. 

OPTIC  #13  -  Enter  KPSILOH  (SIGMA  Weighting  Matrix  Multiplier) 

In  this  option,  EPSILON  is  entered.  This  scalar  parameter  is  the 
SIGMA  weighting  matrix  multiplier .  It  is  multiplied  to  each  diagonal 
element  in  the  SIGMA  matrix. 

OPTIC  #14  -  Run  Design . ■ .Unknown ,  Regular  A  Irregular  Plants 

In  this  option  the  KO  and  K1  matrices  are  formed.  The  user  is 
asked  to  enter  the  type  of  design  to  be  accomplished  and  the  reply  can  be 
either  the  full  word  or  a  single  character. 

When  using  an  unknown  design  with  state  equation  input  from 
OPTION  #3,  the  user  can  opt  to  print  out  the  G(0)  matrix  after  it  is 
formed. 

If  the  regular  design  is  used,  the  program  asks  if  the  user 
wishes  to  see  the  C  *  B  matrix  after  it  is  formed. 

The  first  step  in  using  the  irregular  design  is  to  provide  the 


program  vith  a  measurement  matrix.  If  the  measurement  matrix  is  formed 
in  OPTION  #18  or  previously  entered,  the  user  need  not  re-initialise  the 
values  unless  they  are  to  be  changed.  The  rows  of  the  matrix  are  entered 
after  the  prompt  is  given.  When  all  matrix  values  are  received  the 
program  enters  the  QUERY  Mode.  When  forming  the  KO  and  K1  matrices,  the 
user  has  the  option  of  printing  the  [  ,  F„  ]  matrix. 

When  the  matrices,  KO  and  K1 ,  are  formed,  MULTI  gives  the  user 
the  message: 

KO  &  K1  MATRICES  FORMED 

and  the  specific  matrix  values  can  be  checked  by  using  OPTION  #114. 

OPTIONS  #15  thru  #17  -  Reserved  Options 

OPTION  #18  -  Measurement  Matrix  Formation. . .or 

Bgg  k  Cglam  Operations  on  C  *  i  or  other  Matrix 

The  option  has  two  functions  as  noted  by  the  title. 

If  forming  a  measurement  matrix  to  be  used  in  the  irregular  plant 
controller  design  of  OPTION  #14,  the  user  must  enter  the  CESA  feedback 
gain  matrix,  K1 ,  from  CESA  Option  #38  (Ref.  9).  MULTI  advises  when  the 
measurement  matrix  has  been  formed.  It  can  be  printed  out  by  using 
OPTION  #118. 

This  option  can  also  be  used  to  perform  simple  row  and  column 
operations  on  any  matrix.  There  is  provision  for  the  user  to  perform 
these  operations  on  the  [  C  *  B  ]  matrix  without  entering  or  altering  the 
actual  matrix  values.  Otherwise  the  user  must  enter  the  matrix  size  and 
then  the  row  values  after  each  prompt.  The  matrix  of  interest  can  be  no 
larger  than  10  X  10.  The  a d d i t i o n  /  s  u  b  t  r  a  c  t  i  o n  and 
multiplication/division  operations  are  defined  as  follows: 


180 


Enter  X,  Y,  Z  such  that... 


(X,rov)  *  (Y, factor)  *  (Z,rov)  -  (lew  Z,rov) 

(X, column)  *  (Y, factor)  +  (Z.column)  -  (lew  Z, colon) 

(X.row)  *  (Y, factor)  “  (lev  X,row) 
or  (X,column)  *  (Y, factor)  -  (lev  X.coltaan) 

After  each  operation  the  resulting  matrix  is  printed  so  that  the 
next  operation  can  be  determined.  The  final  matrix  is  not  stored  into 
any  memory  location  for  future  use  in  other  options  and  is  destroyed  at 
the  termination  of  the  option. 

OPTIOK  #19  -  Copt  Design  Parameter  Data  from  Local  File 

This  option  is  uBed  as  an  alternative  to  OPTIONS  #11  thru  #14  and 
copies  design  parameter  data  from  a  local  data  file.  As  it  is  mandatory 
th8t  the  data  file  be  in  the  proper  format,  the  following  example  file  is 
provided : 

100-0  or  1  or  I  Indicates  Unknown,  Regular,  Irregular 

110-  2  ALPHA  value 

120-  .111  SIGMA  matrix  diagonal  elements 

130-  .5  EPSIL0M  value 

140-  4.6  4.6  4.6  KO  matrix  values  by  row 

150-  -4.6  -4.6  -4.6 

160-  1.  1.  6. 

170-  2.3  2.3  2.3  XI  matrix  values  by  row 

180-  -2.3  -2.3  -2.3 
190-  .5  .5  3. 

200-  .25  Measurement  matrix  values  by  row 

210-  0. 

220-  0. 

If  the  ALPHA  value  is  1,  the  file  must  not  contain  K1  matrix 
values  since  K1  is  identical  to  KO.  The  measurement  matrix  values  are 
only  read  if  the  file  contains  data  for  am  irregular  plant  design.  The 
line  numbers  above  indicate  only  that  the  data  file  may  be  created  in  the 
CYBER  "EDITOR".  However,  the  line  numbers  must  not  appear  in  the  actual 
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data  file.  The  file  terminators,  *E0R  and  *EOF,  are  optional. 

Pror  to  using  OPTION  #19,  either  OPTION  #2  or  OPTION  #9  must  he 
used  to  set  the  values  for  the  number  of  inputs  and  outputs.  To  use 
OPTION  #19,  the  user  must  supply  and  verify  the  name  of  the  local  file 
which  holds  the  design  parameter  data.  The  program  reads  the  data  into 
the  proper  memory  locations  and  the  user  can  verify  the  values  by 
choosing  OPTIONS  #111  to  #114.  Refer  to  OPTION  #8  for  a  complete 
description  of  the  restrictions  that  apply  to  local  data  file  names. 

2.3  BLOCK  3  -  Simulation  Options 

The  simulation  options  of  this  block  are  used  to  evaluate  the 
controller  matrices  designed  in  OPTION  #14.  The  simulation  is  performed 
after  the  user  provides  the  initial  values  of  the  states  and  integrators, 
the  command  input  vector  ,  and  the  necessary  time  parameters.  The  actual 
simulation  run  is  accomplished  in  OPTION  #26.  The  simulation  iB  obtained 
by  numerically  solving  a  set  of  ordinary  differential  equations  formed 
from  the  plant,  actuator  and  sensor  state  equation  matrices  and  is  run 
from  time  zero. 

The  ASD  Library  Subroutine  ODE  (Ref.  20)  is  the  basis  for  the 
simulation.  The  user  should  be  familiar  with  this  subroutine,  or  as  a 
minimum,  have  access  to  ODE's  error  descriptions. 

During  the  simulation  the  control  input  is  held  constant  over 
each  sampling  period.  Since  it  may  be  advantageous  to  observe  the  system 
output  between  sampling  periods,  the  user  can  specify  a  step  time  less 
than  the  sampling  time. 

As  the  simulation  is  run,  a  time-sequential  printout  of  input  and 
output  values  for  every  time  interval  of  the  simulation  run  can  be 
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•elected.  This  printout  includes  actual  input  and  controlled  input 
values  if  the  user  has  limited  the  inputs  to  certain  values.  If 
limitations  are  applied  to  the  input  values,  the  program  generates  a 
tabular  listing  of  how  often  the  limits  are  exceeded  during  the 
simulation. 

The  simulation  data  is  then  stored  in  a  matrix  for  later  printout 
or  for  plotting. 

OPTIOB  #20  ~  List  Options  20  thru  29 

This  option  lists  the  simulation  options  from  20  to  29. 

OPTIOH  #21  -  Set  State  &  Integrator  Initial  Values 

In  this  option  the  user  specifies  the  state  and  integrator  values 
at  time  T(0).  At  present,  all  simulations  must  begin  at  zero  time  and 
proceed  forward.  The  inital  condition  values  for  the  states  are  entered 
after  the  first  prompt  message,  and  the  integrator  values  are  entered 
after  the  second  prompt  message.  The  values  for  each  must  be 
sequential. 

OPTIOB  #22  -  Set  Irnrnt  Command  Vector.  ¥ 

This  option  is  used  to  select  the  specific  magnitude  for  each  of 
the  commanded  inputs.  At  present,  all  inputs  are  step  values.  The 
entries  for  the  input  command  vector  are  to  be  made  sequentially. 

OPTIOB  #23  -  Enter  Sampling  Times 

The  sampling  times  are  entered  in  this  option.  There  may  be  up 


to  two  different  sampling  times  entered.  The  program  first  askB  how  many 
sampling  times  will  be  entered,  and  then  prompts  the  user  to  enter  the 
same  number  of  sampling  time  values. 


OPTION  #24  is  used  to  enter  the  time  length  for  the  simulation 
run.  The  user  is  cautioned  to  select  this  time  vith  care  since  the 
combination  of  simulation  time/sample  time  or  simulation  time/step  time 
determines  the  solution  interval  time.  If  the  simulation  time  is  too 
large,  Subroutine  ODE  may  not  be  able  to  finish  a  simulation  run  before 
the  CYBER  CP  time  limit.  In  this  case,  MULTI  will  be  terminated  without 
warning  and  no  memory  files  will  be  generated. 

OPTION  #25  -  Enter  Calculation  8ten  Size.  ST 

As  mentioned  in  the  introduction  to  this  section,  the  system 
input  or  output  values  between  sampling  times  may  be  of  interest.  This 
option  allows  the  user  to  select  the  step  time  between  each  sampling 
period.  If  the  step  time  is  chosen  larger  than  the  sampling  time,  the 
program  proceeds  using  the  sampling  time  as  the  discrete  time  interval. 
For  example,  if  sampling  time  is  0.02  seconds,  a  step  size  of  0.01  would 
give  two  sets  of  data  during  the  sampling  interval.  If  the  step  size  was 
greater  than  .02  seconds,  it  is  set  to  sampling  time  automatically. 

OPTION  #26  "  to"  Simlitinn 

When  OPTION  #26  is  chosen,  MULTI  forms  and  solves  the  set  of 
ordinary  differential  equations  formed  from  the  plant,  actuator  and 
sensor  state  equations.  Very  fast  actuators  and  sensors  are  approximated 
if  no  actual  values  have  been  entered  and  the  user  receives  messages  to 
indicate  that  no  values  have  been  set.  The  user  also  receives  a  message 
that  no  control  limits  have  been  applied  to  the  inputs,  if  such  is  the 


MULTI  runs  a  simulation  for  each  sampling  time  selected  in  OPTION 


#23.  As  the  program  proceeds  through  each  sampling  time,  the  uaer  is 
first  given  the  length  of  the  run,  the  step  time  and  the  sampling  time. 
The  user  is  then  told  how  many  time  increments  are  generated  during  the 
simulation  run.  The  user  can  opt  to  suppress  the  time-sequential  listing 
of  input/output  values.  The  mesage: 

— CALCULATIONS  Dl  PK0G1KSS — 

ib  provided  if  the  time-sequential  listing  is  suppressed,  since  the 
solution,  in  some  cases,  may  take  several  seconds. 

If  the  number  of  time  increments  over  the  simulation  time 
interval  is  greater  than  100,  the  input/output  data  must  be  "packed"  into 
a  matrix  for  later  use  in  plotting  or  providing  a  tabular  listing  of  the 
input/output  values.  When  this  occurs,  missing  time  increments  are  noted 
when  the  tabular  data  of  input/output  values  are  listed.  However,  the 
user  can  always  obtain  a  full  list  of  these  values  by  choosing  the 
time-sequential  listing. 

At  the  conclusion  of  each  simulation  run,  the  program  provides  a 
table  telling  show  many  times  the  input  values  exceeded  the  limitations 
chosen  in  OPTION  #27.  This  listing  is  not  retained  in  memory  and  can 
only  be  accessed  once. 

The  tabular  listing  of  the  input/output  data  points  for  each 
sampling  time  can  be  obtained  from  OPTION  #126. 

OFTIOW  #27  -  8et  Control  *"p"t 

In  this  option  the  user  can  set  input  control  limitations.  The 
program  supplies  the  correct  directive  on  how  to  enter  the  limits  and 
provides  proper  messages  when  the  limitations  are  to  be  entered.  The 
user  tell  which  input  is  to  limited  and  then  gives  the  minimum  and 
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maximum  values  of  tbe  limitation. 

This  option  can  also  be  used  to  eliminate  control  input 
limitations  if  they  are  previously  entered. 

To  obtain  a  listing  of  tbe  current  input  limits.  Option  #127  is 


used. 

OPTIOM  #28  -  Reserved  Option 

OPTIOB  #29  -  Copy  Simulation  Parameters  for  a  Local  File 

This  single  option  is  used  to  enter  state  and  integrator  initial 
values,  command  vector  values,  sampling  times,  simulation  run  length  and 
interval  step  time.  It  can  be  used  as  an  alternate  to  OPTIONS  #21 
through  #25.  Tbe  reader  should  refer  to  the  discussion  of  OPTION  #8  for 
file  name  restrictions. 

The  data  file  must  follow  the  format  shown  below: 


Initial  state  values 

Initial  integrator  values 

Cosaand  vector  values  (transposed) 

#  of  sampling  times,  sampling  times 
Simulation  run  time 
Step  ti 


100-  0.  0.  0. 

110-  0.  0. 

120-  1.  -1. 

130-  2  .01  .02 
140-  3 
150-  .01 

As  with  the  other  data  file  examples,  tbe  line  numbers  above  only 
indicate  that  the  data  matrix  may  be  created  in  the  CYBER  "EDITOR",  and 
must  be  suppressed  when  the  file  is  Baved  as  a  local  file.  The  *E0R  and 
♦EOF  terminators  are  optional. 

Prior  to  using  this  option,  the  number  of  states,  inputs  and 
outputs  must  be  set  by  using  OPTION  #2  or  OPTION  #9.  To  use  OPTION  #29 
the  user  must  enter  and  verify  the  name  of  tbe  local  file  vhich  holds  tbe 
simulation  data.  After  the  program  reads  the  data  into  memory,  the  user 
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can  verify  the  entries  by  using  OPTIONS  #21  to  #25. 

2.4  BLOCK  4  -  Plotting  Options 

An  integral  aspect  of  any  design  is  the  ability  to  quickly 
interpret  and  analyse  the  results.  The  quickest  method  of  analyzing 
results  is  to  view  them  in  either  tabular  or  graphical  form.  MULTI 
provides  the  user  with  two  types  of  graphical  interpretation  tools. 

The  first  and  fastest  tool  is  a  quick  sketch  at  the  user's 
terminal  of  any  combination  of  input/output  responses.  The  second  tool 
is  a  CALCOMP  plot.  Each  has  its  own  advantages  and  limitations.  The  use 
of  both  is  explained  in  this  section. 

...  To  provide  a  plot  MULTI  requires  several  questions  to  be  answered 
concerning  what  data  is  to  be  plotted.  After  using  OPTION  #31  and  OPTION 
#34,  the  user  becomes  aware  that  an  easier  method  of  requesting  plots  is 
required.  In  this  light,  the  short  version  of  requesting  plots  is 
available  by  using  OPTIONS  #32  and  #35.  Finally,  if  the  user  desires  the 
same  information  to  be  plotted  after  another  simulation  (for  example: 
OUTPUT  #2  vs.  INPUT  #2)  and  has  already  entered  these  choices  by  using 
OPTION  #31,  #32,  #34  or  #35,  an  immediate  plot  can  be  generated  by  using 
OPTION  #33  or  OPTION  #36. 

When  OPTION  #26  has  been  completed,  MULTI  has  formed  and  stored  a 
three-dimensional  matrix  of  the  data  values.  For  this  discusson  it  is 
only  necessary  to  explain  that  the  row  dimension  is  equivalent  to  the 
input  or  output  values  at  each  time  increment,  the  column  dimension  is 
equivalent  to  the  input  or  output  number,  and  the  third  dimension  is 
equivalent  to  the  sampling  time.  A  plot  is  generated  after  a  user 
chooses  the  type  of  plot  to  be  generated,  the  sampling  time  of  interest. 
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input  or  output,  the  number  of  inputs  or  outputs  and  the  particular 
inputs  or  outputs  of  interest.  After  this  information  is  received,  MULTI 
can  proceed  with  plotting  the  correct  data  matrix  rows  and  columns. 

There  are  four  basic  types  of  plots  that  can  be  obtained.  For  a 
single  sampling  time,  the  user  can  select  a  plot  of  any  output  versus  any 
input,  a  plot  of  one  to  four  inputs  or  outputs  on  a  single  graph,  or  a 
plot  of  the  same  input  or  output  for  up  to  four  different  flight 
conditions.  For  the  latter  type  of  plot,  the  user  returns  to  the  main 
program  to  run  another  simulation  and  then  re-enters  the  plotting 
routine.  The  last  type  of  plot  is  for  comparing  inputs  and  outputs 
between  different  sampling  times.  When  planning  to  plot  data  from  two 
different  sampling  times,  it  is  mandatory  that  the  step  time,  for  the 
simulation  runs,  be  equal  to  or  less  than  the  smallest  sampling  time.  If 
this  restriction  is  not  followed,  MULTI  is  terminated  early  and  no  data 
files  are  generated  by  the  program. 

OPTIONS  #31  through  #33  pertain  to  terminal  plots,  or  quick 
sketches  at  the  user's  terminal.  OPTIONS  #34  through  #36  are  for  CALCOMP 
plots.  As  mentioned  above,  there  are  shortened  versions  of  both  the 
terminal  plot  and  CALCOMP  plot  requests.  The  same  plots  are  obtained 
with  these  shortened  versions  as  are  produced  by  using  the  full  versions. 
The  difference  lies  in  entering  the  data  for  the  plot. 

QPTIOB  #30  -  List  Options  30  thru  39 

This  option  lists  the  simulation  options  from  30  to  39. 

QPTIOB  #31  -  Quick  8ketch  at  User's  Terminal 

This  option  provides  the  user  with  4  plotting  choices.  The 
choices  are  presented  in  a  menu  at  the  start  of  the  option.  Each  choice 
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has  a  set  of  questions  associated  with  it.  For  example,  choice  #2 
requires  the  user  to  enter  the  sampling  time  number,  the  number  of  inputs 
and  outputs,  and  the  specific  inputs  or  outputs  to  be  plotted.  These 
questions  are  presented  and  need  to  be  answered  each  time  this  option  is 
used.  Since  this  can  be  a  tedious  task  for  more  that  a  few  plots,  the 
user  may  choose  the  shorter  versions  of  this  option,  OPTION  #32  or  OPTION 
#33,  when  applicable. 

OPTION  #32  -  Quick  Sketch  -  Short  Version 

If  the  user  is  experienced  in  the  use  of  MULTI,  he  may  elect  this 
option  rather  than  OPTION  #31.  The  required  plotting  data  is  entered  in 
one  data  string  as  follows; 

TYPE  #1  -  Input/Output  Pairs 

ENTER...#  of  pairs  (OPTION  #35  only),  sampling  time  #,  input  #s, 
output  #6 

TYPE  #2  -  Inputs  or  Outputs 

ENTER. . .sampling  time  #,  #  of  inputs/outputs,  input/output  #s 
TYPE  #3  -  Multiple  Plight  Conditions 

ENTER...#  of  flight  conditions,  sampling  time  #,  input/output  # 

TYPE  #4  -  Multiple  Sampling  Times 

ENTER...#  of  sampling  times,  sampling  time  #s,  input/output  # 

OPTION  #33  -  Quick  Sketch  -  Retaining  Same  Plotting  Choices 

This  option  is  self-explanatory.  Once  the  information  for  a  plot 
has  been  entered  by  using  OPTIONS  #31  or  #32,  the  user  may  change  a 
design  parameter,  run  a  simulation,  and  choose  OPTION  #33  to  obtain  a 
plot.  This  plot  of  the  new  simulation  data  is  a  plot  of  the  same  type  as 
previously  generated. 
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OPTION  #34  -  Calcomp  Plot 


This  option  provides  the  user  with  the  same  4  plotting  choices  as 
in  OPTION  #31.  The  menu  presented  at  the  start  of  the  option  is  the 
same,  along  with  the  questions  associated  with  each  menu  choice.  These 
questions  must  be  answered  each  time  this  option  is  chosen  and  can  become 
a  tedious  task.  The  experienced  MULTI  user  may  choose  the  shorter 
version,  OPTION  #35,  or  when  applicable,  OPTION  #36.  The  program 
requests  two  titles  from  the  user.  The  first  is  a  30  character  maximum 
string  that  is  put  on  the  Y-axis  of  the  plot.  The  second  is  a  60 
character  maximum  strint  that  is  put  below  the  plot  as  a  full  title.  The 
program  then  requests  a  plot  size  factor.  Nominally,  this  is  set  at  0.9 
for  a  Gx9  plot. 

OPTION  #35  ~  Calcomp  Plot  -  Short  Version 

This  options  is  essentially  the  same  as  OPTION  #32  but  is  used 
for  CALCOMP  plots.  The  plotting  data  is  entered  as  a  data  string  with 
the  format  listed  in  OPTION  #32' s  description. 

OPTION  #36  -  Calcomp  plot  -  Retaining  Same  Plotting  Choices 

This  option  is  self-explanatory.  Once  the  information  for  a  plot 
has  been  entered  by  using  OPTIONS  #34  or  #35,  the  user  m*y  change  e 
design  parameter,  run  a  new  simulation,  and  choose  this  option  to  obtain 
a  plot.  This  plot  of  the  simulation  data  is  a  plot  of  the  same  type  as 
previously  generated  by  OPTION  #34  or  OPTION  #35. 

The  CALCOMP  plots  are  Btored  in  a  local  file  called  PLOT  which 
the  user  can  route  to  the  plotter  after  MULTI  is  terminated.  Local 
restrictions  permit  no  more  than  five  plots  on  one  plotting  file,  thus, 
MULTI  advises  the  user  when  five  plots  have  been  generated.  All  CALCOMP 
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plots  are  5  by  8  inches  in  size.  The  user  is  also  required  to  provide  a 
title  £or  the  plot.  This  title  is  placed  along  the  Y-axis  of  the  plot; 
the  X-axis  is  always  labeled  "TIME".  Titles  can  be  a  maximum  of  20 
characters  in  length. 

2.5  BLOCK  5  -  Printing  Options 

All  available  options  in  this  block  print  data  which  the  user  has 
entered  using  the  options  in  BLOCKS  1  to  4. 

OPTION  #100  -  List  Options  100  thru  130 

This  option  lists  the  printing  options  from  100  to  130. 

OPTIOHS  #101  to  #129  -  Print  Data  Entered  in  OPTIONS  #1  to  #29 

As  explained  in  Section  1.2,  the  options  in  this  range  are 
related  to  the  sub-100  numbered  options.  To  print  out  the  data  entered 
by  OPTION  #1,  the  user  selects  OPTION  #101;  to  print  out  the  data  entered 
by  OPTION  #2,  the  user  selects  OPTION  #102;  etc. 

OPTION  #103  prints  out  the  data  from  OPTION  #2  and  then  giveB  the 
user  the  option  of  selecting  the  A,  B,  C  or  D  matrix  from  OPTION  #3. 

ALPHA  and  EPSILON  are  combined  in  one  printout  and  obtained  by 
selecting  OPTIONS  #111  or  #113. 

OPTION  #114  prints  out  the  control  matrices,  K0  and  K1 ,  from 
OPTION  #14  and  also  prints  out  a  heading  explaining  which  type  of  plant 
the  user  chose  in  generating  the  matrices.  If  this  heading  is 
suppressed,  the  program  is  unsure  of  how  the  matrices  are  generated,  but 
the  values  displayed  are  still  correct  KO  and  K1  matrices. 

The  date  from  OPTIONS  #21  and  #22  are  also  combined  into  one 
printout  and  the  user  obtains  both  sets  cf  data  when  selecting  OPTION 
#121  or  OPTION  #122. 
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The  seme  is  true  for  templing  timet,  timulation  time  end 
calculation!  atep  site.  They  form  one  printout  and  can  be  printed  by 
eelecting  OPTION  #123,  #124  or  #125. 

When  OPTION  #126  it  entered,  the  uaer  is  required  to  choose  the 
tabular  listing  for  the  input  or  the  output.  The  full  word,  INPUT  or 
OUTPUT,  must  be  entered  and  spelled  correctly  for  the  program  to  continue 
properly. 

OPTION  #130  -  List  Plotting  Selections 

This  single  option  is  used  to  determine  the  current  plotting 
selections  from  the  last  entry  to  OPTIONS  #31  to  #36.  When  chosen,  the 
program  tells  the  user  if  the  data  is  for  a  terminal  plot  or  a  CALCOMP 
plot  and  then  gives  the  current  plotting  choices. 

2.6  MULTI  Generated  Memory  Files 

When  OPTION  #99  is  used  to  end  MULTI,  there  are  three  data  files 
generated.  A  file  called  MEMO  is  formed  from  the  data  entered  from 
option  BLOCK  1;  a  file  called  MEM10  is  created  from  the  data  entered  from 
option  BLOCK  2;  and  a  file  called  MEM20  is  opened  to  hold  the  data 
entered  in  option  BLOCK  3.  These  files  are  local  files  and  should  be 
rewound  before  they  are  Btore  into  permanent  file  space,  edited,  or 
routed  to  a  line  printer.  Another  filed  called  PLOT  is  also  generated  if 
the  user  selects  any  of  the  CALCOMP  plotting  options. 

The  routing  command  for  MEMO,  MEM10  and  MEM20  is: 

COMMAND  -  ROUTE ,  MEMO/MENl 0/MEM2 0 . DC-PB , TID-9 1 , FID-XXXXX, ST-CSB 
and  the  routing  command  for  PLOT  is: 

COMMAND  -  ROUTE, PLOT, DC-PT, TID-9 1  ,FID-XHXX,  ST-CSB 

These  routing  commands  send  the  files  to  the  facilities  at  AFIT. 
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The  user  can  choose  any  5  letter  "flag  ID"  designated  as  "XXXXX"  above. 

3.  gsaaasi  si  MBSSHm  optioms 

MULTI  contains  four  main  options  BLOCKS  which  have  ten  options 
each,  and  30  printing  options.  The  MULTI  options  are  sumnarized  below. 


BLOCK  1  -  Plant  Input  Options 


Option  #  0. 
Option  #  1. 
Option  #  2. 
Option  #  3. 
Option  #  4. 
Option  #  5. 
Option  #  6. 
Option  #  7 
Option  #  8. 
Option  #  9. 


List  Options  0  thru  9 
Enter  G(0)  Matrix 

Enter  #  of  8tatea,  Inputs  4  Outputs  (>,M,P) 

Enter  Plant  1,  1,  C  1  D  Matrices 
Enter  Actuator  State  Equation  Matrix  Oats 
Enter  Sensor  State  Equation  Matrix  Data 
Option  Reserved 
Option  Reserved 

Copy  G(0)  Info  from  Local  File 

Copy  Plant,  Actuator  &  Sensor  Info  from  Local  File 


BLOCK  2  -  Design  Parameter  Input  Options 


Option  #10. 
Option  #11. 
Option  #12. 
Option  #13. 
Option  #14. 
Option  #15. 
Option  #16. 
Option  #17. 
Option  #18. 

Option  #19. 


List  Options  10  thru  19 
Enter  ALPHA 

Enter  SIGMA  Weighting  Matrix 

Enter  EPSILOM  (SIGMA  Matrix  Multiplier) 

Run  Deaign... Unknown,  Regular  &  Irregular  Plants 
Option  Reserved 
Option  Reserved 
Option  Reserved 

Measurement  Matrix  Formation. . .or 
Row  &  Column  Operations  on  C*B  or  other  Matrix 
Copy  Deaign  Parameters  from  Local  File 


BLOCK  3  -  Simulation  Options 


Option  #20. 
Option  #21. 
Option  #22. 
Option  #23. 
Option  #24. 
Option  #25. 
Option  #26. 
Option  #27. 
Option  #28. 
Option  #29. 


List  Options  20  thru  29 

Set  State  and  Integrator  Initial  Values,  X(0)  6  Z(0) 

Set  Input  Comand  Vector  V 

Enter  Sample  Times 

Enter  Simulation  Time 

Enter  Calculation  Step  Size 

Run  Simulation 

Set  Control  Input  Limits 

Option  Reserved 

Copy  Simulation  Parameter  from  a  Local  File 


BLOCK  4  -  Plotting  Options 


Option  #30.  List  Opitons  30  thru  30 
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Option  #31.  Quick  Sketch  at  User's  Terminal 

Option  #32.  Quick  Sketch  -  Short  Version 

Option  #33.  Quick  Sketch  -  Retaining  Same  Plotting  Choices 

Option  #34.  CALCOMP  Plot 

Option  #35.  CALCOMP  Plot  -  Short  Version 

Option  #36.  CALCOMP  Plot  -  Retaining  Same  Plotting  Choices 

Option  #37.  Option  Reserved 

Option  #38.  Option  Reserved 

Option  #39.  Option  Reserved 

BLOCK  5  -  Printing  Options 

All  100*-Series  Options  Print  Data  Values... 

For  Values  Set  in  Option  #1...  Use  Option  #101 
For  Values  Set  in  Option  #2...  Use  Option  #102 

Etc. 

For  Plotting  Selections........ Use  Option  #130 

To  obtain  a  controller  design  and  simulation  these  option 
numbers,  in  the  following  order,  are  generally  used: 

For  Unknown  Plants... 

1  (or  9);  11,  12,  13,  14U 

or  243  (or  9);  11  to  14U  (or  19);  21  to  25  (or  29),  26;  31  to  36 
For  Irregular  Plants... 

263  (or  9);  11  to  141  (or  19);  21  to  25  (or  29),  26;  31  to  36 

If  a  measurement  matrix  is  to  be  generated  for  the  irregular 
case,  OPTION  #18  should  be  used  prior  to  OPTION  #14.  In  all  cases, 
OPTIONS  #4,  #5  and  #27  which  set  specific  actuator,  sensor  and  input 
control  limits  should  be  used,  if  desired,  prior  to  OPTION  #26. 

4.  nl»  for  Unknown  Plants  (Ref.  1) 

This  appendix  is  concluded  with  a  computer  generated  design  and 
simulation  using  an  unknown  plant  design  defined  by  the  following  state 
equations : 
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This  design  includes  no  actuators,  no  sensors  an'',  no  inpu 


control  limitations. 
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COMMAND-  ATTACH* MULTI * ID=T3201 09 
PFN  IS 
MULTI 

AT  CY=  001  SM=AFIT 
COHMAN [»-  C  ONMECT  >  I NPUT  >  OUTPUT 
C  ClMMANP-  SCREEN* 80 
COMMAND-  MULTI 


WELCOME  TO  MULTI VAR I  ABLE  DESIGN 
C  1981  PORTER* SMYTH >PASCHALL 

THIS  PROGRAM  USES  THE  DESIGN  TECHNIQUES  DES'ELOPED  BY 
PROFESSOR  BRIAN  PORTER*  UNIV.  OF  SALFORD*  ENGLAND 


OPTION,  PLEASE  >  «0 

C'LANT  INPUT  OPTIONS: 

0.  LIST  OPTIONS  0  THRU  9 
1.  ENTER  G<T0>  MATRIX 

3.  ENTER  «  OF  STATES*  INPUTS  &  OUTPUTS  <N*M*P) 

3.  ENTER  PLANT  A*  B>  C*  &  D  MATRICES 

4.  ENTER  ACTUATOR  STATE  EQUATION  MATRIX  DATA 

5.  ENTER  SENSOR  STATE  EQUATION  MATRIX  DATA  ' 

6.  OPTION  RESERVED 

7.  OPTION  RESERVED. 

S.  COPY  GO:o  INFO  FROM  LOCAL  FILE 

9.  COPY  PLANT*  ACTUATOR  8,  SENSOR  INFO  FROM  LDCAL  FILE 


OPTION*  PLEASE  >  «2 

THIS  OPTION  SETS  THE  NUMBER  OF  STATES*  INPUTS  8,  OUTPUTS 
ENTER  NUMBER  OF  STATES,  INPUTS  AND  OUTPUTS  >3,2*2 


OPTION,  PLEASE  >  «3 

this  option  enters  the  plant  a,b*c,  and  d  matrices 

ENTER  "A"  MATRIX... 3  ROWS  WITH  3  ELEMENTS  EACH 
ROW  1  >-l  0  0 
ROW  2  >0  -2  0 
ROW  3  >0  0  -3 

PLANT  MATRIX  A... 

1 QQ0E+01  0.  0. 

0.  -. 2000E+01  0. 

0.  0.  3000E+01 
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IS  THIS  CORRECT. . . YES*NO*$. . .  >YES 


ENTER  "£:**  MA TR I X ...  3  ROWS  WITH  2  ELEMENTS  EACH 
POM  1  >1  -1 
ROW  2  >1  l 
ROW  3  >1  0 

PLANT  MATRIX  B... 

.  100CE+Q1  -. 1 OOOE+Ol 

. 1 000E+01  . 1000E+01 

. 1000E+01  0. 


IS  THIS  CORRECT... YES» NO, *...>Y 

ENTER  “C"  MATRIX... 2  ROWS  WITH  3  ELEMENTS  EACH 
ROW  1  >1  0  0 
ROW  2  >0  1  1 

OUTPUT  MATRIX  C.  .  . 

. 1000E+01  0.  0. 

0.  . 1 G00E+01  . 1 00QE+01 


IS  THIS  CORREC T ...  YES , NO j  * ... > NO 

ENTER  "C"  MATRIX... 2  ROWS  WITH  3  ELEMENTS  EACH 
ROW  1  >1  0  0 
ROW  2  >0  1  0 

OUTPUT  MATRIX  C. . . 

. 1000E+01  0.  0. 

0.  . 100QE+01  0. 


IS  THI  S  CORRECT.  .  .  YES,  NO,  :f .  .  .  >  YES 
IS  THERE  A  “I r  MATRIX. . .YES'  OR  NO...>NO 


OPTION*  PLEASE  >  «103 

THERE  ARE  3  STATES*  2  INPUTS  ?<  2  OUTPUTS 
ENTER. . .A,B,C»D...FDR  PRINTOUT  >A 

PLANT  MATRIX  A... 

-. 1 OOOE+Ol  0.  0. 

0.  -. 2000E+01  0. 

0.  0.  3000E+01 


OPTION*  PLEASE  >  «103 

THERE  ARE  3  STATES »  2  INPUTS  t.  2  OUTPUTS 
ENTFP. . .  A, B» C» D. . • FOR  PRINTOUT  >B 

PLANT  MATRIX  B. . . 

. 1 OOOE+Gl  1 000E+01 

. I 000E+01  . 1 000E+01 

. 1000E+01  0. 


OPTION,  PLEASE  >  «103 

THERE  ARE  3  STATES,  2  INPUTS  *  2  OUTPUTS 
ENTER.  .  .A,B,C,  D.  ...FDR  PRINTOUT  >C 

OUTPUT  MATRIX  C... 

. 1000E+0I  0.  0. 

0.  . 1 000E+01  0. 


OPTION,  PLEASE  >  «:  10 

DESIGN  PARAMETER  INPUT  *  OPTIONS: 

10.  LIST  OPTIONS  10  THRU  13 

11.  ENTER  ALPHA 

12.  ENTER  SIGMA  WEIGHTING  MATRIX 

13.  ENTER  EPSILON  OSIGMA  MATRIX  MULTIPLIER) 

14.  RUN  DESIGN. . .UNKNOWN,  REGULAR  &  IRREGULAR  PLANTS 

15.  OPTION  RESERVED 
IE.  OPTION  RESERVED 
17.  OPTION  RESERVED 

IS.  MEASUREMENT  MATRIX  FORMATION. '..  OR 

ROW  &  COLUMN  OPERATIONS  ON  C>B  DR  OTHER  MATRIX 
19.  COPY  DESIGN  PARAMETERS  FPOM  LOCAL  FILE 


OPTION,  PLEASE  >  «11 

THIS  OPTION  SETS  THE  PROPORTION  OF  INTEGRAL  AND  DIRECT  FEEDBACK 
ENTER  ALPHA  >0.9 


OPTION,  PLEASE  >  »12 

THIS  OPTION  SETS  THE  WEIGHTING  BETWEEN  OUTPUT  CHANNELS 

ENTER  SI6MP  WEIGHTING  M ATP IX... 2  DIAGONAL  ELEMENTS  ONLY 
DIAGONAL  ELEMENTS  >2,1 


198 


OPTION.  PLEASE  >  «13 

THi:  OPTION  SETS  THE  WEIGHTING  MATRIX  MULTIPLIER 
ENTER  SIGMA  MATRIX  SCALAR  MULTIPLIER,  EPSILON  >60 


OPTION*  PLEASE  >  «14 

ENTER  DESIGN  METHOD. . . UNKNOWN  *  REGULAR  *  IRREGULAR. . . > UNKNOWN 

THIS  OPTION  COMPUTES  KO  &  K1  FOR  UNKNOWN  PLANTS 

ENTER  ”0"  TO  SKIP  G<0>  MATRIX  PRINTOUT 
ENTER  "1"  TD  OBTAIN  THIS  DATA  PRINTOUT. .. >1 

G<0>  MATRIX... 

. 1 000E+01  -. 1 000E+0I 

. 5000E+00  . 5000E+00 


K'O  5=  K1  MATRICES  FORMED 


OPTION,  PLEASE  >  «114 

CONTROL  MATRICES  ARE  FOR  PLANTS  WHICH  ARE  UNKNOWN 

KO  MATRIX... 

. 54  00E+02  . 5400E+02 

-.340 0E+02  . 54  0  0E+ 02 


K 1  MATRIX. . . 

. 6000E+02  . 6000E+02 

-.6000E+02  . 6000E+02 


OPTION,  PLEASE  >  «20 

SIMULATION  OPTIONS: 

20.  LIST  OPTIONS  20  THRU  £9 

21.  SET  STATE  &  INTEGRATOR  INITIAL  VALUES,  X<0>  &  Z<0> 

22.  SET  INPUT  COMMAND  VECTOR,  V 

23.  ENTER  SAMPLE  TINES 

24.  ENTER  SIMULATION  TIME 

£5.  ENTER  CALCULATION  STEP  SIZE 

26.  RUN  SIMULATION 

27.  SET  CONTROL  INPUT  LIMITS 

28.  OPTION  RESERVED 

29.  COPY  SIMULATION  PARAMETERS  FROM  LOCAL  FILE 
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□rrinn*  please  >  **2i 

TMi:  OPTION  SETS  THE  INITIAL  CONDITION  VECTGPS  FDP  THE  STATES  8.  INTEGRATORS 

CNTEP  THE  X<0>  VECTOR  OF  3  ELEMENTS 
>0*0*0 

ENTER  THE  Z<0>  VECTOR  OF  2  ELEMENTS 
>0*0  • 

OPTION*  PLEASE  >  **22 

THIS  OPTION  SETS  THE  INPUT  COMMAND  VECTOR*  V 

ENTER  THE  V  VECTOR  OF  2  ELEMENTS 
*V"  COLUMN  >3* -l 

OPTION*  PLEASE  >  *>23 

THIS  OPTION  SETS  THE  SAMPLING  TIME  FOP  EACH  RUN 

ENTER  NUMBER  (MAX  OF  2>  OF  SAMPLING  TIMES  >1 
ENTER  1  SAMPLING  TIME(S)  >.02 

OPTION*  PLEASE  >  **24 

THIS  OPTION  SETS  THE  TOTAL  SIMULATION  TIME 
ENTER  TOTAL  TIME  >4 

OPTION,  PLEASE  >  **25 

THIS  OPTION  SETS  THE  CALCULATION  STEP  SIZE 
ENTER  STEP  SIZE  >.02 

SeTIDN*  PLEASE  >  <*2£ 

THIS  OPTION  RUNS  THE  SIMULATION  IF  ALL  DATA  IS  AVAILABLE 

SIMULATION  INCLUDES  NO  ACTUATORS 
SIMULATION  INCLUDES  NO  SENSORS 
SIMULATION  INCLUDES  NO  CONTFOL  LIMITS 

RUN  TIME=4.  STEP  TIME®. 02  SAMPLE  TIME®. 02 

THERE  APE  200  TIME  INCREMENTS 

ENTER  "0"  TO  SKIP  SEQUENTIAL  LISTING 

ENTER  "1"  TO  OBTAIN  THIS  BATA...>0  ! 

- CALCULATIONS  IN  PROGRESS - 

SIMULATION  FOP  SAMPLING  TINE  **1  COMPLETE 
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OpTION»  PLERSE  >  «126 

ENTFP. .. INPUT* OUTPUT. . .FOR  DPT*  MRTRIX  *> OUTPUT 
ENTER  SAMPLING  TIME  OF  INTEREST  RS...1*2»3*  ETC...>1 

OUTPUT  DRTR  MRTRIX... 


TIME 

. 2000E-01 
. GOODE- 01 
. 1 OOOE+OO 
. 1400E+00 
. 1300E+00 
. 2200E+00 
.  2800E+00 
. 3 OOOE+OO 
. 3400E+00 
. 3900E+00 
. 43002+00 
. 4600E+00 
. 5000E+00 
. 5400E+0Q 
. 5800E+00 
. 6200E+00 
. 66 0 0E+ 0  0 
.  7 OOOE+OO 
. 7400E+00 
. 7800E+00 
. 3200E+00 
. 3600E+00 
. 9 OOOE+OO 
. 9400E+00 
.  9800E+G0 
. 1 020E+01 
.  1 060E+01 
. 1100E+01 
.  1 14  0E+01 
.  1 180E+01 
.  1220E+01 
. 1260E+01 
.  13002+01 
.  134  0E+01 
.  1330E+01 
.  1420E+01 
.  1460E+01 
.  1500E+01 
.  1540E+01 
.  1530E+O1 
.  1820E+01 
.  1880E+01 
.  1 700E+01 
.  1740E+01 
.  1730E+01 
.  13203+ 01 
.  1360E+01 


1 

1 277E+00 

3793E+00 

6113E+00 

3236E+00 

1 01SE+01 

1 196E+01 

1  360E+01 

1509E+01 

1646E+01 

1 772E+01 

18868+01 

1991E+01 

20S7E+01 

2175E+01 

2255E+01 

23282+01 

2394 E+ 01 

2455E+01 

2511E+01 

2561E+01 

2607E+01 

2649E+01 

2687E+01 

2721E+01 

2753E+01 

2731 E+ 01 

2307E+01 

2331E+01 

2852E+01 

2S71E+01 

2SSSE+01 

29 04 E+ 01 

2918E+01 

2931E+01 

2942E+01 

2952E+01 

2961E+01 

2969E+01 

2977E+01 

2933E+01 

2989E+01 

2994E+01 

2998E+01 

3002E+01 

3006E+01 

3009E+01 

3011E+01 


2 

421 4E-01 
1229E+00 
1942E+00 
2572E+00 
-.3129E+00 
3622E+00 
-. 4061E+00 
-. 4453E+00 
4302E+00 
-.51 16E+00 
5397E+00 
5651E+00 
-.5331E+00 
6090E+00 
— . 6280E+00 
6454E+G0 
6613E+00 
6760E+00 
-.68952+00 
7021E+00 
-. 7138E+00 
7246E+00 
7348E+00 
7444E+00 
7534E+00 
7619E+00 
— • 7699E+ 00 
7775E+00 
7848E+00 
7917E+00 
7933E+00 
8045E+00 
3106E+0Q 
-.S163E+00 
8219E+00 
3272E+00 
-. 8324E+00 
3373E+00 
8421E+00 
-. 8467E+00 
-.851 1E+00 
-. 8554E+00 
-.S596E+00 
-. 8636E+00 
-. 3675E+00 
-. 3713E+00 
-. 3749E+00 
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19006+01 

. 30146+01 

-.87856+00 

19406+01 

. 3015E+01 

-.88196+00 

19S0E+01 

. 3017E+01 

-. 88536+00 

2020E+01 

.30186+01 

88856+00 

2060E+01 

. 30196+01 

-.89166+00 

2100E+01 

. 30206+01 

-.  8947E+00 

2140E+01 

. 3Q21E+01 

-.89766+00 

21S0E+01 

. 3022E+01 

-.  9005E+00 

2220E+01 

. 30226+01 

-. 9033E+00 

2260E+01 

. 3022E+01 

-.90606+00 

2300E+01 

. 3022E+01 

-.90886+00 

2340E+01 

.3022E+01 

-.91116+00 

2380E+01 

. 3022E+01 

-.9136E+00 

2420E+01 

.30226+01 

-.  9160E+G0 

2460E+01 

.30226+01 

-.91846+00 

2500E+01 

. 30226+01 

-.92076+00 

2540E+01 

. 3021E+01 

-.92296+00 

25806+01 

. 30216+01 

-.  92506+00 

£6 20E+01 

. 3 02 IE +01 

-. 9271E+00 

2660E+01 

. 3020E+01 

-.  9291E+00 

2700E+01 

.30206+01 

-.931 1E+00 

27406+01 

. 30196+01 

-.93316+00 

2780E+01 

. 3019E+01 

-.93496+00 

2820E+01 

. 30186+01 

-.93676+00 

28606+01 

.30186+01 

-.93856+00 

2900E+01 

. 30176+01 

-. 9402E+00 

2940E+01 

.  3  0 1 7E+  0 1 

-.  9418E+00 

£9806+01 

.30166+01 

-.94356+00 

30206+01 

■  30 1 66+L* 1 

-. 9451E+00 

3060E+01 

.  30156+01 

—  • 94666+  0  0 

31006+01 

. 3015E+01 

-. 9481E+00 

3140E+G1 

. 3014E+01 

-.94956+00 

318QE+01 

. 30146+01 

-.95096+00 

3220E+01 

.30136+01 

-.95236+00 

32606+01 

.30136+01 

-. 9537E+00 

33006+01 

. 3012E+01 

-. 9549E+00 

33406+01 

. 301 2E+ 01 

-.95626+00 

33806+01 

.3011 E+ 01 

-. 9574E+00 

34.306+01 

.3011 E+ 01 

-.95866+00 

34606+01 

.30116+01 

— • 95976+00 

35006+01 

.30106+01  • 

-.96096+00 

35406+01 

.30106+01 

-. 96206+00 

35806+01 

. 30096+01 

-.96306+00 

36206+01 

. 30096+01 

— . 964  06+00 

36606+01 

.30096+01 

-.96506+00 

37006+01 

. 30086+01 

-.96606+00 

37406+01 

.30086+01 

-.96696+00 

37306+01 

.30086+01 

-.96796+00 

3820E+01 

.30076+01 

-.96886+00 

38606+01 

. 30076+01 

9696E+00 

39006+01 

. 3007E+01 

9705E+00 

39406+01 

.30066+01 

-.97136+00 

39306+01 

. 30066+01 

-. 9721E+00 
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OPTION*  PLEASE  >  «30 

PLOTTING  OPTIONS* 

30.  LIST  OPTIONS  30  THRU  39 

31.  f'iJICK  SKETCH  AT  USER'S  TERMINAL 

32.  QUICK  SKETCH— SHORT  VERSION 

33.  QUICK  SKETCH— RETAINING  SAME  PLOTTING  CHOICES 

34.  CALCOMP  PLOT 

3*5.  CALCOMP  PLOT— SHORT  VERSION 

?*.  CALCOMP  PLOT— RETAINING  SAME  PLOTTING  CHOICES 

37.  OPTION  RESERVED 

38.  OPTION  RESERVED 

39.  OPTION  RESERVED 


OPTION*  PLEASE  >  «3l 

THIS  OPTION  PRODUCES  A  PLOT  AT  YOUR  TERMINAL 
c'!. EASE  CHOOSE  ONE  OF  THE  FOLLOWING* 

FOP  A  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  UP  TO  £  INPUT  AND  OUTPUT  PAIRS 

2.. .  A  PLOT  OF  UP  TO  4  INPUTS  OR  GUTPIJTS 

3.. . A  PLOT  OF  UP  TO  4  DIFFERENT  FLIGHT  CONDITIONS 

<FQR  ANY  SINGLE  INPUT  OR  OUTPUT) 

OP  FOR  UP  TO  4  DIFFERENT  SAMPLING  TIMES 

4.. . A  PLOT  OF  ANY  SINGLE  INPUT  DR  OUTPUT 

ENTER  CHOICE  DESIRED  >2 

CHOICE  «2... YOU'VE  CHOSEN  TO  PLOT  INPUTS  OR  OUTPUTS 

YOU  MIGHT  HAVE  SELECTED  TO  PUN  SEVERAL  SAMPLING  TIMES... 
DESIGNATE  THESE  SAMPLING  TIMES  AS  ...  1  *  2*3*  ETC... 

ENTER  THE  1  SAMPLING  TIME'S)  OF  INTEREST  >1 

ENTER. .. INPUT  OP  OUTPUT... FOR  YOUR  PLOT  > INPUT 

HOW  MANY  INPUTS  DO  YOU  WANT  TO  PLOT  >2 

Wu I C H  INPUT'S)  DO  YOU  WANT  TO  PLOT... 

ENTER  THE  2  CORRESPONDING  COLUMN  NUMpEP t;S>  >1*2 

POP  NO  GRID  ON  PLOT  ENTER  ”0"*  FOR  A  GRID  ENTER  " 1 M  >1 
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4 
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CURVE  X  ABOVE  IS  INPUT  1 
CURVE  0  ABOVE  IS  INPUT  2 


OPTION?  PLEASE  >  ©31 

THIS  OPTION  PRODUCES  A  PLOT  AT  YOUR  TERMINAL 
PLEASE  CHOOSE  ONE  OF  THE  FOLLOWING: 

P  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  UP  TO  2  INPUT  AND  OUTPUT  PAIRS 

2.. . A  PLOT  OF  UP  TO  4  INPUTS  OR  OUTPUTS 

3 .. .  A  PLOT  OF  LIP  TO  4  DIFFERENT  PLIGHT  CONDITIONS 

(FOR  ANY  SINGLE  INPUT  OR  OUTPUT ) 

OR  FOP  I LIP  TO  4  DIFFERENT  SAMPLING  TIMES 

4.. . A  PLOT  OF  ANY  SINGLE  INPUT  DP  OUTPUT 

ENTER  CHOICE  DESIRED  >2 

CHOICE  «2... YOU'VE  CHOSEN  TO  PLOT  INPUTS  DR  OUTPUTS 

YOU  MIGHT  HAVE  SELECTED  TO  RUN  SEVERAL  SAMPLING  TIMES... 
DESIGNATE  THESE  SAMPLING  TIMES  ETC... 

ENTER  THE  1  SAMPLING  TIME<S>  OF  INTEREST  >1 

ENTER.  ..  INPUT  OP  OUTPUT...  FOR  YOUR  PLOT  >0UTPIJT 
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HQt.i  MANY  OUTPUTS  DO  YOU  WANT  to  PLOT  >2 


WHICH  OUTPUT <S>  I'D  YOU  WANT  TO  PLOT... 

L'UTEP  THE  £  CORRESPONDING  COLUMN  NUMBER <S>  >1»2 


FOR  NO 

GRID  ON 

PLOT 

ENTER 

“0“. 

FOR  A 

GRID  ENTER  “ 

1"  >1 

"»  “i 

-*  •  i-  L_ 

I 

- +- 

> 

iXXXXXXX X XXX X XX 

— +- 

xxxx 

xxxxxx 

XXXXXX  X 

iXXXXXXXXXX 

2.78 

-+ 

+ 

+  xxxxx 

+ 

+ 

+ 

+ 

+ 

+  + 

I 

XXX 

I 

2.  34 

-+ 

+ 

XX  + 

+ 

+ 

+  ' 

♦ 

+ 

+ 

+  + 

I 

X 

I 

1.90 

-+ 

XX 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

I 

XX 

I 

1.46 

-+ 

XX+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

I  X 

I 

1. 03 

-+  X 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

I  X 

I 

.  586 

-+  X 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

IX 

I 

.  146 

-X 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

00 

I 

-.293 

-+  00 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  + 

I 

00000 

I 

-.732 

-+ 

+ 

□00000000000000+ 

+ 

+ 

+ 

+ 

+  + 

I 

OGODOOOOOOOOCO 

cnonooaooocoaooocooooao 

-1.17 

-+ - 

— + — 

- +_. 

- .+. 

— +- 

- +- 

- + — 

0.  3.93 


CURVE  X  ABOVE  IS  OUTPUT  1 
CURVE  0  ABOVE  IS  OUTPUT  2 


OPTION*  PLEASE  >  «:34 

THIS  OPTION  PRODUCES  A  CALCOMP  PLOT 

PLEASE  CHOOSE  ONE  OF  THE  FOLLOWING:  ’  . 

ROc'  A  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  IJP  TO  2  INPUT  AND  OUTPUT  PAIRS 

2. . .A  PLOT  OF  UP  TO  4  INPUTS  DR  OUTPUTS 

3.. . A  PLOT  OF  UP  TO  4  DIFFERENT  FLIGHT  CONDITIONS 

<FOR  ANY  SINGLE  INPUT  OR  DUTPUT> 

OP  FOP  UP  TO  4  DIFFERENT  SAMPLING  TIMES 

4.. . A  PLOT  OF  ANY  SINGLE  INPUT  OR  OUTPUT 

ENTER  CHOICE  DESIRED  >2 

CHOICE  «2. . . YOU'VE  CHOSEN  TO  PLOT  INPUTS  OR  OUTPUTS 
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ynu  night  HAVE  'ELECTED  TO  PUN  SEVERAL  SAMPLING  TINES... 
de* ignatf  these  sampling  times  as...i»2»3»  etc... 

ENTER  THE  1  SAMPLING  TIME (S>  OF  INTEREST  >1 

ENTER. .. INPUT  OR  OUTPUT. ..FOP  YOUR  PLOT  > INPUT 

HOW  MANY  INPUTS  DO  YOU  WANT  TO  PLOT  >2 

WHICH  INPUT <S>  DO  YOU  WANT  TO  PLOT... 

ENTER  THE  2  CORRESPONDING  COLUMN  NUMBER <S>  >1»2 

ENTER  TITLE  FOR  Y-AXIS  OF  PLOT  (MAX  OF  30  CHARACTERS) 

> INPUT" 1  AND  INPUT  «2  < 

TITLE  IS:  INPUT«1  HMD  INPUT  «2 

IS  THIS  CORRECT. . . YES>  NO>  $. . . >Y 

ENTER  MR  IN  PLOT  TITLE  (MAX  60  CHARACTERS) 

>  IRREGUC.BR  PLRNT  DESIGN-  TEST  CASE 
MAIN  TITLE:  IRREGULAR  PLANT  DESIGN-  TEST  CASE 

I S  THIS  CORRECT . . . YES » NO » $. . . >Y 

ENTER  PLOT  SIZE  FACTOR... >.3 


LOCAL  FILE  "PLOT"  CONTAINS  THE  CALCOMP  DATA  FDR  THIS  ENTRY  TO  OPTION  a 
34 

YOU  HAVE  GENERATED  A  TOTAL  OF  1  PLOT (S> 

PE  SURE  TO  ROUTE  "PLOT"  TO  THE  PRINTER  BEFORE  LOGOUT 


OPTION*  PLEASE  >  ©34 

THIS  OPTION  PRODUCES  A  CALCOMP  PLOT 


PLEASE  CHOOSE  ONE  OF  THE  FOLLOWING: 

FOR  A  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  UP  TO  2  INPUT  AND  OUTPUT  PAIRS 

2.. .  A  PLOT  OF  Uc*  TO  4  INPUTS  OR  OUTPUTS 

3.. . A  PLOT  OF  UP  TO  4  DIFFERENT  FLIGHT  CONDITIONS 

(POP  ANY  SINGLE  INPUT  OP  OUTPUT) 

OP  POP  UR  TO  4  DIFFERENT  SAMPLING  TIMES 

4.. . A  PLOT  OF  ANY  SINGLE  INPUT  OR  OUTPUT 

ENTER  CHOICE  DESIRED  >2 

CHOICE  «2... YOU'VE  CHOSEN  TO  PLOT  INPUTS  OR  OUTPUTS 

YOU  MIGHT  HAVE  SELECTED  TO  RUN  SEVERAL  SAMPLING  TIMES... 
DESIGNATE  THESE  SAMPLING  TIMES  AS...1»2>3»  ETC... 

ENTER  THE  1  SAMPLING  TIME(S>  OF  INTEREST  >1 

ENTER. . . INPUT  OP  OUTPUT... FOP  YOUR  PLOT  >OUTPUT 
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Hiju  many  outputs  r<ci  you  mant  to  plot  )2 

MMICH  OUTPUT  OS'*  DO  YOU  MANT  TO  PLOT... 

ENTER  THE  £  CORRESPONDING  COLUMN  NUMBER <S)  )t»2 

ENTER  TITLE  FOR  Y-AXIS  OF  PLOT  <MAX  OF  30  CHARACTERS) 
> OUTPUT?*  1  AND  OUTPUT »2  < 

TITLE  IS*  OUTPUT**  1  AND  DUTPUT«2 

IS  THIS  CORRECT... YES»ND» *. . . >Y 

ENTER  MAIN  PLOT  TITLE  WAX  €>0  CHARACTERS) 

> IRREGULAR  PLANT  DESIGN  -  TEST  CASE 
MAIN  TITLE:  IRREGULAR  PLANT  DESIGN  -  TEST  CASE 

IS  THIS  CORRECT. . . Y£S> ND> $. . . >YES 

ENTER  PLOT  SIZE  FACTOR.  ..>0.9 


LOCAL  PILE  "PLOT"  CONTAINS  THE  CALCDMP  DATA  FOR  T«IS  ENTRY  TO  OPTION 

•34 

YOU  HAVE  GENERATED  A  TOTAL  OF  £  RLOT<S> 

BE  SURE  TO  ROUTE  "PLOT"  TO  THE  PRINTER  BEFORE  LOGOUT 


OPTION.  PLEASE  >  «99 

ALL  PLANT  IN»UT  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 

called  -memo- 

all  DESIGN  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 
CALLED  “ MEM 1 0 " 

all  SIMULATION  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 
CALLED  “MEM20" 


HAVE  A  NICE  DAY! 


STOP 

56600  MAXIMUM  EXECUTION  FL. 

£•3.1!  3  CP  SECONDS  EXECUTION  TIME. 

CCMMAHf'—  REMIND* PLOT 

COMMAND-  ROUTE. PLOT* DC=RT » TID-31 »FID*SMYTH» ST=CSB 


COMMAND-  LOGOUT 
CPA  56.600  SEC. 

10  59.670  SEC. 

CRUS 

CONNECT  TIME  0  H PS. 
l£/Q3/St  LOGGED  OUT  AT 
< 


46.124  ADJ. 
17.662  ADJ. 
74.970 
43  MIN. 
12.29. 13. 


a 


< 
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SECURITY  CLASSIFICATION  OF  This  PAGE  (When  Dai*  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 
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